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Project  MAC  was  organized  at  the  Massachusetts  Institute  of  Technology- 
in  the  spring  of  1963  for  the  purpose  of  conducting  a  research  and  develop¬ 
ment  program  on  Machine-Aided  Cognition  and  Multiple-Access  Computer 
systems.  It  operates  under  contract  with  the  Office  of  Naval  Research, 
acting  on  behalf  of  the  Advanced  Research  Projects  Agency  of  the  Depart¬ 
ment  of  Defense. 

The  broad  goal  of  Project  MAC  is  the  experimental  investigation  of  new 
ways  in  which  on-line  use  of  computers  can  aid  people  in  their  individual 
intellectual  work,  whether  research,  engineering  design,  management, 
or  education.  One  envisions  an  intimate  collaboration  betwQen  man  and 
computer  system  in  the  form  of  a  real-time  dialogue  where  both  parties 
contribute  their  best  capabilities.  Thus,  an  essential  part  of  the  research 
effort  is  the  evolutionary  development  of  a  large,  multiple-access  computer 
system  that  is  easily  and  independently  accessible  to  a  large  number  of 
people,  and  truly  responsive  to  their  individual  needs.  The  MAC  computer 
system  is  a  first  step  in  this  direction  and  is  the  result  of  research  initiated 
several  years  ago  at  the  MIT  Computation  Center. 


Project  MAC  was  organized  in  the  form  of  an  interdepartmental,  inter¬ 
laboratory  "project"  to  encourage  widespread  participation  from  the  MIT 
community.  Such  widespread  participation  is  essential  to  the  broad,  long¬ 
term  project  goals  for  three  main  reasons;  exploring  the  usefulness  of 
on-line  use  of  computers  in  a  variety  of  fields,  providing  a  realistic  com¬ 
munity  of  users  for  evaluating  the  operation  of  the  MAC  computer  system, 
and  encouraging  the  development  of  new  programming  and  other  computer 
techniques  in  an  effort  to  meet  specific  needs. 

Faculty,  research  staff,  and  students  from  fourteen  academic  departments 
and  four  interdepartmental  research  laboratories  are  participating  in 
Project  MAC.  For  reporting  purposes,  they  are  divided  into  sixteen  groups, 
whose  names  correspond  in  many  cases  to  those  of  MIT  schools,  depart¬ 
ments  and  research  laboratories.  Some  of  the  groups  deal  with  research 
topics  that  fall  under  the  heading  of  computer  sciences;  others  with  research 
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topics  which,  while  contributing  in  a  substantive  way  to  the  goals  of  Project 
MAC,  are  primarily  motivated  by  objectives  outside  the  computer  field. 

.  J 

The  purpose  of  this  Progress  Report  is  to  outline  the  broad  spectrum  of 
research  being  carried  out  as  part  of  Project  MAC.  Internal  memoranda 
of  Project  MAC  arc  listed  in  Appendix  A,  and  MAC-rclatcd  theses  are 
listed  in  Appendix  B.  Some  of  the  research  is  cosponsored  by  other 
governmental  and  private  agencies,  and  its  results  arc  described  in 
journal  articles  and  reports  emanating  from  the  various  MIT  depart¬ 
ments  and  laboratories  participating  in  Project  MAC.  Such  publications 
arc  listed  in  Appendix  C  of  the  report.  Project  MAC  Technical  Reports 
are  listed  in  Appendix  D. 


Robert  M,  Fano 
Cambridge,  Massachusetts 
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A  Visually-Controlled  Manipulator  -  Marvin  L.  Minsky 

Since  the  completion  of  a  computer-controlled  mechanical  hand  in  1961 
(H.  Ernst,  Ph.  D.  Thesis,  M.I.T.  ,  l  y6l)  little  has  been  done  in  the  area 
of  design  of  autonomous  manipulators.  We  have  been  working  the  develop¬ 
ment  of  a  versatile  ,  computer-operated,  visually-oriented  machine  for 
handling  objects  in  complicated  spatial  situations.  A  real-time,  on-line, 
television-camera  interface  has  been  built  for  the  PDP-6  computer  and 
has  been  used  to  track  the  motion  of  simple  visual  objects,  including  the 
human  eye.  This  is  of  interest  in  studying  the  mechanism  of  human 
visual  preception.  Also,  an  AMF  Versatran,  Model  C,  mechanical 
remote  manipulator  has  been  coupled  to  the  PDP-6. 

Unrecognizable  Sets  of  Numbers  -  Marvin  L.  Minsky  and  Seymour  Papert 

The  problem  here  is,  given  a  set  A  of  non-negative  integers,  is  there 
a  finite-state  device  that  will  discriminate  between  members  and  non¬ 
members  of  A  when  the  numbers  are  presented  in  binary,  or  some 
other  similar  numbe r- representation  system  with  a  different  radix  (such 
as  ternary  of  decimal)?  A  technique  for  showing  that  certain  sets  are 
not  recognizable  has  been  written.  The  technique  deals  with  the  asymptotic 
behavior  of  the  function  ir^  (n),  which  is  the  number  of  integers  less  than 
n  in  the  set  A;  unless  tt^  (n)  has  a  specified  behavior,  A  is  not  recogniza¬ 
ble.  However,  not  all  sets  A  with  a  satisfactory  asymptotic  behavior  of 
-A  («)  recognizable.  Another  technique  developed  by  McNaughton 
handles  some  of  these  cases.  The  goal  for  research  in  this  area  is  to 
determine  a  necessary  and  sufficient  condition  that  A  be  recognizable 
in  r-ary  (r  2  2)  by  a  finite  state  device.  A  paper  will  appear  soon  in  the 
Journal  of  the  ACM. 
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The  Algebraic  Approach  to  Finite  Automata  -  Seymour  Paper!  and 
Robert  NcNaughton 

The  work  by  Schutzenberger,  Rhodes  and  Krohn,  and  some  of  our  own 
results  have  convinced  us  of  the  fruitfulness  of  the  algebraic  approach 
to  finite  automata,  and,  in  particular,  the  concept  of  the  semi-group  of 
the  machine.  The  two  problems  that  have  been  most  impressively  dealt 
with  are  the  problem  of  decomposition  of  machines,  and  the  problem  of 
general  star  height.  It  is  our  intention  to  continue  to  focus  on  this  area. 
Research,  however,  will  lean  toward  application  of  already  existing 
knowledge  to  problems  about  machines,  rather  than  advancing  knowledge 
in  the  field  of  algebra.  We  have  been  concerned  with  simplifying  some 
of  the  proofs  of  known  theorems  (applying  algebra  to  machines)  both  to 
make  fundamental  notions  clear  enough  for  future  research  and  thinking, 
and  also  to  bring  these  concepts  within  the  purview  and  jjarlance  of  begin¬ 
ning  students  in  the  field. 

It  is  felt  that  use  of  the  computer  will  be  fruitful  in  this  area.  Some  of 
the  programmers  in  the  Artificial  Intelligence  Group  have  been  writing 
LISP  programs  that  will  do  such  things  as  find  the  semigroup  of  a  given 
finite-state  machine.  Use  of  a  time-shared  computer  in  this  way  can  be 
v'ery  useful  in  checking  out  hunches  in  the  search  for  valid  generalities 
or  in  providing  counter-examples. 
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The  Design  of  LISP  2  -  Daniel  G.  Bobrow 

The  design  of  a  new  programming  language,  called  LISP  2,  is  now  being 
cariied  out  jointly  by  Systems  Development  Corporation  and  Information 
International  Incorporated,  with  aid  from  members  of  the  Project  MAC 
Aitificial  Intelligence  Group.  LISP  2  will  be  an  ALGOL-style  language  and 
will  have  the  desirable  ALGOL  features  of  dynamic  allocation  of  storage, 

1  ecui  sive  functions,  and  dynamic  arrays.  LISP  2  will  have  list  process¬ 
ing  as  in  the  current  LISP  1,  5  language;  however,  the  control  language 
will  be  in  ALGOL  rather  than  in  LISP.  In  this  list  processing,  useful 
list  structures  will  be  compacted  to  provide  room  for  the  allocation  of 
space  for  arrays  aitd  for  new  list-structure  elements. 

In  addition  to  the  LISP  1.5  function-oriented  list  processing,  LISP  2  will 
embody  format-directed  string  processing  similar  to  that  of  the  COMIT 
programming  language.  With  this  processing  a  user  can  specify  a  pattern 
against  which  a  string  is  matched.  If  it  is  matched  successfully,  the 
stiing  is  parsed,  that  is,  subdi\ided  into  segments  which  match, the 
individual  patterns  in  the  given  pattern.  These  segments  can  then  be 
used  to  construct  a  new  segment  according  to  a  given  formal.  This 
processing  is  useful  in  the  manipulation  of  natural  language  strings. 

Since  the.  front  end  of  the  LISP  2  system  will  contain  a  finite-state  machine 
as  an  input-string  grouper  and  a  syntax-directed  compiler,  the  user  will 
be  able  to  change  the  language  easily.,  specify  new  syntactic  forms,  define 
new  data  types,  and  make  semantic  additions  to  the  languag-e-.-  -  - 

The  compilei  for  LISP  2,  first  compiled  w'ith  a  bootstrap  compiler  written 
in  LISP  1.  5,  w'ill  be  written  almost  completely  in  LISP  2  itself;  LISP  2 
will  thus  be  relatively  independent  of  machine  types,  and  we  expect  to 
create  versions  of  it  on  a  number  of  machines,  such  as  the  Q-32  at  SDC, 
the  IBM  360,  the  DEC  PDP-6,  and  the  GE  645.  The  compiler  will  w'ork 
in  a  two-stage  process.  First,  the  external  communication  language  will 
be  translated  into  an  intermediate  language,  which  is  compatible  upward 
from  the  current  LISP  1.  5  notation.  Then  a  list -structure  assembly- 
language  program  will  be  generated  and  compiled  by  a  simple  machine - 
code  assembly  program.  A  preliminary  version  of  the  LISP  2  system 
should  soon  be  in  operation  on  the  Q-32  at  SDC. 
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Man-Compul.er  Symbiosis  -  Warren  Teilelman 

A  LISP  programming  system  is  being  constructed  to  facilitate  the  develop¬ 
ment  of  sophisticated  problem-solving  systems.  Emphasis  is  on  easing 
the  familiar  debugging  sequence:  write  some  code,  run  the  program, 
make  some  changes,  write  more  code,  run  the  program  again,  etc.  As 
a  system  becomes  more  complex,  making  changes  becomes  harder  and 
harder. 

Our  goal  is  to  make  the  computer,  via  a  programming  system,  play  an 
active  role  in  this  modification  process  by  providing  the  means  whereby 
changes  can  be  effected  immediately,  in  ways  tliat  seem  natural  to  the 
user.  The  flavor  of  the  system  is  such  lliat  the  user  feels  he  is  giving 
advice,  or  making  suggestions  to  the  computer  about  the  operations  of 
his  programs;  the  necessary  work  to  achieve  llie  desired  effect  is  per¬ 
formed  by  tiie  symbiotic  system.  The  system  thus  acts  as  an  interface 
between  the  user  and  his  programs. 

A  "format  list-processing  language",  FLIP,  has  been  developed  and 
implemented.  FLIP  performs  operations  of  the  type  done  by  COMIT, 
but  is  considerably  more  sopliisticated  and  operates  from  within  a  LISP 
system.  This  language  will  be  used  to  process  the  advice  given  by  the 
user  and  will  ena'Te  the  symbiotic  system  to  accept  suggestions  in  a 
very  relaxed  format.  It  has  already  been  used  to  construct  a  set  of 
editing  functions  for  LISP  programs. 

Consistent  with  emphasis  on  the  importance  of  change  and  flexibility,  the 
system  itself  may  easily  be  modified  by  giving  it  advice  about  its  own 
operation. 
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Mathematical  Assistant  -  William  A.  Martin 

I  have  embedded  a  picture-display  language  in  the  LISP  language  of  the 
PDP-6.  Now  complex  display  programs,  such  as  the  ARGUS  character- 
recognition  scheme  of  W .  Teitelman,  are  easily  written.  I  have  also 
developed  a  LISP  program  that  displays  mathematical  expressions  and 
equations.  At  present,  the  displayed  equations  have  substantially  the 
same  form  as  those  appearing  in  mathematical  publications.  This 
mathematical  display  program  in' PDP-6  LISP  communicates  with  the 
rest  of  the  mathematical  assistant  programs  in  the  time-shared  MAC 
LISP  through  an  input/ output  dataphone  system. 

The  mathematical  assistant  system  is  a  new  technique  for  the  comparison 
and  simplification  of  symbolic  mathematical  expressions.  It  uses  finite- 
field  arithmetic  for  hash-coding  functions  of  a  complex  variable.  By  con¬ 
structing  a  large  but  finite  field  (p=8,  589,  949,  373)  within  which  there  are 
representations  of  the  integers  and  tt,  in  an  appropriate  sense,  one  can 
use  exponent  arithmetic  to  map  symbolic  expressions  onto  integers.  If 
two  expressions  are  equivalent  with  respect  to  rational,  exponential,  and 
some  trigonometric  identities,  they  will,  with  high  probability,  map  onto 
the  same  integer.  In  the  recognition  of  common  subexpressions,  this 
system  appears  to  be  faster  and  more  accurate  than  those  which  use 
canonical  orderings. 

Heuristics  of  Theorem  Proving  in  Group  Theory  -  Lewis  M.  Norton 

We  are  developing  a  system  of  heuristics  of  theorem  proving  in  group 
theory.  It  is  coded  in  LISP  and  handles  problems  that  involve  the  simple 
consequences  of  definitions.  We  are  now  using  it  to  test  methods  of  re¬ 
ducing  irrelevant  effort  in  proofs  already  produced.  Rather  than  use 
explicit  foims  of  the  predicate  calculus,  it  uses  forms  more  similar  to 
stands  id  mathematical  practice.  Thus,  the  definition  of  abelian  for  a 
group  G  is  written: 

((ABELIAN  G  G)  IMPLIES  (AND  (MEMBER  AI  G)  (MEMBER  A2  G)) 
(EQUAL(=:=PROD  Al  A2  G)  (-'PROD  A2  Al  G))). 

Since  logical  axioms  are  embedded  in  the  flow  of  the  program,  the 
involved  theorems  can  be  proven  more  quickly  than  has  been  done  by 
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other  methods.  For  e.xample,  the  system  proves: 

Gi\'en:  G  a  group,  J  a  subgroup  of  G,  M  a  normal 
subgroup  of  G. 

Prove:  .7  P  II  i.s  normal  in  J. 

In  this  work,  we  are  trying  to  simulate  ihe  analytical  and  rational 
methods  of  human  problem  solv'ing. 

Computer  Expe  rimeitl  s  iit  Finite  Algebra  -  \V.  Douglas  .Maurer 

Any  set  of  mathematical  methods  and  iechni(|ues  may  be  mec’uinir.ed  on 
a  digital  computer.  Even  infinite  mathematical  structures  may  be 
studied,  provided  that  the  e.xprcssions  used  to  speak  about  them  aic 
finite,  since  it  is  these  symbols  which  the  computer  manipulates.  Ihe 
first  mechanization  is  for  finite  algebra.  This  system  for  Unite  algebr.-i 
consists  of  cin  English-like  language  with  80  coirrmarids  Ji.>r  the  ntanipu- 
lation  of  finite  groups  and  semigroups,  sitbsets,  maps,  and  labeled 
constants.  Rings  and  fields  are  to  be  added  later. 

The  system  accepts  an  arbitrary  semigroup  from  the  console,  il  the 
subgroup  is  gi\'en  by  its  Cayley  table  {multiplication  table).  It  can  do 
the  following: 

1.  construct  the  Cayley  tables  of  lour  special  kinds  of 
group  and  three  special  kinds  ol  semigroup: 

2.  take  the  direct  product  of  any  two  Cayley  tables  and 
produce  a  third; 

5.  generate  the  subsemigroup  of  any  semigroup  generated 
by  a  single  element  or  subset  of  that  semigroup  and  then  proceed  to 
build  the  Cayley  table  of  the  subsemigroup; 

4.  find  all  the  subsemigroups  of  a  semigroup  and  all  the 
normal  and  ma.ximal  subgroups  of  a  group; 

5.  find  all  the  left,  right,  two-sided,  or  only  the  maximal 

ideals  of  a  semigroup: 

6.  test  whether  a  Cayley  table  is  associative  and  whether 
a  semigroup  has  an  identity,  is  a  group,  and  is  abelian; 

7.  tell  whether  a  map  from  one  semigroup  is  homomorphic 
to,  one-to-one,  or  onto  another  semigroup; 
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8.  prim  out,  erase  from  its  storage  area,  or  rename  any 
Cayley  table,  subset,  map,  and  labeled  constant; 

9.  add  a  zero  element  and  a  unit  element  to  a  semigroup; 

10.  produce  the  natural  map  of  a  group  onto  the  set  of  left, 
or  right,  cosets  modulo  a  subgroup,  or  the  Cayley  table  of  the  factor 
group  modulo  a  subgroup. 

I  have  recently  added  a  fuirction  to  the  system  which  determines  whether 
any  two  given  Cayley  tables  are  isomorphic  and,  if  so,  demonstrates  the 
isomorphism.  This  function  runs  in  less  than  thirty  seconds  for  two 
groups  of  twelftit  order.  I  plan  to  introduce  a  program  which  accepts 
mathematical  statements  in  a  source  language  and  produces  object  codes 
to  verify  or  disprove  the  statemeiits  over  a  restricted  class  of  test  objects, 
such  as  a  collection  of  counterexamples. 

A  Chess -Playing  Program  -  Burton  H.  Bloom 

,An  elaborate  static-position  evaluator  has  been  substantially  debugged 
for  the  CHESS  program.  It  includes  material  balance,  development, 
center-control,  and  space  e\'aluations.  Current  plans  are  to  also  in¬ 
clude  king-safety  and  pawn-weakness  and  some  other  criteria.  Several 
plausible -move  generators  are  being  tested.  These  look  for  captures 
for  gain,  checks,  threads  to  capture  for  gain,  retreating  attacked  units, 
and  safe  moves  that  avoid  immediate  loss  of  material.  Plans  are  under¬ 
way  to  provide  for  liandling  defensive  moves,  interposition,  development, 
and  center  control. 

Strategic  heuristics  arc  being  debugged  for  complicating  situations  to 
exploit  possible  opponent's  errors,  and  for  controlling  the  use  of  the 
different  generators  in  different  situations.  An  elaborate  system  of 
Macro-operations  is  used  for  defining  the  evaluators  and  generators. 

With  debugging  of  current  features,  the  program  is  expected  to  match 

low-to-averagc  club  players.  ,  ^ 
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MATHLAB:  On-Line  Symbolic  Computation  -  Carl  Engleman  (MITRE 
Corporation) 

MATHLAB  is  a  LISP  program  that  is  intended  to  provide  the  services 
of  an  on-line,  time-shared  computer  to  a  mathematical  scientist,  such 
as  a  physicist,  applied  mathematician,  or  electronics  engineer,  for  the 
daily  performance  of  his  tedious,  mechanical,  mathematical  tasks. 

In  its  present  state  of  development,  the  program  is  primarily  oriented 
towards  symbolic,  rather  than  numerical,  computation.  Employing 
simple  commands,  such  as  "solve"  or  "differentiate",  the  program  has 
facilities  for  symbolic  manipulations,  such  as  addition,  multiplication, 
and  differention  of  e.xpressions  and  equations.  It  is  also  capable  of  the 
simplification  of  e.xpressions  and  the  substitution  of  one  e.xpression 
within  another. 

Finally,  it  can,  using  the  INTEGRATE  program  of  M.  Manove,  succeed 
in  the  symbolic  integration  of  certain  e.xpressions  and  the  symbolic  solu¬ 
tion  of  certain  e.xpressions. 

A  detailed  discussion  of  the  MATHLAB  program  may  be  found  in  MITRE 
Corporation  Technical  Memorandum  TM-04258,  from  Department  73, 
dated  26  May  1965. 
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INTEGRATE;  On-Line  Indefinite  Integration  -  Michael  Manove 
(MITRE  Corporation) 

INTEGRATE  is  a  LISP  program  designed  to  symbolically  integrate 
rational  functions  in  one  variable  over  the  field  of  rational  numbers. 

Any  rational  function  is  admissible  as  input,  and  such  functions  may 
be  represented  as  the  sums,  differences,  products,  quotients,  or 
integral  powers  of  constants,  polynomials,  and  other  rational  functions. 
The  antiderivative  of  a  rational  function  will  be  expressed  as  the  sum 
of  a  rational  function  with  rational  coefficients,  and  the  logarithm  of 
a  rational  function  possibly  with  complex  coefficients. 

The  program  will  also  integrate  many  rational  functions  of  trigonometric 
functions,  of  the  form  f(sin  x,  co.x  x,  tan  x,  cot  x,  sec  x,  esc  x)  where 
"f"  is  a  rational  function,  but  this  feature  is  limited  by  the  fact  that 
answers  are  always  expressed  in  terms  of  "tan  (l/2)x". 

The  coding  of  INTEGRATE  is  divided  into  tw'O  major  sections.  The 
first  section  will  find  the  rational  part  of  the  antiderivative  and  determine 
rational  functions  whose  integrals  are  the  summands  of  the  logarithmic 
part.  The  second  section  will  find  the  logarithmic  parts  of  the  anti- 
derivatives  of  a  broad  range  of  rational  functions.  This,  of  course,  is 
subject  to  the  size  limitations  of  computer  memory  and  other  hardware 
variables. 

This  program  has  been  incorporated  in  the  MATHLAB  program  of 
C.  Engleman.  A  detailed  discussion  of  INTEGRATE  may  be  found  in 
MITRE  Corporation  Technical  Memorandum  TM-04204,  from  Depart¬ 
ment  73,  dated  22  April  1965. 
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Molecular  Model  Buildlim  -  Cyrus  Levinthal 

Many  of  the  recent  developments  in  molecular  biology  have  depended  on 
an  understanding  of  the  three-dimensional  structure  of  large  molecules. 
Such  an  understanding  is  usually  reflected  in  the  ability  to  construct  a 
three-dimensional  model  of  a  particular  molecule.  For  really  large 
molecules,  however,  such  construction  is  difficult  and  time  consuming 
and,  if  complete,  would  entail  an  impossible  enumeration  of  the  many 
small  interactions  which  contribute  to  the  molecular  stability.  Thus,  if 
a  structure  is  known,  physical  models  can  be  built,  but  construction 
complexity  prevents  the  use  of  such  models  in  the  examination  of  a  large 
numbter  of  different  configurations. 

During  the  last  si.x  months  we  have  written  a  set  of  programs  with  which 
we  can  construct,  display,  and  analyse  macromolecules  on  the  MAC 
system's  ESL  display  unit.  With  computer-controlled  display  and  real¬ 
time  control  of  the  molecule,  an  observer  can  obtain  a  true  three-dimen¬ 
sional  visulization  of  a  particular  molecule:  and  by  interacting  with  the 
program  which  generates  the  molecule,  he  can  indicate  the  way  in  which 
particular  parts  of  the  structure  are  to  be  altered. 

The  first  program  written  in  this  project  calculates  the  coordinates  of  the 
atoms  in  a  protein.  Only  those  angles  about  which  rotation  is  possible 
are  entered  as  input  variables  in  the  program'  all  other  rotation  angles 
and  chemical  bond  lengths  are  entered  as  rigid  constraints.  The  calcula¬ 
tion  treats  each  chemical  bond  along  the  linear  peptide  backbone  as  a 
translation  and  a  rotation  of  a  coordinate  system  attached  to  individual 
atoms.  In  this  way,  the  step  from  one  bacltbone  atom  to  the  next  involves 
a  translation  and  a  rotation-matrix  multiplication.  The  updating  of  the 
rotation  matrix  is  determined  either  by  the  fixed  rotation  angles,  which 
are  introduced  as  constraints  in  the  program,  or  as  input  data  which  are 
added  by  an  investigator.  The  basic  data  for  these  constraints  on  bond 
angles  and  lengths  have  been  obtained  over  a  number  of  years  from  the 
X-ray  crystallographic  studies  of  Pauling,  Perutz,  Kendrew,  and  many 
others,  and  they  are  now  firmly  based  on  experimental  information. 
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The  th  1  e e - di mens ioiial  structure  of  the  protein  is  determined  entirely  by 
the  physico-chemical  interactions  which  lead  to  the  lowest  energy  con- 
figuiation.  Since  the  bulk  of  the  free  energy  is  contributed  by  many  small 
inte  1  ac tion s I  which  for  the  most  part  have  extremely  short  range,  we 
should  like  to  determine  for  each  pair  of  atoms  whether  the  two  members 
lie  close  together  in  real  space  and  thus  contribute  substantially  to  the 
interaction  energy.  Since  even  a  small  protein  molecule  may  contain  1500 
atoms,  the  total  number  of  pairs  to  be  tested  is  very  large. 

I 

To  a\oid  the  enumeration  of  all  possible  pairs,  w^e  have  wuutten  a  set  of 
programs  that  first  divides  space  into  cubes  and  then  makes  a  list  of  all 
atoms  in  each  cube.  By  searching  through  the  cubes,  and  listing  those 
pairs  in  which  one  mentber  is  in  a  cube  and  a  second  member  is  in  the 
same  cube  or  in  one  of  the  twenty-six  surrounding  ones,  we  are  able  to 
enumeiate  all  pairs  of  tteighboring  atoms.  This  list-processing  procedure 
finds  the  holes  in  a  molecule  by  locating  empty  cubes  surrounded  by  filled 
ones;  and,  in  addition,  it  determines  the  inside  and  outside  of  the  molecule 
by  finding  cubes  which  have  filled  neighbors  on  one  side  and  empty  neigh¬ 
bors  on  the  other  side.  The  knowledge  of  these  topological  aspects  of  a 
folded  protein  is  necessary  for  the  assessment  of  interactions  of  particular 
amino  acid  residues  with  the  water  normally  surrounding  all  protein  mole¬ 
cules.  These  programs  used  SLIP  subroutines  written  by  Professor 
Weizenbaum,  wdiose  cooperation  W'as  invaluable  in  our  work. 

When  an  investigator  wishes  to  alter  the  structure  of  a  molecule,  he  inserts 
a  ps eudo - ene  1  gy  between  a  particular  pair  of  atoms.  An  energy-minimizing 
routine,  which  looks  at  the  sum  of  all  interactions,  including  the  pseudo- 
eneigy,  calculates  the  appropriate  change  required  for  the  variable  angles. 
If  the  change  leads  to  a  violation  of  van  der  Waals  contacts,  another  routine 
eliminates  thi°  violation.  The  latter  routine,  as  well  as  those  for  energy 
minimization,  require  the  calculation  of  the  partial  derivatives  of  the  dis¬ 
tances  between  pairs  of  atoms  with  respect  to  all  of  the  variable  angles 
which  can  effect  these  distances.  We  have  used  tw'o  procedures  in  altering 
the  structure  and  eliminating  van  der  Waals  violation:  the  first  uses  one 
of  several  different  minimization  routines  in  a  sequence  of  small  steps;  the 
second  uses  a  linear  approximation  and  calculates  appropriate  angular 
change  by  solving  linear  equations,  subject  to  the  additional  constraint 
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that  the  sum  of  all  angle  changes  shoiilcl  be  a  minimum.  When  we  have 
tested  and  run  programs  and  subroutines  together  in  a  smooth  system,  we 
can  attempt  to  solve  for  the  unknown  structure  of  a  real  protein. 

By  making  use  of  molecular  models  constructed  in  this  way,  we  can 
evaluate  the  importance  of  electrostatic  interactions  in  the  stability  of 
protein  configurations.  These  interactions  occur  because  of  the  e.xistence 
of  a  large  number  of  small  electric  dipoles.  Although  each  of  the  dipole 
interactions  is  small,  they  have  a  longer  range  of  interaction  than  the 
\-an  del-  Waal  forces  and  arc  so  numerous  that  they  contriliutc  appreciably 
to  overall  molecular  stability.  The  results  of  our  calculations  indicate 
that  dipole  interactions  produce  a  stabilizing  energy,  if  the  helical  regions 
of  the  proteins  are  antiparallel  to  each  other,  and  a  destabilizing  energy, 
if  the  helices  are  parallel  to  each  other. 

The  structure  of  the  myoglobin  protein  molecule  has  been  completely 
determined  by  X-ray  crystallographic  methods.  When  we  used  the  appro.\i- 
matc  coordinates  of  this  molecule  in  our  model,  we  found  the  net  stabilizing 
effect  of  these  dipole -dipol e  interactions  to  be  significant.  W'e  have  also 
used  our  model -builditig  program  to  refine  coordinates  of  the  myoglobin 
model,  by  supplemeitting  the  X-ra\-  data  with  data  on  the  chemical  bonds 
derived  from  otiter  sources.  Prior  to  tliis  computer  refinement,  our 
attempts  to  carry  out  a  calculation  of  this  kind  iit’.'oK’ed  the  construction 
and  measurement  of  physical  models  made  of  brass  wires.  This  calcula¬ 
tion  was  both  tedious  and  inaccurate  for  a  molecule  as  latcjje^^a^'myoglobin. 

W'e  have  written  a  program  that  uses  Project  MAC'S  PDP-6  as  a  display 
unit  for  structures  whose  coordinates  are  generated  on  the  709d.  'Ihis 
program  was  designed  to  test  the  usefulness,  as  aids  in  the  visualization 
of  a  three-dimensional  structure,  of  pefspectit’c  and  modulation  of  line 
intensity  with  depth.  It  is  written  in  such  a  way  that  one  cait  test  the 
feasibility  of  a  three-dimensional  visualization  scheme  that  requires  fewer 
rotation  calculations  than  the  ESL  display.  These  tests  will  determine  the 
required  characteristics  of  additional  off-line  display  equipment. 

^V c  have  written  a  set  of  programs  which  calculate  the  strain  energy  of 
displacement  of  chemical  bond  angles  and  lengths  from  their  equilibrium 
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positions.  When  the  rigi'd' constraints  on  bond  length  and  angles  in  the 
programs  described  previously  are  relaxed,  the  number  of  interactions 
in  a  protein  becomes  prohibitively  large.  However,  studies  of  small- 
molecule  structure  can  be  carried  out  in  a  reasonable  amount  of  computer 
time.  We  have  used  this  procedure  to  refine  models  of  the  sugars  which 
are  a  part  of  nucleic  acid  molecules,  and  we  are  studying  the  effect  of 
this  on  the  overall  structure  of  nucleic  acids.  Since  we  can  cany  out 
this  refinement  for  any  small  organic  molecule,  we  can,  in  principle, 
display  the  three-dimensional  model  of  a  small  molecule  when  our  only 

knowledge  of  it  is  what  atoms  compose  it  and  which  ones  are  covalently 
bonded  together. 

[Editor's  Note:  Figures  1  and  Z  are  stereographs  of  a  DNA  molecule  and 
several  protein  molecules,  that  were  photographed  from  the  ESL  Display 
Console  at  Project  MAC.  To  ,use  the  stereo  pairs  without  a  viewer, 
ignore  the  double  image  print  for  the  moment  and  focus  your  eyes  on  a 
point  about  10  feet  away.  Bring  the  figure  slowly  up  into  your  line  of 
sight,  maintaining  the  same  focus  of  10  feet.  Four  images  will  be  within 
your  field  of  view,  and  then  the  two  middle  images  should  merge  into  a 
central  image.  You  now  have  three  images  side  by  side.  The  center  one 
will  be  the  stereo  composite  of  the  two  printed  images.  You  may  need 
practice  to  hold  this  image  until  your  eyes  focus  on  the  stereo  image. 

If  you  have  difficulty  with  the  multiple  images,  try  holding  your  hand, 
edge  on,  in  front  of  your  nose,  so  that  the  left  image  cannot  be  seen  by 
your  right  eye,  and  the  right  image  cannot  be  seen  by  your  left  eye.  This 
will  form  a  crude  stereoscope  and  allow  you  to  concentrate  on  a  single 
image,  the  stereo  composite. 


Figure  1  shows  three  views  of  the  same  portion  of  a  DNA  molecule.  ' 

These  were  produced  by  photographing  stationary  images  that  differ  by 
approximately  10  degrees  of  rotation  about  the  Z  axis.  Normally,  a 
person  at  the  ESL  Display  Console  is  able  to  obtain  a  solid-appearing  view 
of  the  molecule  by  real-time  rotation  of  the  image  about  three  axes.  The 
views  shown  in  Figure  1  are  approximately  10  million  times  life  size. 
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Deoxyribonucleic  Acid  (DNA)  is  the  carrier  of  genetic  information  in 
living  systems,  and  consists  of  two  phosphate  chains  that  run  in  opposite 
directions  and  are  coiled  into  a  right-hand  helix.  Nitrogen  bases,  held 
together  by  hydrogen  bonds,  link  the  ester  chains  like  steps  of  a  spiral 
staircase  and  provide  the  sequence  which  carries  the  genetic  coding. 

Figure  2  shows  three  views  of  complete  and  abbreviated  polypeptide 
protein  chains  arranged  in  the  a  -  helix  configuration.  The  a  helix,  capable 
of  forming  hydrogen  bonds  between  all  of  its  amide  groups,  is  a  substantial 
constituent  of  protein  molecules,  w'hich  perform  vital  enzymic  functions, 
such  as  controlling  chemical  reactions,  in  living  systems. 

The  top  view'  in  Figure  2  is  a  fairly  complete  polypeptide  structure, 
while  the  other  two  views  are  simplified  versions  of  the  same  molecule. 
The  bottom  image  is  the  same  molecule  as  the  middle  one,  showing  the 
e.xaggerated  stereo  separation  obtained  with  approximately  20  degrees  of 
rotation.] 
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Computer-Aided  Teaching  -  Daniel  Roos 

The  Department  of  Civil  Engineering  has  used  the  Project  MAC  time-sharing 
system,  in  its  undergraduate  and  graduate  education  programs.  Three 
different  types  of  educational  activities  relied  heavily  upon  CTSS:  class¬ 
room  demonstrations,  homework  and  term  paper  assignments,  and  an 
experimental  time-sharing  section  of  course  1.  15. 

The  time-sharing  system  has  been  used  very  successfully  in  classroom 
demonstrations  to  assist  in  the  teaching  of  engineering.  In  a  classroom 
equipped  with  closed-circuit  television,  an  IBM  1620  computer,  CalComp 
and  Gerber  plotters,  and  an  IBM  1050  remote  console  connected  to  the 
Project  MAC  computer,  the  instructor  has  all  the  computing  facilities 
available  to  demonstrate  points  that  he  is  making.  If  a  student  asks  the 
perennial  question,  ''What  happens  if  ?"  the  teacher  can  immediately 

demonstrate  the  implications  using  on-line  computer  facilities. 

Several  problem-oriented  languages  have  been  developed,  with  which  an 
engineer  can  communicate  with  the  computer  in  engineering  terminology. 

Two  of  these  languages,  STRESS  for  structural  problems  and  COGO  for 
geometric  problems  are  being  used  in  conjunction  with  time-sharing  for 
educational  applications  by  both  instructors  and  students.  A  student  does 
not  need  intimate  computer  knowledge  to  effectively  use  these  languages. 

With  minimum  instruction,  generally  one  hour  or  less,  he  can  use  the 
problem-oriented  language  and  the  time-sharing  system  to  obtain  solutions 
to  engineering  problems. 


i 


I 
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.  The  Structural  Design  Language  -  Robert  D.  Logcher,  Gerald  M.  Sturinan, 

Richard  V.  Goodman,  Lewis  C.  Teague,  Salvatore  G.  Mazzotta,  and 
John  R.  Roy- 

Structural  design,  the  determination  ol  structural  elements  and  their  lay'out 
and  makeup,  is  a  task  that  involves  innumerable  variables  and  decision 
parameters  and  requires  extensive  computation  and  decision-making.  The 
nature  of  the  computations,  and  the  manner  in  which  they  are  executed, 
makes  classical  forms  of  computer  use  inadequate  and  unsatisfactory-.  The 
purpose  of  the  work  described  here  is  to  determine  the  lorm  and  organiza¬ 
tion  required  for  a  useful  design  system. 

Structural  design  requires  interactions  with  other  technologies  involved  in 
the  design  product,  such  as  mechanical  and  electrical  engineering  and 
architecture.  The  engineer  must  accept  and  recognize  the  conflicting  aims 
of  these  technologies  and  evaluate  compromises  analytically.  He  is  there¬ 
fore  engaged  in  an  iterative  process  in  which  he  must  remain  tlexible.  The 
iteration  consists  usually  of  a  series  of  trial  solutions  together  with  ait 
intuitive  searching  process  that  leads  to  an  optimum  solution. 

Generally  the  amount  of  computation  performed  between  decisions  is  small, 
but  varies  greatly  with  the  problem,  the  engineer,  and  the  various  parts  ol 
the  design  task.  Problems  vary  widely  in  size  and  complexity,  from  a  few 
hundred  variables  and  parameters  to  thousands.  Depending  on  the  problem 
'  and  the  engineer's  sophistication,  the  design  task  may  be  automated  with 
nil  decisions  made  prior  to  solution  or  may  be  broken  into  many  processes 
with  decisions  inserted  during  solution.  A  process  ot  design  such  as 
elastic  analysis  involves  a  great  deal  of  computation;  while  others,  such  as 
checking  performance  against  criteria,  require  very  little. 

A  ti  me  -  sha  r  ing  environment  is  required  to  join  the  creative  design  process 
to  a  computer.  The  user  may  thus  perform  small  and  large  tasks  in  an 
unpredetermined  order  and  e.xpress  decisions  between  tasks.  I  he  problem- 
oriented  language  STRES.S  represents  the  first  phase  in  computer-aided 
'  structural  design.  This  language  contains  a  data-input  mechanism  and 

programs  for  linear  elastic -frame  analysis.  .Since  this  analysis  is  the  most 
difficult  computational  problem,  but  the  best  delined  mathematically,  its 
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solution  was  the  first  undertaken.  General  data- structure  capabilities  for 
engineering  data  were  developed  for  STRESS  which  make  its  logical  inclu¬ 
sion  in  a  broader  language  possible. 

Over  the  past  year,  the  Structural  Design  Language  (STRUDL)  has  been 
developed  in  an  open-end  form.  The  major  commands  of  this  system  are; 


1. 

PROBLEM  INPUT 

6. 

SELECTIVE  OUTPUT 

2. 

PRELIMINARY  ANALYSIS 

7. 

WEIGHT  DETERMINATION 

3. 

DETERMINATE  ANALYSIS 

8. 

CHECK  PROCESS 

4. 

MEMBER  SELECTION 

9. 

SPECIFICATION  CHECK 

5. 

STIFFNESS  ANALYSIS 

10. 

REVISE  TABLES 

Each  of  these  calls  a  particular  operation  commonly  used  in  the  design 
piocess.  The  engineer  can  decide  which  operations  are  to  be  used  and  in 
what  order  they  are  to  be  executed. 

For  example,  PRELIMINARY  ANALYSIS  permits  him  to  insert  assumptions 
regarding  the  behavior  of  an  indeterminate  structure  and  thereby  to  make 
the  problem  determinate.  He  may  then  use  DETERMINATE  ANALYSIS  to 
obtain  a  complete  solution.  MEMBER  SELECTION,  based  on  previous^ 
design-criteria  input  and  permanently-stored  tables  of  standard  members, 
may  follow,  and  in  turn  be  followed  by  STIFFNESS  ANALYSIS,  which  pro¬ 
vides  a  rigorous  solution  for  the  indeterminate  structure.  With  SELECTIVE 
OUTPUT,  the  user  can  obtain  only  that  output  which  is  of  immediate  interest, 
lathei  than  voluminous  data  which  is  actually  generated. 

As  modifications  to  the  structure  are  made,  either  automatically  or  at  the 
user's  judgment,  new  analyses  are  executed  until  the  engineer  feels  that  an 
optimum  design  has  been  achieved.  This  type  of  on-line  operation  provides 
a  new  dimension  to  structural  design  and  gives  the  engineer  an  ability  for 
optimization  not  previously  possible. 

Seven  of  the  commands  listed  have  been  coded  and  implemented.  The 
remaining  three,  as  well  as  others  which  further  extend  system  capability, 
require  further  work. 
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Analysis  -  J.  Melvin  Biggs  and  Daniel  Belt  ran -Maldonado 

In  the  design  of  Civil  Engineering  structures,  the  problem  of  analyzing 
those  subjected  to  dynamic  actions  arises  frequently.  The  laborious  nature 
of  the  analysis  on  the  one  hand,  and  the  possibility  of  more  refined  methods 
on  the  other,  make  the  solution  of  these  problems  by  a  computer  very 
desirable.  Indeed,  the  technical  literature  indicates  extensive  use  of 
computers  in  this  connection,  mostly  through  special-purpose  programs 
coded  to  solve  individual  problems.  The  development  of  a  programming 
system  to  analyze  a  variety  of  structures  for  free  or  forced  vibration,  such 
as  space  frames,  space  trusses,  plane  grids,  plane  frames,  and  plan 
trusses,  is  in  progress.  This  system  is  a  problem-oriented  language 
called  DYNALS.  similar  in  many  respects  to  STRESS.  DYNALS  is  unique, 
not  only  because  it  is  applicable  to  a  wide  range  of  structural  configurations, 
but  also  because  it  encompasses  a  large  part  of  the  total  design  problem. 

Given  the  geometry  of  the  structure  and  the  stiffness  and  mass  of  individual 
elements,  it  forms  the  total  mass  and  stiffness  matrices,  obtains  eigen¬ 
values  and  eigenvectors,  and  computes  dynamic  response  of  the  structure 

in  minute  detail  for  any  specified  distribution  and  time  variation  of  applied 
forces. 

Three  aspects  of  the  problem  make  implementation  in  the  time-sharing 
mode  highly  desirable: 

1.  Iterative  nature  of  the  design  process  calls  for  a  simple 

and  rapid  way  of  modifying  a  problem  and  analyzing  the  effects  of  such 
modification  on  the  structural  response, 

2.  The  wide  variety  and  great  volume  of  possible  output  available 
from  a  dynamic  analysis  make  the  ability  to  request  at  will  selected  forms 
ol  output  very  desirable, 

3.  Iime-sharing  offers  unique  advantages  in  the  debugging  and  test¬ 
ing  process,  which  are  quite  valuable  to  a  user  who  is  working  with  such  a 
large  and  complex  program  block. 


DYNALS  will  be  incorporated  into  STRUDL,  which  deals  with  a  broader 
portion  of  the  design  process  for  Civil  Engineering  structures.  To  make 
the  system  efficient  for  both  large  and  small  structures,  the  various  data 
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arrays  are  dynamically  allocated  so  that  only  needed  arrays  are  required 
in  core  at  any  given  lime. 

DYNALS  accepts  and  decodes  data  pertaining  to  the  continuous  structure 
and  performs  a  "discretization"  of  the  stiffness  and  mass  characteristics. 
The  discrete  system  consists  of  a  number  of  degrees  of  freedom,  equal 
to  the  number  of  kinematic  degrees  for  a  joint,  times  the  number  of  joints. 
The  analyst,  presented  with  various  alternatives  regarding  the  formation 
of  the  mass  matrix,  chooses  \\ith  regard  to  the  nature  of  the  problem  at 
hand.  He  may  accordingly  reduce  the  number  of  degrees  of  freedom  to  a 
specified  selected  set.  A  wide  variety  of  output,  available  upon  rec]uest 
from  the  console,  ranges  from  eigenvalues  and  eigenvectors  of  the  free 
vibrating  structure  to  the  complete  response  of  one  and  all  of  the  mechani¬ 
cal  or  geometrical  quantities  of  the  structure  while  it  vibrates  under  a 
'time -varying  load  system. 

DYNALS  may  be  subdivided  into  the  following  programming  blocks: 

1.  INPUT  phase,  decoding  of  input  statements,  storage  of  data 
and  process  parameters,  and  consistency  checks, 

2.  Formation  of  the  logical  stiffness  matrix, 

3.  Formation  of  the  logical  mass  matrix  for  the  case  where  member 
masses  are  specified, 

4.  Reduction  of  the  number  of  degrees  of  freedom  to  a  selected  set 
when  specified, 

5.  Reduction  of  loading  data, 

6.  Solution  of  the  eigenvalue  problem, 

7.  Using  the  characteristic  shapes  obtained  in  6.  as  input  joint 
displacemeri.ts ,  the  preparation  of  sets  of  static  solutions  for  forces, 
reactions  and.member  distortions,  (n  is  the  number  of  modes  to  be 
combined.  ) 

I 

8.  Using  the  results  from  5.  and  6.  ,  the  numerical  integration  of 
the  modal  equations  of  motion,  the  results  of  which  are  used  to  combine 
the  static  solutions  developed  in  7.  ,  (Record  is  kept  of  maxima  and  time  of 
their  occurence  for  forces,  reactions  displacements  and  distortions.  ) 

9.  OUTPUT  phase. 

As  of  June  1965,  the  system  design  and  a  substantial  portion  of  the  program¬ 
ming  is  complete. 
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52iLE.e.„»Hna.l^len,-0,^„.ed  Lan.ua„  -  Ro,.„,  l,  Schlfto.„  and 

Laurence  N.  Beckreck 

SEPOL  (Soil  Englnaerlng  Pnoblem-Orlenled  Langu.ga,  m  a  sy.„„  „f 
mi.rr.lated  programs  „hrch  .er.e.  a  .hr.e-fold  purpoaa.  Flr.t,  SEPOL 
.»  Hosigned  as  a  communication  ay., am  i„  „hich  the  communication  be».,n 
man  and  machine  is  in  an  interrogation  form  tna,  elicits  the  program  from 
e  user  and  channels  the  program  towards  the  desired  objective.  Second 
It  1.  a  computer-aided  design  system  that  so  widens  the  area  of  choice 
and  breadth  o,  calculation  that  by  sheer  speed  and  utility,  a  greater  range 
of  calculation  ..  open  to  the  soil  engineer.  Third.  1,  performs  a  basic 

aid  Lr  analysis.  By  subdividing  the  system 

and  ^  »“■'  kyPhlha.as 

and  analytical  units. 

SEPOL  I,  covered  in  this  report,  i,  a  system  designed  to  calculate,  via 
time-.har.ng,  the  magnitude  and  progress  of  settlement  of  an  earth  mas, 

W  en  the  so.l  surface  is  subjected  to  a  specified  loading.  I,  was  designed 
as  a  prototype  to  determine  the  feasibility  and  usefulness  of  this  type  of 
computing  system,  and  also  to  provide  useful  informalion  for  soil  engineers. 

SEPOL  I  is  composed  of  many  small  subroutines,  each  of  which  represents 
a  unique  task  that  is  always  pe.formed  as  a  unit  and  can  be  performed  any- 
Wher.  and  many  times  in  an  analytical  sequence.  The  control  function  of 
IS  interrogatory.  A  se,  of  questions  1,  posedi  and,  since  the 
lespon.e  to  the  questions  provides  the  linkage  belween  bores,  the  set  o' 
questions  and  answers  builds,  for  each  situation,  a  specific  main  program, 
SEPOL  was  designed  as  an  easily-erpandable  system.  To  add  a  new 
method  of  analysis  requires  fits,  the  programming,  in  FORTRAN,  of  tasks 
or  calculation,  no,  already  in  the  system,  and  then  the  augmentation  of  a 
c  ionar>  of  vaiiables.  Access  to  the  new  system  is  accomplished  by 

adding  choices  at  the  appropriate  decision  points. 
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l/ O  System  Research  -  Ken  Reinschmidt,  D.  A.  Gardner,  and  John  W.  Weber 

The  hardware  which  connects  the  IBM  1620  to  the  MAC  time-sharing 
system  consists  of  two  logically-different  parts.  First,  a  general  l/O 
interface  to  the  1620  makes  it  easier  to  connect  non-IBM  devices,  and 
second,  a  dataphone  interface  which  goes  between  the  general  interface 
and  the  dataphone  which  is  connected  to  MAC.  Because  of  certain  idio¬ 
syncrasies  in  the  1620  l/ O  control,  such  as  the  lack  of  interrupt,  time¬ 
sharing  operation  must  be  under  software  control.  Subroutines,  written 
in  1620  assembler  language  and  available  to  the  FORTRAN  programmer, 
allow  him  to  use  the  1620  typewriter  and/or  card  read/punch  as  a  CTSS 
console,  pass  numeric  data  between  programs  in  7094  and  the  1620,  pass 
alphanumeric  data,  and  perform  alphanumeric  comparisons  on  data  in  the 
1620.  The  subroutines  were  developed  in  a  modular  form  which  facilitates 
modification  or  expansion  of  the  present  system.  Mr.  Weber  also 
developed  a  supervisor  subroutine  for  the  1620,  which  gives  a  more  auto¬ 
matic  operation.  The  main  uses  of  the  system  have  been  for  obtaining 
plotted  results  of  data  developed  at  the  ESL  Console  and  for  card  l/O  of 
lengthy  I/O  operations. 

Another  improvement  on  the  1620  v/as  made  by  the  addition  of  the  Gerber 
VP-600  plotter  as  an  output  device.  Tnis  plotter  was  put  on-line  with  the 
1620  in  August,  1964,  through  a  modification  of  the  old  CalComp  Plotter. 

Several  plotting  Subroutines  were  developed  for  the  Gerber  unit:  line, 
number  and  character  plotting  routines  by  Mr.  Gardner;  and  then  curve 
and  circle  routines. 

An  on-line,  date -acquisition  system  has  been  developed  jointly  by  members 

of  the  System's  Laboratory  and  the  Structural  Model's  Laboratory.  The 

hardware  part  of  the  system,  connected  to  the  1620  through  the  general 

interface  mentioned  above,  consists  of  electro -mechanical  switches, 

digital  logic,  signal  conditioning  and  amplification,  and  an  analog-to- 

digital  converter.  The  system  allows  the  1620  to  switch  to  any  one  of  50 

strain  gauges  and  read  the  converted  voltage  into  storage.  ^ 

1 


32 


CIVIL  ENGINEERING  DEPARTMENT 


A  large  set  of  software,  which  facilitates  the  use  of  the  system,  is  a 
problem-oriented  language  for  the  data  acquisition  system.  With  this 
language,  the  user  can  do  any  of  the  following: 

1.  Balance  any  or  all  gauges, 

2.  Select  how  many  and  which  gauge  he  desires  to  read, 

3.  Read  each  gauge'  as  many  times  as  desired, 

4.  Do  simple  calculations, 

5.  Convert  his  data  to  several  different  forms, 

6.  Plot  results, 

7.  Store  intermediate  results  on  the  disk, 

8.  Change  loading  and  repeat  operations. 

The  system  was  developed  by  Mr.  Reinschmidt,  has  been  used  in  class¬ 
room  demonstrations  and  on  some  plate  studies,  and  is  now  being  used 
extensively  on  models  of  a  dam  under  construction  in  Spain. 
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Bridge  Design  -  John  F.  Brotchie,  Charles  A.  Cornell,  Jose  M.  Roesset, 
and  Guillermo  Guzman  Barron-Torres 

The  purpose  of  this  research  is  to  develop  an  integrated  and  comprehensive 
system  for  the  directed  design  and  drafting  of  highway  bridge  structures. 

It  is  similar  in  several  respects  to  the  Structural  Design  Language  reported 
elsewhere,  but  is  oriented  specifically  toward  the  design  of  highway  bridges. 

The  design  system  proposed  consists  of  a  series  of  logical  blocks,  self- 
sufficient  in  themselves.  These  blocks,  each  one  consisting  of  a  large  and 
expandable  number  of  links,  are  to  be  interconnected  and  each  is  amenable 
to  direct  input  and  output.  Individual  blocks  perform  the  following: 

1.  An  internal  definition  of  geometry; 

2.  A  determination  of  suitable  alternate  forms,  span  arrangements 
and  materials; 

3.  A  design  of  superstructure,  piers,  abutments,  approaches,  or 
parts  thereof; 

4.  An  analysis  of  the  structure  or  components  above; 

5.  Optimization; 

6.  Detailing; 

7.  Cost  analysis  and  calculation  of  quantities  to  various  degrees  of 
detail,  depending  on  the  stage  of  the  design  at  which  the  link  was  entered; 

8.  Documentation,  including  drafting  and  specification  preparation. 
The  point  of  entry  to  the  system,  the  flow  between  the  blocks,  and  the  kind 
of  exit  may  be  chosen  by  the  designer  at  the  beginning  of  the  process,  but 
are  preferably  chosen  during  it,'  Time-sharing  facilitates  this  choice,  and 
a  large  computer  facility  provides  the  storage  necessary  to  retain  all 
required  information  throughout  the  process  and  communicate  it  to  the 
designer  in  the  most  efficient  way,  whether  by  plotter,  printer,  teletype 

or  oscilloscope.  Thus,  the  logic  of  the  design  process  and  the  basis  for 
determining  design  parameters  is  controlled  by  the  designer. 

Some  phases  of  the  system  are  already  completed  and  were  recently  linked 
by  Guillermo  Guzman  to  demonstrate  a  time-sharing  version  of  the  program. 
These  were  the  input  phase,  the  definition  of  internal  geometry,  the  selection 
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of  alternates  for  straight  bridges,  and  the  design  of  simple-span, 
composite,  precast-prestressed,  concrete  superstructures  together 
with  their  optimization  on  the  basis  of  a  comprehensive  minimum  cost. 

A  pilot  version  of  the  steel  design  phase  will  also  be  added  shortly. 

Optimal  Synthesis  of  Road  Networks  -  Alan  M.  Hershdorfer 

This  study  was  originally  an  investigation  into  the  properties  of  optimal 
traffic  flo'  •  patterns  on  metropolitan  road  networks.  A  total-travel-time 
function  to  be  minimized  is  expressed  as  a  piecewise,  linear,  convex 
approximation  to  the  experimental  travel-time  vs.  traffic-volume  relation¬ 
ship.  The  linear  programming  code  RSMFOR  was  used  to  determine  cost 
flow  solutions  for  a  multi -commodity  flow  system.  The  model  was  then 
extended  to  a  mixed-integer  programming  formulation.  .  The  integer 
variables  represented  decision  variables  related  to  the  provision  or 
deletion  of  road  capacity.  In  this  model,  numerous  related  problems  of 
^t’affic  network  synthesis  are  studied.  The  function  to  be  minimized  is  a 
linear  function  of  construction  cost  and  user  cost. 

The  simplest  synthesis  problem,  that  in  which  construction  cost  is  nil, 
occurs  in  the  problem  of  optimally  orienting  a  set  of  one-way  streets. 
Integer-linear  constraints  on  lane  allocation  require  all  lanes  in  a  road 
link  to  be  oriented  in  one  of  the  two  directions  or  to  be  split  equally  between 
the  two  directions.  This  one-way  street  problem  was  studied  extensively 
with  a  time-sharing  version  of  the  Land  and  Doig  Algorithm  for  mixed- 
integer  programming,  which  prescribes  values  for  integer  variables  in  an 
ordered  sequence  of  linear  programs.  With  time-sharing,  the  algorithm 
was  monitored  during  computation  and,  thus,  trial  solutions  were  very 
rapidly  evaluated.  Since  the  structure  of  the  flow  constraints  forced 
certain  variables  to  near  integer  values  after  a  small  number  of  integer 
decisions  had  been  made,  a  good  primal  feasible  integer  solution  was  very 
rapidly  determined.  Time -sharing,  used  as  a  monitor  on  the  algorithm, 
has  enabled  us  to  develop  a  heuristic  method  for  quickly  obtaining  solutions 
to  problems  that  involve  the  synthesis  of  flow  networks. 
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Sequences  and  the  Four-Color  Problem  -  David  A.  Huffman 

Methods  for  the  generation  of  binary-valued  pseudo-random  sequences, 
which  are  more  general  than  shift-register  sequences,  have  been  studied. 
These  sequences  have  the  interesting  property  that  they  remain  invariant 
to  a  type  of  transformation  which,  in  a  manner  reminiscent  of  tag 
systems,  -substitute  one  string  of  symbols  for  another. 

The  four-color  problem  (a  famous,  as  yet  unsolved  problem  in  topology) 
has  been  studied  with  the  view  that  it  furnishes  a  difficult  and  realistic 
example  of  a  problem  for  which  a  man,  aided  by  appropriate  program¬ 
generated  computer  displays,  may  be  able  to  gain  insight  more  rapidly 
than  a  man  working  with  pencil  and  paper.  The  maps  associated  with  this 
problem  are  nicely  suited  to  representation  by  certain  simple,  symmetric 
list  structures. 

A  Design  Language  for  Digital  Systems  -  Gerald  J.  Burnett  and 

Jack  B.  Dennis 

This  research  was  promoted  by  the  need  for  a  more  consistent  and  con¬ 
venient  method  of  designing  digital  systems  and  by  the  desire  for  computer 
aid  in  this  design  process.  We  have  developed  a  design  language  to  meet 
a  number  of  requirements.  Some  of  the  more  basic  requirements  are  the 
following; 

1.  The  language  must  be  complete,  that  is,  it  must  be  able  to 
describe  the  logic  and  timing  of  any  digital  system; 

2.  The  language  must  be  convenient  for  the  designer  and  must  be 
simple  enough  to  express  all  of  the  common  hardware  logic  and  timing 
functions ; 

3.  The  language  must  have  features  analogous  to  macros  and 
repeats  in  assembly  languages  in  order  to  describe  readi  1  y -ite rated  logi¬ 
cal  structures; 

4.  The  language  must  Ijc  alile  to  describe  a  system  in  sections  or 
components  that  are  designed  by  different  individuals. 
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As  an  example  of  the  description  of  a  familiar  logical  structure  in  the 
design  language,  we  describe  the  logic  required  to  perform  parallel 
binary  addition  of  a  number  contained  in  register  B  to  a  number  con¬ 
tained  in  register  A  and  leave  the  result  in  register  A. 

The  addition  method  used  for  the  example  involves  two  steps.  When  the 
pulse  named  "partial-add"  occurs,  each  bit  of  register  A  is  complemented 
if  the  corresponding  bit  of  register  B  is  one.  The  pulse  named  "carry" 
complements  a  stage  of  A  if  there  is  a  carry  signal  into  that  stage.  The 
description  of  one-  stage  of  the  adder  is  as  follows: 

begin  partial-add;  B[  i  ]  T  A[  i  ]  ; 

carry;  if  C[  i  ]  then  f  A[  i  ]  ; 
whe re:  i  4  0  them 

C[i-1  I  =  C[  i  I  A  A[i  ]  V  m  C[i  ]  A  B[i  ]  ; 


The  combinational  logical  operators  A,  V  and  m  have  their  usual  signifi¬ 
cance,  while  the  upward  arrow  denotes  the  action  of  complementing,  as  a 
unary  operator,  or  of  conditional  complementing,  as  a  binary  operator. 

The  first  two  statements  of  the  example  show  ways  of  indicating  that  cer¬ 
tain  actions  arc  to  take  place  on  the  occurrence  of  definite  pulses,  cither 
absolutely  as  in  the  first  statement  or  conditioned  by  a  logic  level  as  in 
the  second.  The  third  statement  illustrates  several  features  of  the 
"language.  The  where  clause  indicated  that  the  statement  following  then 
is  to  be  considered  part  of  the  description  only  if  the  Boolean  expression 
i  ^  0  is  true.  In  this  instance,  the  following  statement  represents  combina- 
tional'logic  that  generates  level  C[i-l  1  according  to  the  expression  at  the 
right  of  the  equivalence  sign. 

The  complete  description  of  the  adder  involves  ifi  ration  of  the  above 
logic  for  each  position  in  the  two  registers  and  is  handled  by  a  for 
statement  as  sliown  below.  Ciertain  declarations  are  included  so  that  the 
types  of  named  objects  are  unambiguous.  The  integer  quantity  n  is  the 
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number  of  bits  in  registers  A  and  B  and  is  presumed  to  be  assigned  a 
definite  value  outside  this  portion  of  the  description. 


begin 

integer  n; 

register  A[  0,  n-I  ]  ,  B[  0,  n-l]  ; 
level  '  C[  0,  n-I  ] 
add:  partial  add; 

delay  500n; 
carry; 

C[ n-I  ]  =  0; 

for  i  =  0  through  n-l  then 

begin  partial  add:  B[  i  ]  T  A[  I  ]  ; 

carry:  if  C[  i  ]  then  t  A[  i  ]  ; 
where  i  ?!  0  then 

C[  i-1  ]  =  C[  i  ]  A  A[  i  ]  V  “i  C[  i  ]  A  B[  i  ]  ; 

end; 

end; 


The  pulse  add  performs  the  partial  addition  and  initiates  a  500-nano- 
sccond  delay.  At  the  end  of  the  delay,  the  carry  pulse  occurs.  The 
meaning  of  the  for  statements  is  intrepreted  as  if  the  (compound)  state¬ 
ment  following  then  were  written  out  successively  for  integral  values  of 
i  from  0  through  n-l. 


Work  related  to  the  design  language  will  be  devoted  next  to  the  develop¬ 
ment  of  suitable  simulation  techniques  and  means  whereby  a  designer  can 
employ  the  language  on-line  in  the  process  of  creating  digital  system 
plans. 
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Man -Machine  Communication  -  Joseph  VVcizenbauiTi 


SLIP 

SLIP,  a  list-processing  system  capable  of  being  imbedded  in  a  large 
class  of  algebraic  compiling  languages,  was  imbedded  in  the  MAD  system. 
MAD  is  the  compiling  system  most  frequently  used  in  CTSS,  and  there¬ 
fore,  SLIP  has  been  made  available  to  large  sections  of  the  MAC  user's 
community.  A  number  of  students  and  faculty  members  have  made  use 
of  SLIP  for  research  and  thesis  purposes. 


PPL 

The  on-line  SLIP-based  programming  system  OPL  was  rewritten  in  the 
MAD-SLIP  system  discussed  above.  It  has  served  as  a  test  bed  for  a 
number  of  e.xperiments,  particularly  te.xt  manipulation,  and  ha's  been  an 
invaluable  tool  in  clarifying  issues  underlying  interactive  progfamming 
systems,  their  design  and  implementation.  The  problem  of  the  program 
ming  system-user  interface  has  been  illuminated,  but  not  solved. 


ELIZA  ? 

A  major  effort  in  the  design  of  a  natural-language,  on-line,  man/machinc 
conversational  system  was  begun  and  carried  forward  during  the  report¬ 
ing  pc,  riod.  At  this  time,  ELIZA  is  capable  of  conversing  with  users  in 
natural  language  in  certain  restricted  conte.xts.  We  arc  beginning  to 
expand  ELIZA  so  that  it  can  form  the  basis  for  both  information  gathering 
and  a  retrie''al  system  operating  on  a  natural-language  basis.  ELIZA 
was  also  used  for  experiments  on  two-person  conversational  interaction 
at  the  Massachusetts  General  Hospital.  The  development  of  ELIZA 
required  implementation  of  powerful  text  manipulation  functions  within 
MAD-SLIP.  rhese  latter  are  of  general  use  as  external  MAD  functions. 
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Research  on  the  Theory  of  Automata  -  Robert  McNaughton 

A.  CLASSIFICATION  AND  COMPLEXITY  OF  FINITE-STATE  OPERATIONS 
The  mos!  commonplace  distinction  in  the  theory  of  automata  has  been 
that  between  finite-state  operations  and  operations  requiring  potentially- 
infinite  automata.  Partial  motivation  of  our  research  has  been  the 
rontention  that  finite-state  operations  should  themseives  bo  classified  in 
some  interesting  manner,  Stinaulated  by  some  woric  of  Ivl.  P.  Schutx.cnberger , 
we  have  proposed  a  notion  that  we  think  gives  an  important  classification, 
namely  the  notion  of  a  "topological"  event  (a  term  which  we  do  not  feel  is 
optimal) . 

Interestingly,  the  notion  can  bo  characterized  in  several  distinct  ways: 

1)  as  an  event  whose  reduced-state  graph  contains  no  loop  that  counts  any 
word  modulo  m,  for  m  >;  2)  as  an  event  expressible  as  a  regular  expression 

containing  free  use  of  all  the  Boolean  operators,  but  without  star  (denot¬ 
ing  closure,  or  iteration);  3)  as  an  event  expressible  in  a  certain  first- 
order  language  of  symbolic  logic;  4)  as  an  event  realized  by  a  nerve  net 
in  which  feedback  is  confined  to  the  case  in  which  output  of  a  neuron 
drives  one  of  its  own  inputs;  and  5)  as  an  event  whose  semi -group 
contains  no  non-lrivial  subgroups.  The  fact  that  this  class  of  events  can 
bo  characterized  in  so  many  different  ways  -  structurally,  behaviorally, 
and  algebraically  -  each  interesting  in  its  own  right  -  is  evidence  enough 
of  the  importance  of  the  concept.  A  long  paper  on  this  matter  is  planned, 
which  will  include  proofs  establishing  the  equivalence  of  the  five  charac¬ 
terizations  . 

There  are  two  interesting  ways  of  defining  the  loop  complexity  of  a 
regular  event,  each  in  terms  of  a  class  of  regular  expressions.  The 
first  class  is  that  of  restricted  regular  expressions,  i.  e.  those  made  up 
from  the  operators  union,  concatenation,  and  star.  The  .°ocond  class 
consists  of  those  made  up  of  those  operators  and,  in  addition,  inter¬ 
section  and  complementation;  such  expressions  are  called  "general 
regular  expressions".  Thus  general  regular  expressions  contain  all  the 
Boolean  operators,  while  restricted  regular  expressions  contain  only 
union.  In  each  case,  the  loop  complexity  of  an  event  is  defined  as  the 
minim.'il  star  height  of  all  regular  expressions  in  the  class  that  represent 
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the  event.  (The  star  height  of  a  regular  expression  is  the  maximum 
length  of  a  nested  chain  of  the  stars.  )  Thus,  for  example,  the  loop  com¬ 
plexity  of  a  topological  event  with  respect  to  general  regular  expressions 
IS  zero.  But  there  are  topological  events  with  arbitrarily  large  star 
height  with  respect  to  restricted  regular  expressions. 

It  seems  that  the  most  useful  tool  for  investigating  star  height  with 
respect  to  restricted  regular  expression  (which  is  the  problem  studies 
originally  by  L.  Eggan)  is  graph  theory.  On  the  other  hand,  the  most 
useful  tool  for  investigating  star  height  with  respect  to  general  regular 
expressions  is  the  theory  of  finite  semigroups.  Both  concepts  seem 
interesting  and  important,  but  both  have  unsolved  problems.  Thus,  wo 
do  not  even  know  whether  there  exists  an  event  of  star  height  3  with 
respect  to  general  regular  expressions;  and,  although  wo  do  know  the 
existence  of  events  with  arbitrarily  large  star  height  with  respect  to 
restricted  regular  expressions,  wo  do  not  know  whether  there  exists  an 
algorithm  for  determining  this  star  height  when  a  regular  event  is  given. 

B,  CONSTRUCTION  OF  WELL-TIMED  NETS  FROM  BADLY-TIMED 
ELEMENTS 

Research  in  this  area  has  been  a  continuation  of  research  at  the 
University  of  Pennsylvania.  A  seminar  was  held  consisting  of  some 
members  of  the  Machine  Structures  Group  of  Project  MAC  and  some 
design  engineers  from  Honeywell.  The  number  attending  was  small,  but 
the  discussions  were  intensive  and  fruitful.  Some  improvements  were 
made  in  the  report  written  at  the  University  of  Pennsylvania  (Badly -Timed 
j:iemcnts  and  Well-Timed  Nets,  by  R.  McNaughton).  This  report  will  be 
issued  as  a  monograph  next  year. 

C.  INFINITE  HISTORY  OF  FINITE -STATE  MACHINES 

Two  widely-differing  studies  in  the  literature  have  led  into  this  area. 

J.  R.  Buchi  studied  a  certain  system  of  symbolic  logic,  and  developed  a 
decision  procedure  by  showing,  in  effect,  that  every  formula  in  the 
system  says  something  about  the  iiifinitc  history  of  a  finite-state  machine. 
And  D.  Muller  analyzed  the  behavior  of  an  asynchronous  circuit  by  regard¬ 
ing  the  circuit,  which  is  subject  to  varying  input  conditions,  and  in  which 
the  relative  timing  of  the  parts  is  unpredictable,  as  a  finite  automaton 
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with  an  infinite  duration  in  time:  a  buzzer,  for  example,  would  be  con¬ 
ceived  of  as  a  two-state  automaton  going  from  one  state  to  the  other  ad 
infinitum. 

The  proof  of  an  analogue  to  the  so-called  Kleene  Myhill  theorem  for  the 
infinite-duration  problem  has  been  written  up,  and  will  be  submitted  for 
publication.  Also,  there  is  an  interesting  connection  between  this 
problem  area,  and  the  topic  of  infinite  games;  i.  e.  ,  games  which  two 
players  play  in  turn  an  infinite  number  of  times. 

It  would  appear  that  this  problem  area  is  more  abstruse  than  the  other 
problem  areas.  In  spite  of  the  fact  that  the  area  was  approached  by  two 
widely  separated  studies,  no  new  application  seems  to  be  forthcoming. 
We  are  far  from  the  position  that  this  area  is  worthless;  however,  we 
do  think  it  deserves  less  attention  than  the  other  two. 

Program  Segmentation  -  Richard  Y.  Kain  and  David  J,  Kuck 

A  question  of  major  importance  is  segmentation  of  users  programs  to 
decrease  memory  space  and  swap  times.  There  are  several  problems 
which  immediately  arise.  The  major  one  relates  to  efficiency  of  the 
algorithm  which  might  be  developed  to  perform  such  segmentation.  In 
particular,  is  the  time  spent  segmenting  a  program  greater  than  the 
time  saved? 

Four  approaches  to  this  problem  are: 

1.  The  compiler  examines  a  program's  gross  structure  and 
makes  simple  decisions  about  where  to  divide  the  program; 

2.  The  compiler  divides  a  program  into  small  "atoms"  and 
then  attempts  to  combine  the  "atoms"  into  segments; 

3.  The  compiler  divides  programs  into  small  segments,  and 
then  adds  to  each  segment  a  small  program  which  determines  those 
segments  which  might  be  brought  into  memory  at  the  same  time; 

4.  The  compiler  describes  possible  paths  through  a  program  in 
terms  of  regular  expressions,  which'.are  examined  when  allocation 
decisions  are  made. 
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The  first  approach  is  difficult  to  implement  without  some  subdivision  of 
a  program;  therefore,  it  will  be  considered  as  a  special  case  of  the 
second  method. 

The  second  method  "reduces"  the  segmentation  problem  tc  a  clumping 
problem;  given  a  set  of  atoms,  which  are  most  closely  related  to  ea^h 
other?'  A  number  of  approaches  to  this  problem  have  been  considered, 
but  all  unfortunately  assume  that  there  is  some  way  to  measure  "inter¬ 
relationships"  between  pairs  of  atoms;  thus  reducing  the  interrelation¬ 
ship  to  a  number.  We  will  assume  that  these  numbers  have  been  placed 
in  a  matrix.  A,  where  A.^  is  the  interrelationship  between  atoms  i  and 

j.  Some  clumping  strategies  are: 

a)  Build  clumps  by  adding  "most  suitable"  neighbors  until 
some  criterion  of  maximum  size  or  minimum  suitability  is  met; 

b)  Partition  matrix  A  until  the  sum  of  the  elements  not  in 
the  matrices  along  the  main  diagonal  is  small; 

c)  Take  random  walks  through  the  "atoms,  "  letting  the 
probability  of  a  path  be  proportional  to  the  interrelationships  of  the 
elements  at  the  ends  of  the  path,  with  termination  by  length  of  the  walk; 

d)  Find  the  largest  eigenvector,  whose  components  are  the 
probabilities  of  visits  to  the  corresponding  states  in  the  long  run,  and 
clump  according  to  the  size  of  the  components  of  the  eigenvector. 

The  third  approach  is  a  crude  attempt  towards  "dynamic  segmentation.  " 
By  including  a  simple  program  which  determines  neighbors  (possibly  by 
examining  the  parameters  of  the  running  program,  such  as  loop  limits) 
with  each  atom,  only  those  atoms  which  are  likely  to  be  used  during  the 
next  quantum  need  be  fetched.  If  the  decision  is  too  complex,  a  standard 
set  of  atoms  would  be  used. 

The  fourth  approach  is  another  attempt,  less  sophisticated  than  the  third, 
to  dynamically  segment  the  program.  Program  data  are  not  checked  in 
the  choice  of  atom  set  and  less  economy  is  achieved.  This  method  is 
superior  to  a  priori  segmentation,  since  the  beginning  atoms  in  a  segment 
might  not  be  needed  if  the  previous  execution  had  halted  in  the  middle  of  a 
segment. 


a 


COMPUTER  COMMUNICATION  STRUCTURES 


45 


Progress  has  been  made  in  finding  interrelationships  among  methods  2a 
through  2d,  though  no  useful  generalizations  are  immediately  possible. 
Unfortunately,  evaluation  must  depend  upon  detailed  experience  with 
each  scheme  or  simultations  which  are  necessarily  very  dependent  upon 
the  models  of  user  programs. 

Waveform  Transformation  and  Graphical  Display  -  Thomas  G.  Stockham,  Jr. 

A  subroutine  for  the  plotting  of  numerical  data  on  the  E.S.  L.  display  was 
completed  by  Henry  Ledgard.  The  plotting  process  gives  freedom  from 
details  of  scale-factor  selection  and,  at  the  same  time,  produces  graphs  which 
are  labeled  decimally  in  convenient  round  numbers  and  may  have  logrithmic 
or  linear  scales  on  either  axis.  Subprograms  that  perform  various  sub¬ 
tasks  for  the  logical  specification  of  completed  plots  are  available  as  a 
separable  package.  These  routines  are  free  from  constraints  imposed  by 
the  specific  display  unit  for  which  they  are  intended  and,  thus,  provide  a 
general  service  for  any  plotting  equipment.  Their  main  subtasks  are 
generation  of  numerically-convenient  plot  boundaries  and  associated 
scaling  parameters,  generation  of  normalized  grid-line  parameters, 
generation  of  text  for  dimensioning,  and  scaling  of  data  to  fit  the  normal¬ 
ized  plotting  space.  Finished  plots  are  produced  from  normalized, 
floating-point  array  pairs  of  arbitrary  length.  The  total  generation  time 
is  between  one  and  two  seconds  on  the  IBM  7094. 

An  algorithm  has  been  developed  by  W.  R.  Chiodi  with  sufficient  flexibility  to 
implement  a  variety  of  linear  transformations,  such  as  Fourier,  Hilbert, 
and  convolution  transformations.  H.  G.  Murray,  Jr.  has  realized  the 
proposed  algorithm  in  a  computer  program.  This  system  operates  on  a 
standard-parts  technique  by  which  a  required  variability  is  concentrated 
in  two  building -block  functions,  two  simple  sets  of  rules  for  combining 
them,  a  set  of  parameters,  and  a  simple  transformation  on  these  parame¬ 
ters.  Since  these  variables  can  be  specified  through  suitable  numerical 
tables  and  brief  sets  of  algebraic  statements,  the  bulk  of  the  computation 
^  is  provided  by  a  static  framework  into  which  they  are  embedded.  This 
iiwork  will  become  part  of  a  problem-solving  system  in  electrical  linear - 
system  theory.  (See  Chiodi,  and  also  Murray,  Appendix  B,) 
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A  Table-Directed  Translator  -  Chung  L.  Liu' 

Research  in  the  table -directed  translator  system  has  continued  for  the 
past  year.  In  order  to  further  investigate  and  evaluate  some  of  the  ideas 
conceived  in  the  course  of  our  research,  such  a  system  is  being  designed 
and  implemented. 

It  is  also  hoped  that  this  system,  when  completed,  will  e.Ktend  a  great 
deal  of  our  programming  powers.  While  a  table -directed  translator 
system  provides  users  with  the  flexibility  of  designing  their  own  program¬ 
ming  languages,  it  is  most  important  from  the  users'  point  of  view  to 
learn  to  use  it  with  ease.  This  is  what  we  stress  in  our  design.  Clear 
boundary  lines  between  modules  of  the  system  and  simple  input  languages 
are  emphasized. 

It  is  e.xpected  that  most  parts  of  the  sysirom  will  be  programmed  in  a 
higher -level  language,  so  that  reprogramming  work,  due  to  hardware 
transition,  will  be  minimized.  We  intend  to  use  the  MAD  language  for 
this  purpose, 

On-Line  Braille  -  Edward  L.  Glaser 

It  should  be  understood  that  Braille  is  not  a  simple  one-to-one  code  with 
normal  characters  of  the  English  language.  Rather,  it  is  reminiscent  of 
shorthand  in  that  many  standard  contractions  and  abbreviations  are  used. 
Although  these  contractions  and  abbreviations  are  not  phonetic,  the  rules 
for  their  use  are  often  context  dependent.  As  a  consequence,  a  computer 
is  necessary  to  convert  standard  teletype  code  to  a  modified  teletype  code 
that  is  adequate  to  drive  an  electric  Brailler. 

Although  this  covild  have  been  done  within  the  IBM  7094  computer  at 
Project  MAC,  it  was  decided  to  use  a  separate  satellite  computer,  so 
that  Braille  would  also  be  available  in  the  earliest  phase  of  the  new 
MULTICS  time-sharing  system  on  the  General  Electric  635  computer. 

As  a  consequence,  a  Digital  Equipment  Corporation  PDP-8  computer  was 
acquired  this  year,  since  it  is  a  small,  high-speed,  versatile  machine 
that  is  completely  adequate  for  the  job.  To  produce  Braille  output,  we 
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selected  an  electric  Braille  Writer  that  was  designed  and  developed  in 
the  M.l.  T.  Mechanical  Engineering  Department. 

The  program  for  the  conversion  from  teletype  code  to  Braille  has  been 
written  and  checked  out.  Also,  the  interface  between  the  PDP-8  and  a 
Data  Set  has  been  checked  out  so  that  the  Braille  system  can  communi¬ 
cate  with  the  MAC  system  as  though  it  were  a  standard  teletype  station. 
Still  remaining  are  final  check-out  of  the  interface  between  the  PDP-8 
computer  and  the  electric  Brailler  and  check-out  of  the  Brailler  itself. 
(See  Luconi,  Appendi.x  B.) 

Future  plans  include  use  of  the  PDP-8  machine  as  a  base  for  additional 
studies  of  man-machine  communication.  Some  of  the  possible  areas  of 
investigation  are:  auditory  coding,  interaction  of  auditory  and  visual 
stimuli,  and  visual  and  tactile  stimuli.  To  date,  a  method  for  producing 
arbitrary  auditory  wave  forms  has  been  designed  and  fitted  to  the  machine 
It  consists  of  two  6-bit  digital-to-analog  converters  that  are  driven  by  the 
machine.  The  output  analog  wave  form  is’  sent  through  low -pass  filters 
to  stereo  amplifiers  and  loud  speakers.  Programming  the  equipment  for 
this  application  has  not  been  completed, 

:^alysis  of  Time-Shared  Computer  Systems  -  Allan  L.  Scherr 

Some  aspects  of  the  operation  of  time -shared,  interactive  computer 
system  ■  were  analyzed  by  making  extensive  user  and  performance 
measurements  of  the  Project  MAC  system.  Emphasis  was  on  the 
reaction  of  hardware  systems  to  the  demands  that  users  make.  Simply 
stated,  the  problem  was  to  characterize  both  time-shared  systems  and 
their  users  in  order  to  predict  the  performance  of  the  two  operating 
together. 

Portions  of  the  problem  included  specification  and  measurement  of  user 
characteristics,  the  development  and  verification  of  both  simulation  and 
mathematical  models  for  time-shared  systems,  and  the  specification  and 
measurement  of  performance  metrics  for  such  systems. 
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First,  simulation  models  were  used  to  study  the  effects  of  changing  small 
details  in  the  operation  of  CTSS-like  systems.  Then,  a  continuous -time 
Markov  process  model  was  derived  to  predict  the  performance  of  a  broad 
class  of  systems.  Three  lengthy  programs  were  written  to  gather  CTSS 
data.  The  first  measured  interaction  parameters;  the  second  measured 
input,  output,  and  idle  times  of  the  system;  and  the  third,  statistics  for 
interaction  parameters  of  five  system  commands..  This  data  was  gathered 
on  the  MAC  system  between  October  1964  and  March  1965. 

The  CTSS  data  were  used  as  a  basis  for  comparison  with  model  pre¬ 
dictions,  and,  in  order  to  take  measurements  and  build  models,  many 
definitions  of  commonly  used  time-shared  terminology  were  made  precise. 
This  research  was  used  as  the  basis  for  a  doctoral  thesis.  (See  Appendix 
B  and  MAC-TR-18.  ) 

Queueing  Models  for  File  Memory  Operation  -  Peter  J.  Denning 

A  problem  which  has  received  much  attention  at  Project  MAC  is  the 
scheduling  of  user  programs.  The  Scheduling  Algorithm  assigns  each 
waiting  program  a  burst  of  processor  time  during  which  it  runs;  and  each 
program  continues  to  receive  bursts,  at  intervals,  until  it  has  run  to  com¬ 
pletion.  While  not  running,  it  is  stored  on  a  drum  or  disk.  The  process 
of  placing  a  program  in  core  memory  or  of  removing  it  from  core  and 
placing  it  on  the  drum  or  disk  is  called  swapping.  An  optimum  Scheduling 
Algorithm  results  in  the  most  efficient  system  operation  and  yet  treats  the 
users  fairly.  Now  that  computing  systems  are  about  to  be  built  in  which 
segments  of  many  users'  programs  are  simultaneously  present  in  core 
and  several  processors  operate  simultaneously,  another  problem  arises; 
scheduling  of  the  necessary  swaps  and  of  any  requests  generated  by  the 
various  user  processes  for  file  memory  use.  At  present,  requests  are 
made  singly  to  access  file  memory;  and  a  request  does  not  occur  until  the 
preceding  one  has  been  processed.  In  a  system  with  many  simultaneous 
processes,  requests  for  file  memory  use  will  be  placed  in  a  queue.  This 
research  attempts  to  find  a  reasonable  model  of  the  way  in  which  user  pro¬ 
cesses  demand  file  memory  use  and  a  reasonable  and  optimal  method  of 
handling  the  queued  requests. 
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It  is  not  possible  to  construct  a  detailed  model  for  a  computational  pro¬ 
cess  in  a  segmented,  multi-programmed  system,  since  no  such  system 
exists.  However,  we  can  make  reasonable,  probabilistic  assumptions: 
for  example,  that  segment  lengths  are  Poisson  distributed  about  some 
mean  number  of  pages;  that  inter-request  times  are  Poisson  distributed; 
and  that  a  given  process  will  be  able  to  generate  requests  while  its  write 
requests  are  in  se rvice, .  but  will  be  suspended  while  any  read  request  is 
in  service. 

Having  constructed  a  reasonable  model  for  user  processes,  we  had  to 
consider  its  request  for  file  usage.  The  methods  of  queueing  theory  were 
applied,  in  particular,  the  following  queue  disciplines: 

1.  fi  rst-come-first  served, 

2.  shortest-job-first, 

3.  sho  rtest-access -time-first. 

The  access  time,  defined  as  the  rotational  positioning  delay  of  a  drum  or 
disk  file,  is  measured  from  the  time  a  channel  becomes  free  until  the 
desired  starting  location  of  a  request  is  in  position;  it  is  meaningless  if 
there  are  no  requests  waiting  in  the  queue.  A  shortest-acces s -time- 
first  queue  was  found  to  be  the  most  efficient  of  the  three  queue  discipline,? 
or  any  combination  of  them.  Expressions  were  derived  for  the  average 
wait  in  queue  and  the  average  number  of  waiting  requests,  and  the  most 
efficient  queue  discipline  was  taken  to  be  the  one  for  which  over-all  idle 
time  was  minimized.  Over-all  idle  time  is  the  idle  time  of  all  user  pro¬ 
cesses,  plus  idle  time  of  all  processors,  plus  idle  time  of  all  file  memory 
channels.  A  method  was  assumed  to  exist  for  optimal  core  memory 
allocation,  that  is,  the  allocation  problem  for  core  is  irrelevant  to  the 
discussion  of  file  memory  scheduling,  although  we  considered  allocation 
schemes.  In  the  analysis,  simulation  was  found  to  be  a  most  useful  tool, 
because  it  showed  which  assumptions  were  valid  and  which  were  not.  It 
also  showed  the  form  of  probability  distributions  for  waiting  times  in 
queue,  average  number  in  queue,  and  so  on,  that  should  be  expected. 
Finally,  simulation  was  used  to  verify  the  predictions  of  the  mathemati¬ 
cal  model.  (See  Denning,  Appendix  B.  ) 
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Input/Output  Subsystems  -  Arthur  A.  Smith 

The  execution  of  input/output  instructions  in  a  multi- access  computer 
system  may  he  considered  as  either  simple  ordinary  instructions,  with 
long  execution  times  when  the  input/output  equipment  is  part  of  the  central 
processor,  or  as  a  form  of  communication  between  the  program  being 
executed  by  the  main  arithmetic  processor  and  a  program  written 
especially  for  the  I/O  subsystem  involved,  when  the  I/O  equipment  is  a 
separate  entity.  The  former  consideration  ignores  the  possibility  and 
the  desirability  that  the  central  processor  execute  an  entirely  separate 
program  while  the  I/O  is  occurring.  Execution  is  implicit  in  the  second 
viewpoint,  which  therefore  is  considered  more  fruitful. 

My  studies  regard  a  process  as  able  to  force  its  execution  by  a  particular 
processor,  in  particular  by  a  processor  suited  to  the  properties  of  the  1/ 0 
device.  That  is,  upon  execution  of  a  certain  instruction,  an  image  of  the 
state  word  of  the  process  is  transmitted  to  the  I/O  processor,  which 
modifies  the  state  word  by  the  execution  of  instructions  written  in  its  own 
order  code.  This  I/O  processor,  by  executing  an  instruction  analogous 
to  the  one  which  started  its  operation,  causes  resumption  of  the  usual 
sequence  of  instructions,  with  the  appropriately-modified  state  word. 

Thus,  it  executes  an  instruction  which  forces  the  subsequent  instructions 
to  be  executed  by  a  main  (arithmetic)  processor.  I  plan  to  publish  a 
!  report  on  possible  realizations  of  this  system  arid  alternatives  to  it. 

Non-Repeatability  of  Multi-Process  Computations  -  Earl  C.  Van  Horn 

A  MAC  system  consists  of  several  processors,  al  main  memory,  several 
file-storage  devices,  and  several  input/output  devices.  Tlie  supervisor 
program  of  such  a  system  schedules  available  processors  among  the 
various  jobs  to  be  done.  We  have  defined  the  notion  of  a  process  as  a 
precise  molding  of  the  job  to  be  done.  A  process  is  a  locus  of  control  or 
activity  within  an  instruction  sequence.  We  have  been  studying  the 
relationships  among  processes  and  various  other  system  entities  and  have 
concentrated  on  those  relationships  that  arc  invariant  with  respect  to  the 
way  the  supervisor  schedules  available  processors  among  the  processes. 

A  collection  of  processes  working  together  for  a  single  user  on  a 
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paniciilar  problem  forms  a  separate  computation.  All  of  the  processes 
of  a  separate  computation  may  reference  the  same  set  of  system  entities. 

I  his  set  constitutes  the  sphere  of  protection  associated  with  the  computa¬ 
tion. 

The  possibility  exists  in  a  MAC  system  that  different  runs  of  the  same 
eoniputation  from  the  same  initial  state  with  the  same  input  data  will  yield 
different  output  data.  .Such  a  computation  is  said  to  be  non- repeatable . 
Computations  running  in  the  present  MAC  system,  CTSS,  are  repeatable 
because  they  always  consist  of  a  single  process.  An  attempt  to  adapt 
CTSS  for  multi -process  computations  would  encounter  the  problem  of 
non- repeatability.  Because  in  practice  no  program  is  ever  debugged, 
repeatability  cannot  be  guaranteed  to  a  \iser  unless  sufficient  conditions 
for  repeatability  are  enforced  by  the  hardware  and  supervisor  of  the 
system.  These  sufficient  conditions  may  be  published  as  a  set  of  rules 
for  tisers  to.  follow  in  the  writing  of  programs  or  they  may  be  followed  by 
a  compiler  in  the  generation  of  multi-process  programs.  In  such  cases, 
the  burden  of  establishing  repeatability  lies  on  the  person  debugging  the 
program  or  compiler,  respectively.  In  order  to  study  the  interactions 
of  processes  i,n  multi-process  computations,  we  have  defined  a  Multi- 
Process  Automaton  of  Type  A  as  a  model  of  a  multi-process  computation. 

\\  e  have  postulated  one  set  of  conditions  sufficient  to  make  a  run  of  such 
an  automaton  repeatable.  A  rigorous  proof  of  this  hypothesis  is  forth¬ 
coming.  Later,  we  hope  to  formulate  and  prove  the  necessary  conditions 
for  repeatability  in  such  an  automaton. 

Semantics  for  Multiprog  rammed  Computations  -  Jack  B.  Dennis  and 

Earl  C.  Van  Horn 

The  description  of  computations  perlormed  by  a  multiprogrammed  computei* 
system  requires  functions  not  available  in  existing  programming  languages. 
In  particular,  functions  related  to  parallel  programming,  protection  of  data 
and  procedure  information,  and  communication  among  independent  compu¬ 
tations  are  usually  absent.  IV  e  have  .studied  the  problem  posed  by  these 
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semantic  aspects  of  computations  within  time-shared  computer  systems 
and  have  proposed  a  set  of  meta -instructions  through  which  these  essential 
facilities  may  be  realized.  Some  of  the  meta-instructions  defined  already 
exist  in  multiprogrammed  computer  systems  as  monitor  call  operations; 
however,  many  of  the  meta-instructions  are  yet  to  be  made  part  of  any 
computer  system.  The  operations  performed  by  meta-instructions  are 
based  on  the  program  structure  developed  in  our  previous  work.  A  com¬ 
putation  is  defined  by  a  collection  of  capabilities  represented  by  entries  in 
a  C-list.  The  process  active  in  a  computation  may  invoke  capabilities 
represented  in  its  C-list,  for  instance,  the  addressing  of  a  procedure  or 
data  segment  or  the  executing  of  an  I/O  function. 

A  group  of  the  proposed  meta-instructions  relates  to  the  control  by  one 
computation  of  other  inferior  computations.  This  hierarchical  relation¬ 
ship  of  computations  is  essential  to  program  debugging  within  time-sharing 
systems  and  to  the  proper  handling  of  faults  or  exceptional  conditions 
encountered  by  a  computation  process.  Other  meta -instructions  relate  to 
the  accessing  and  control  of  a  directory  structure  for  retained  objects, 
principally  segments  of  procedure  and  data.  We  suggest  a  hierarchical 
directory  structure  associated  with  each  intelligence  that  uses  the  computer 
system  in  order  to  facilitate  the  choice  of  unique  identifiers  for  retained 
objects  and  make  the  sharing  of  objects  among  computations  feasible. 

Optimal  Allocation  of  System  Resources  -  Robert  L.  Potter 

The  problem  of  optimal  memory  and  processor  allocation  in  a  multi¬ 
processor,  time-shared  computer  is  being  studied.  It  is  formulated  as  a 
two-space  mapping  problem  with  emphasis  on  the  use  of  probabilistic 
methods.  A  computer  has  a  certain  number  of  memory  and  of  processor 
units  available.  A  two-dimensional  space  can  be  visualized.  I  constructed 
tvvo'spaces  with  the  same  normal  axes,  a  memory  unit  axis  and  a  processor 
unit  axis;  one  is  called  "resource"  space,  the  other,  "needs"  space.  The 
resource  needs  of  users  are  distributed  in  needs  space.  An  allocation 
algorithm  assigns  each  user  to  a  portion  of  resource  space.  'This  is  a 
mapping  from  needs  space  onto  resource  space,  and  the  mapping  function 
must  assign  resources  equitably  and  efficiently. 
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The  probability  densities  of  memory  and  processor  needs  per  user  will 
not  be  known  until  measurements  can  be  made  on  an  appropriate  system. 
We  expect  the  densities,  which  are  time-dependent,  to  have  the  same 
general  appearance;  they  are  zero  at  zero  need,  rise  to  some  maximum, 
and  then  fall  off  toward  zero  as  the  need  becomes  infinite. 


The  mapping  function  is  not  easy  to  determine;  it  must  include  cost, 
efficiency,  and  user-satisfaction  factors  reflected  in  the  way  distribution 
in  needs  spaces  adapts  to  the  system.  Cost  factors  may  be  defined  in 
many  ways;  as  the  allocation  scheme  which  minimizes  cost,  in  one  case; 
or  wasted  time  or  resource  units  in  another.  User-satisfaction  factors 
are  elusive  of  definition  until  a  required  system  is  built. 

The  mapping  must  be  time-dependeni,  because  the  present  MAC  system 
considers  allocation  primarily  in  terms  of  program  length,  which  dis¬ 
criminates  against  long  programs.  We  need  a  more  sophisticated  approach 
in  which  the  running  status  of  a  user  depends  not  only  on  program  length, 
but  also  on  the  nature  of  the  resources  employed  since  the  beginning  of 
computation.  Thus,  a  mapping  function  must  consider  a  user's  position 
in  both  resource  space  and  needs  space  for  the  previous  few  moments  and 
the  present  moment.  (See  Potter,  Appendix  B.  ) 

Automatic  Flowcharting  -  Daniel  U.  Wilde 

The  objective  of  this  research' he  development  of  a  tool  to  aid  advanced 
programmers  in  analyzing  and  debuggin  digital  computer  programs 
written  for  general-purpose,  single-address  machines.  Input  to  the 
analysis  program  is  a  BCD  output-listing  tape  from  the  7094  FAP 
assembler.  Output  from  the  analysis  program  will  be  a  flow  chart  in 
symbolic  form  showing  the  effects  of  control  flow  and  data  flow.  The 
following  development  work  has  been  completed: 

A.  DATA  GATHERING 

On  the  first  pass,  the  analysis  program  reads  a  BCD  assembly  tape  and 
produces  the  following  tables  as  a  function  of  the  instruction  location: 

a.  Entry  and  Exit  transfer  tables, 

b.  Active  and  Passive  reference  tables, 
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c.  Data  and  Storage  pseudo-operation  tables, 

d.  A  Symbol  table, 

e.  An  Entry-Point  or  Starting-Location  table,  and 

f.  A  First-Instruction  and  Last-Instruction  location  table. 


B.  DATA  REDUCTION 


The  second  pass  of  the  analysis  program  processes  the  tables  constructed 
during  the  first  pass.  Examples  of  the  processing  are: 

a.  The  Entry  and  Exit  tables  arc  used  to  divide  the  program  into 
"blocks".  A  "block"  is  a  set  of  instructions  between  a  transfer 
entry  point  and  the  next  transfer  exit  point.  Thus,  a  "block"  is 
completely  processed  if  its  first  member  is  executed. 

b.  The  Active  and  Passive  reference  tables  are  used  to  construct 
the  Constant  and  the  Result  tables. 


c.  The  Data  and  Storage  tables  are  used,  along  with  the  Constant 
and  Result  tables,  to  determine  if  the  input  program  changes 
or  modifies  any  of  its  own  instructions. 


C.  MICRO  DATA  COMPRESSION 

The  third  pass  of  the  analysis  program  processes  the  tables  of  the  second 
pass.  Examples  of  the  processing  arc: 

a.  The  Block  table  is  used  to  determine  if  all  non-data  blocks  can 
be  reached  from  either  the  starting  location  or  the  various 
entry  locations.  If  a  block  cannot  be  reached,  then  the  topology 
of  the  surrounding  blocks  is  studied  to  determine  what  con¬ 
nection  assumptions  should  be  made. 

b.  The  connected  Block  table  and  the  Active  and  Passive  tables 
are  then  used  to  construct  the  Latest  reference  table  as  a 
function  of  program  topology  for  each  passive  reference. 

D.  MACRO  DATA  COMPRESSION 

a.  Each  active  reference  represents  data  changes  that  are  made 
by  the  program  being  analyzed.  By  using  the  Latest  table  entry 
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for  each  Passive  table  entry,  a  symbolic  expression  for  each 
Active  table  entry  can  be  constructed  as  a  function  of  program 
topology. 

b.  Cross  reference  information  showing  Latest  values  of  used 
registers,  e.g.,  index  registers,  is  also  provided. 

The  final  flow-chart  output  can  be  set  to  show  either  a  topology-only  level 
or  topology-and-data-processing  level.  Further  refinement  of  these 
output  formats  is  necessary,  and  this  requires  further  consideration  of 
possible  data-flow  representations. 

Visual  Information  Processing  -  Herbert  M,  Teager 

Major  effort,  under  ONR  Contract  Nonr-184  1(69)  and  JSEP  Contract 
DA36-039-AMC-03200(E),  has  b  een  devoted  to  the  area  of  human  sensory 
processing,  primarily  visual.  The  object  of  this  effort  has  been  three¬ 
fold,  First,  from  a  pragmatic  viewpoint,  to  provide  for  handwritten  and 
drawn  symbolic  input  to  digital  systems  with  a  latitude  and  scope  approach¬ 
ing  normal  human  practice.  Second,  to  develop  a  testable  theory  of  the 
visual  processing  that  occurs  in  living' systems,  since  our  theories  with 
respect  to  pattern  recognition  are  generally  inadequate  for  almost  any  of 
the  commonly  encountered  facets  of  human  visual  capability,  such  as 
visual  memory,  depth  perception,  face  and  object  recognition,  et  al. 

Third,  to  understand  the  "what"  and  "how"  with  respect  to  sensory  pro¬ 
cessing,  which  is  a  fascinating  and  potentially  fruitful  field  of  scientific 
research  in  its  own  right. 

The  work  that  was  accomplished  has  been  presented  as  a  paper,  entitled 
Multidimensional  Visual  Information  Processing,  to  the  New  York  Academy 
of  Sciences.  There  is  not  the  space  available  to  duplicate  the  material  of 
that  paper  here.  However,  the  results  with  respect  to  development  of  a 
useful  system  are  summarized  below. 

It  is  possible  to  build  a  simple  hardware  system  to  accomplish  real-time 
character  recognition,  when  connected  to  a  "Teager  table/Rand  tablet", 
for  a  known  individual's  writing  with  the  following  constraints  and 
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capabilities:  Printed-writing  script  must  be  disconnected  and  at  the 
general  level  of  complexity  of  Roman  character  fonts,  with  no  limitation 
on  writing  speed.  The  system  is  invariant  to  size,  location,  rotation 
skew,  and  all  other  affine  transformations  generated  by  oblique  projec¬ 
tion.  Within  the  limits  of  character  complexity  and  delay-line  storage 
capacity,  the  system  is  unlimited  with  respect  to  character  font.  In 
practice,  a  font  of  no  more  than  256  different  characters  is  optimal  with 
respect  to  the  hardware.  These  conclusions  have  been  verified,  with 
simulated  hardware,  by  an  extended  period  of  experimentation  on  writing 
from  a  varied  group  of  subjects,  during  which  the  very  low  error  rate 
that  occurred  (less  than  3%)  could  be  traced  almost  exclusively  to  the  use 
of  symbols  that  were  completely  new  to  the  system. 

Further  work  is  being  expended  to  build  the  hardware,  as  well  as  extend 
its  capability  to  non- real -time  recognition,  unknown  source,  connected 
writing,  and  a  character  font  at  the  level  of  comple.xity  of  Chinese  and 
shorthand. 

With  regard  to  the  development  of  a  tentative  theory  for  visual  processing 
of  two-dimensional  information,  conclusions  have  been  reached  with 
respect  to  the  form  of  measurement  and  processing  which  may  well  be 
occurring  in  living  systems  if  the  known  experimental,  psychological,  and 
artistic  evidence  of  line,  shape,  and  form  incoding  is  considered.  These 
measurements,  in  fact,  form  a  basis  for  the  processing  system  that  has 
been  developed.  Evidence  for  such  measurements  and  analogous  ones  in 
other  sensory  modalities  are  being  sought,  and  consideration  is  being 
given  to  means  of  artificially  accomplishing  such  processing  with  con¬ 
temporary  processing  and  sensing  hardware. 

One  tentative  theoretical  conclusion  with  respect  to  visual  processing  is 
that  optic  nerve  connections  are  far  too  few  in  number  to  be  carrying  in¬ 
formation  to  the  brain  in  the  assumed  one-way,  binary  fashion.  Obviously, 
analogous  speculations  are  for  future  exploration  in  conjunction  with 
appropriate  life  and  medical  scientists. 

At  the  operational  level,  a  multiplexing  system  of  hard-copy  graphical  out¬ 
puts  and  its  attendant  programming  was  turned  over  to  the  ESL  group;  and 
Eric  Westerfeld  developed  a  program  to  translate  hand-drawn  flow  charts 
into  machine-compilable  MAD  programs.  (See  Westerfeld,  Appendix  B) 
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CTSS  Operation  and  Equipment  -  M.  V.  Solomita 


During  this  reporting  period  the  IBM  7094  was  scheduled  for  operation 
twenty  four  hours  a  day,  seven  days  a  week  with  the  exception  of  a  shut- 
dov,n  for  the  Christmas  and  New  Years  Holidays.  The  computer  was  also 
not  in  operation  during  the  installation  of  the  IBM  7289  channel  and  7320A 
high  speed  drums  in  September,  and  the  replacing  of  the  IBM  1301  Disk 
storage  file  with  the  IBM  2302  Disk  storage  file  in  January.  These  changes 
have  helped  increase  the  speed  and  storage  capacity  of  the  Compatible 
Time-Sharing  System. 

In  June  of  1965  the  GE  Model  635  began  arriving  and  by  June  15th  the 
installation  of  the  equipment  began.  Wo  expect  operation  of  the  computer 
to  start  during  the  third  quarter  of  1965. 


The  IBM  7094  Data  Processing  System  configuration  now  consists  of  the 
following : 


Type 

7109 

7110 

2  7302 
7151 

3  7607 

3  761  7 

2  7909 

7606 

7608 

7618 

7631 

2302 

7320 

7750 

2  7320A 

12  729VI 

716 
711 
721 
7289 


Description 

Control  Processing  Unit 
Control  Processing  Unit 

Magnetic  Core  Storage  (32,  768  words  each) 

Console  Control  Unit 

Data  Channels 

Data  Channel  Consoles 

Data  Channels 

Multiplexor 

Power  Converter 

Power  Control 

File  Control  Unit 

Disk  Storage  Unit  (234  million  characters) 

Drum  Storage  Unit  (1,1  18,400  characters) 

Data  Communication  Channel 

Drum  Storage  Units  (1,118,400  characters  each) 

Magnetic  Tape  Units 

Alphabetic  Printer 

Punched  Card  Reader 

Card  Punch 

Data  Channel 
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Peripheral  (Off-Line)  Equipment: 


Type 

1401 

Central  Processing  Unit  (4,  000  characters 

1402 

Card  Read-Punch  . 

2 

1402 

Printers  (132  character/line  each) 

4 

729-V 

Magnetic  Tape  Units 
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Research  on  the  Compatible  Time-Sharing  System 
MULTICS'  Hardware  System  Design 
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Research  on  the  Coinijatible  Tinge -Sharing  System  -  Fei’nando  J.  Corbalo 

'■J'he  system  programming  group  has  continued  developing  the  Compatible 
1  inie-Sharing  System  (CTSS)  for  tli_e  IBM  709*1  computer  of  Project  MAC. 
This  evolution  has  allowed  the  constant  modification  and  testing  of  new 
idttas  in  a  framework  of  man-machine  interaction.  In  fact,  the  ability  c  ' 
the  system  to  evolve  with  time  has  been  a  basic  design  goal.  Hardware 
1  ecjui  1  ements  in  the  past  year  have  also  caused  further  progi'dlmniing  re¬ 
quirements.  Two  IBM  7i20A  high-speed  drums,  attached  to  the  MAC 
computer  in  August,  have  significantly  improved  response  times. 

Ihus,  time-sharing  at  Project  MAC  has  been:  1)  a  service  facility  which 
has  allowed  individual  users  to  e.xplore  the  consequences  of  man-machine 
interaction,  and  2)  a  laboratory  in  which  system  programmers  can  develop 
the  tools  needed  for  effective  man-machine  interaction.  The  7094  system 
is  still  in  the  process  of  being  tuned,  though  it  is  incapable  of  being  tuned 
to  any  high  degree  because  the  basic  multi -programming  ability  and 
storage  allocation  techniques  are  very  difficult.  For  this  reason,  a  second 
genet  ation  time-sharing  system  called  MULTICS  (for  Multiplexed  Informa¬ 
tion  and  Computing  Service)  will  be  implemented  initially  on  a  GE  645  as  a 
development  project  with  the  Bell  Telephone  Laboratories  and  the  General 
Electric  Computer  Department. 

The  Multics  Project  was  started  this  year,  and  papers  describing  the 
system  are  to  be  presented  at  the  Fall  Joint  Computer  Conference,  1965. 
The  major  design  decisions  which  have  been  made  are  in  the  area  of 'hard¬ 
ware.  A  separate  report  by  E.  L.  Glaser  covers  hardware  design.  Soft¬ 
ware  plans  have  been  harder  to  firmly  establish. 

4  uning  anc  improvement  of  the  7094  systenr  has  been  done  w'ith  two  motives 
there  has  been  the  desire  to  develop  an  accurate  pilot  model  of  possible 
future  systems,  and  there  has  been  great  interest  to  learn  from  mistakes. 
On  both  counts  it  is  felt  that  considerable  success  has  been  achieved  during 
the  past  year.  A  detailed  description  of  the  7094  system  may  be  found  in 
the  Compatible  Time-.Sharing  System,  2nd  edition,  M.I.  T.  Press,  1965. 
What  follows  is  a  brief  review  of  the  major  elements  of  the  evolvement  of 
the  7094  system. 
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It  has  been  recognized  that  the  disk  file  system  is  a  major  and  central 
part  of  a  time-sharing  system.  Already  CTSS  has  one  of  the  most 
elaborate  file  systems  developed.  Nevertheless,  it  is  also  recognized 
that  the  first  version  of  the  file  system  is  inadequate  for  a  more  highly- 
tuned  time-sharing  system.  Some  of  the  changes  that  have  been  improve¬ 
ments  are  the  salvager  and  an  improved  input/output  editor  for  placing  or 
removing  material  from  the  disk  file.  These  changes  have  improved  the 
reliability  and  smoothness  of  system  operation  inasmuch  aq  the  features 
involved  occupy  pivotal  software  roles  and  are  expected  to  always  function 
without  mishap. 

However,  of  major  importance  has  been  the  design  of  the  second-generation 
file  system.  This  major  redesign  is  expected  to  correct  four  areas  of 
difficulty.  One,  it  should  be  possible  to  share  files  among  many  users, 
thereby  allowing  multiple  reading  and  suitable  interlocking  upon  writing. 
Two,  multiprogramming  should  be  possible  such  that  more  than  one  user 
can  be  using  the  file  system  simultaneously.  Three,  it  should  be  possible 
to  use  magnetic  tapes  from  within  the  time-sharing  system.  Four,  it 
should  be  possible  to  incrementally  dump  only  newly-created  material,  on 
a  continuous  basis,  rather  than  have  a  massive  dump  periodically  (e.  g., 
once  a  day)  of  the  entire  contents  of  the  disk.  This  process  is  being 
developed  under  the  program  name  DAEMON.  Basic  design  was  concluded 
in  the  fall  of  1964,  Implementation  was  begun  and  is  now  approaching  the 
terminal  stages.  The  implementation  of  this  new  file  system  is  considered 
of  prime  importance,  as  it  gives  the  programming  staff  an  opportunity  to 
shake -down  the  multitude  of  ideas  which  are  involved  in  anticipation  of 
■MULTI  CS. 


Other  improvements  to  CTSS  have  been  the  introduction  of  a  class  of  new 
editing  programs  for  typing  in  programs,  text,  etc.  '  Of  notable  interest 
has  been  the  development  of  the  TYPSET  and  RUNOFF  commands  for 
English  text  and  the  corresponding  ED  and  EDL  commands  for  program 
text.  These  programs  are  based  on  a  synthesis  of  many  editing  techniques 
by  J.  Saltzer  and  have  mostly  replaced  the  need  for  the  previous  editors 
such  as  MEMO,  MODIFY,  DITTO,  INPUT,  EDIT,  and  FILE.  The  new 
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editors  are  distinguished  by  a  high  degree  of  sophistication  for  man- 
machine  interaction  as  well  as  context  editing  techniques  which  avoid 
fixed  format  requirements. 

In  addition,  nearly  all  commands  of  the  system  have  had  some  corrections 
or  improvements.  Particularly  notable  are:  inclusion  of  a  new  version  of 
the  MAD  compiler;  cleaning  up  and  fixing  up  of  the  FAP  command  to  allow 
use  of  more  compacted  input  files:  improvement  of  the  LOAD  command  to 
allow  three  different  program  libraries  and  thus  speed  up  library  search¬ 
ing:  and  reprogramming  of  FILE.  RUNCOM,  SPLIT,  COMBIN,  SNOBAL 
and  other  commands.  Further  new  features  have  been  added  in  the  form 
of  interconsole  message  facilities  and  the  ability  to  "attach"  other  consoles. 
In  addition,  new  debugging  techniques  have  been  introduced,  the  most 
notable  being  the  MADBUG  program  of  R,  Fabry  for  the  FAPDBG-like 
debugging  of  MADBUG  programs.  The  MADBUG  program  is  considered 
only  a  prototype  of  future  debugging  techniques,  as  the  present  CTSS  im¬ 
plementation  forces  it  to  be  relatively  clumsy.  Nevertheless,  the  design 
goal,  that  the  MAD  language  user  need  not  know  machine  language,  has 
largely  been  met. 

Finally,  there  have  been  numerous  attempts  to  find  solutions  for  the  diffi¬ 
cult  problem  area  of  system  documentation.  There  have  been  many  at¬ 
tempts  to  compact  disc  space  requirements  and  allow  larger  utilization  of 
the  available  secondary  storage  capacity.  In  particular,  the  ARCHIVE  and 
CRUNCH  commands  have  been  developed.  In  addition,  the  LOG  command 
has  been  developed  to  notify  system  users  and  allow  user  teams  to  keep 
posted  all  system  changes  of  current  interest.  Moreover,  the  system 
programmers  themselves  have  developed  conventions  and  procedures  in  an 
attempt  to  systematize  record  keeping  for  a  large  system.  In  addition, 
facilities  have  been  developed  for  semi-automatic  "carry"  of  programs  back 
and  forth  between  the  two  7094  computers  to  allow  smoother  and  easier 
maintenance  of  two  slightly  different  CTSS  systems. 


Of  major  importance  in  the  orderly  development  of  the  new  system  has  been 
a  complete  revision  of  the  present  CTSS  Programmer's  Guide  to  include 
information  from  many  CTSS  Bulletins  and  important  MAC  memos  des¬ 
cribing  changes  in  the  features  of  the  system.  This  revised  manual  will  be 
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e  generality  of  design  is  intended  to  make  modular  software  in  the  new 

firres”"' 

to  facil  ta"r  ^quip-ent  control  unit 

faciLtate  processing  of  many  interrupts  in  a  short  period  of  time,  which 

minair'""''  connected  to  a  large  number  of  remote  ter- 


.he  third  major  area  of  design  is  that  of  a  new  high-speed  drum  controller 

this  nin  "  ^  large-capacity  drum.  The  availability  of 

his  new  drum  controller  with  its  associated  drum  will  facilitate  time-sharin, 
wo  reasons:  first,  the  transfer  rate  between  the  drum  and  core 
memory  is  approximately  fifty  percent  of  memory  speed;  second,  because 
the  way  the  drum  controller  operates,  it  is  possible  on  a  statistical 
to  Ignore  any  latency  due  to  the  rotary  access  to  the  drum. 

Details  on  all  of  these  design  areas  will  be  available  in  papers  to  be 
presented  at  the  Fall  Joint  Computer  Conference,  1965, 
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Inti  Qcluction  -  J.  Francis  Reintjes 

Several  research  groups  of  the  Electronic  Systems  Laboratory  are  parti- 
cipating  in  Project  MAC.  Display  systems  research  has  been  conducted 
in  support  of  the  Laboratory  effort  on  computer-aided  design  and  for  the 
purpose  of  enhancing  the  output-display  capabilities  of  the  MAC  time- 
computer  system. 


The  Computer  Applications  Group  continued  its  research  in  fundamentals 
of  man-machine  problem-solving,  with  emphasis  upon  the  development  of 
generalized  computing  techniques  which  are  expected  to  be  applicable  to  a 
wide  variety  of  engineering  design  problems.  Emphasis  in  this  project  is 
on  the  use  of  computers  as  an  aid  in  mechanical  design. 

Computer-aided  electronic  circuit  design  was  also  investigated  under  a 
grant  from  the  National  Aeronautics  and  Space  Administration.  The  main 
efforts  of  the  group  were  devoted  to:  analysis  techniques  for  linear  and 
nonlinear  networks,  logic  design  and  synthesis  through  use  of  the  MAC 
computer,  design  of  threshold  logic,  and  the  development  of  a  low-cost 
teletype-operated  graphical  display. 

The  Project  MAC  time-sharing  system  has  a.' so  been  very  useful  fpr  in¬ 
vestigating  simulation  of  force-rebalance  accelerometer  output  data.  The 
objective  of  this  work  is  to  obtain  maximum  output  information  by  using 
specialized  data  processing.  In  another  aspect  of  the  same  research  pro¬ 
gram,  the  group  investigating  computations  peculiar  to  strapped-down  navi¬ 
gation  systems  found  the  Project  MAC  PDP-6  to  be  highly  useful. 

A  total  of  16  staff  members  and  12  students  in  the  Electronic  Systems 
Laboratory  have  benefited  in  their  research  through  the  use  of  Project 
MAC  facilities. 
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Display  Systems  Research  -  John  E.  Ward  and  Robert  H.  Stotz 
A.  ESL  DISPLAY  CONSOLE 

The  ESL  Display  Console  was  designed  to  satisfy  needs  of  the  M.I.T. 
Computer-Aided  Design  Project,  and  its  construction  was  supported  by 
the  Air  Force  Materials  Laboratory,  Wright-Patterson  Air  Force  Base 
under  Contract  AF-33(657)-10954.  The  unusual  features  of  this  display’ 
system,  which  have  permitted  its  routine  use  in  the  Project  MAC  time¬ 
sharing  system  since  January.  1964,  are:  its  ability  to  time  share  the 
7094  memory  through  the  Direct  Data  Channel;  its  incremental  digital 
line  generation,  scaling,  and  rotation  system;  and  its  automatic -hardware 
pen-tracking  system.  During  the  past  year,  the  system  capabilities  have 
been  extended,  with  Project  MAC  support,  by  adding  a  second  operator 
station  (consisting  of  a  slave  display  scope  and  associated  manual  inputs) 
which  time-shares  the  display  generation  electronics,  thus  allowing  two 
operators  to  work  simultaneously  on  different  problems.  Typical  current 
usages,  which  are  reported  in  other  sections  of  this  progress  report, 
include:  three-dimensional  drawing  programs,  parametric  surface  Studies 
and  portrayal  of  stress  contours  by  the  Computer-Aided  Design  Project; 
three-dimensional  models  of  protein  molecules  by  the  Biology  Department; 
modeling  of  the  human  vocal  tract  by  the  RLE  Speech  Group;  design  of 
highway  interchanges  by  the  Civil  Engineering  Department;  and  general 
graphical  output  by  a  number  of  other  groups  who  use  the  display  in  an 
interactive  way.  Usage  currently  averages  about  16  hours  per  day  for  each 
of  the  two  console  stations. 


The  M.I.T.  Computation  Center  has  recently  decided  to  install  a  complete 
duplicate  of  the  ESL  Console  on  its  time-shared  7094  to  extend  the  avail¬ 
ability  of  display  facilities  to  a  larger  segment  of  the  M.I.  T.  community. 
The  systems  will  be  compatible  and  Project  MAC  display  software  can  be 
used  directly.  Digital  Equipment  Corporation  is  building  the  electronic 
portions  of  this  second  system  from  M.I.T.  drawings,  and  the  manual  in¬ 
put  devices  are  being  fabricated  in  the  Electronic  Systems  Laboratory. 

This  second  console,  which  is  supported  by  an  NSF  Grant,  is  scheduled  for 
installation  in  November,  1965. 
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B.  STORAGE-TUBE  DISPLAYS 

Experiments  with  driving  "one-shot"  storage  tube  displays  from  analog 
deflection  voltages  of  the  ESL  Display  Console  have  been  quite  successful, 
although  the  quality  of  the  present  storage  tubes  leaves  something  to  be 
desired.  Two  such  displays,  connected  to  the  ESL  console  by  three 
coaxial  cables,  have  been  in  experimental  use  during  the  past  few  months. 
Programming  and  control  circuit  modifications  to  permit  the  addressing 
of  up  to  ten  remote  storage  scopes  for  transmission  of  "one  shot"  pictures 
liave  been  accomplished.  Figure  3  shows  typical  storage  displays.  These 
require  only  about  one  half  second  of  console  time  to  create,  but  are  re¬ 
tained  for  periods  up  to  one  hour,  unless  erased  manually  or  by  the  com¬ 
puter.  Experiments  are  continuing  with  the  use  of  analog  data  sets  on 
three  telephone  circuits  (for  X  and  Y  deflection  voltages  and  intensification 
pulses)  for  remote  operation  of  storage  scopes  in  this  mo.de'at  greater 
distances.  Noise  and  linearity  characteristics  of  the  telephone  circuits 
appear  satisfactory,  though  the  difference  in  time  delays  on  the  three 
circuits  is  proving  troublesome. 

1 

C.  DISPLAYS  FOR  THE  GE  645  COMPUTER 

Study  of  future  display  requirements  for  the  Project  MAC  GE  645  Computer 
has  resulted  in  identification  of  two  goals;  1)  a  modest  number  (10-15)  of 
remote  sophisticated  displays,  with  capabilities  similar  to  the  present  ESL 
Console;  and  2)  a  large  number  (perhaps  500)  of  low-cost,  limited- 
capability  displays  which  can  be  located  at  every  MAC  terminal. 

Planning  for  low-cost  remote  display  scopes  currently  centers  around  a 
digitally-controlled  unit  (rather  than  the  analog  transmission  system  dis¬ 
cussed  in  B,  previously),  with  line  and  character  generation  being  performed 
locally.  A  set  of  desirable  characteristics  for  a  unit  which  could  operate 
from  standard  2000 -bit-pe r -second  Dataphones  has  been  prepared  and  is 
being  discussed  with  the  manufacturers  of  display  equipment.  Also, 
research  work  is  being  conducted  on  circuits  for  use  in  such  a  device.  One 
master's  thesis  project  by  Mr.  T.  B.  Cheek  concerns  new  approaches  to 
character  generation,  with  the  objective  of  exploiting  low-speed  requirements 
to  achieve  very  low  cost.  A  low-cost,  hybrid  digital/analog,  line -generation 
system  is  also  being  designed  to  operate  from  a  Dataphone.  This  approach 
combines  a  binary  rate  multiplier  with  an  analog  integrator,  thus 
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eliminaling  the  need  for  digital -to -analog  converters.  A  prototype  unit, 
with  parts  costing  about  $1,000,  should  be  ready  for  initial  tests  by 
September,  1965. 

Operating  more  sophisticated  remote  displays  (like  the  ESL  Console) 
directly  from  the  central  computer  (in  the  manner  of  the  present  ESL 
Console)  will  present  definite  communications  problems,  unless  a  fle.xible, 
low-cost,  multi-megacycle,  two-way  link  becomes  available.  Also,  the 
load  on  the  central  computer  must  be  considered.  Our  e.xperience  indicates 
that  the  best  approach  is  to  use  a  small  satellite  computer  with  each  dis¬ 
play  (or  group  of  d’splays)  to  buffer  display  data  transmitted  over  medium- 
speed  telephone  circuits  (9  to  10.x  lO"^  bits  per  second).  Tiie  satellite 
computer  would  also  perform  real-time  display  processing,  such  as  pen 
tracking,  rotation  calculations,  and  interrupt  analysis.  One  such  com¬ 
puter-buffered  display  (the  DEC  338)  will  be  available  at  Project  MAC 
dining  the  coming  year  for  e.xperimental  evaluation  and  study  of  pro¬ 
gramming  problems.  (Note  that  the  small  computer  program  can  be 
considered  part  of  the  supervisor  program  of  the  central  machine.) 

\ 

D.  TV  SCAN  CONVERSION 

The  display  group  has  completed  the  design  for  an  image -maintaining 
display  system  based  on  a  TV  scan  converter  for  the  U.  S.  .Naval  Underwater 
Ordnance  Station  under  Contract  N140(  1  22) -76 1  18B .  Computer-generated 
pictures  are  drawn  on  the  storage  surface  of  t  scan -con ve rsion  tube  by 
slaving  its  deflection  inputs  to  an  e.xisting  DEC  Type  340  display  system. 

Tlien,  by  switching  the  deflection  irputs  to  a  TV  raster,  the  stored 
picture  may  be  viewed  for  periods  up  to  five  minutes  on  standard  TV 
monitors.  A  1203-line  scan  is  used  to  provide  high  resolution.  Selective 
erasure  or  writing  of  up  to  three  characters  will  be  possible  during  each 
IV  vertical  retrace  time,  so  that  stored  pictures  may  be  altered  and  updated 
without  requiring  a  complete  regeneration  by  the  computer.  This  system 
will  be  e.xperimentally  tested  on  Project  MAC's  PDP-6  computer  before 
being  delivered  to  the  Navy  later  in  the  year. 
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E.  IMPROVED  DISPLAY  TECHNOLOGY 

A  number  of  other  aspects  of  display  technology  are  being  studied  as 
thesis  projects.  Mr.  W.  D.  Stratton  is  working  on  a  new  high-speed 
pen-tracking  system,  based  on  the  amplitude  comparison  of  light-pen 
responses  to  displayed  points  in  a  tracking  pattern,  similar  to  radar 
conical-tracking  schemes.  By  replacing  the  present  digital  point-counting 
schemes,  it  is  hoped  that  pen-tracking  time  can  be  reduced  by  a  factor  of 
perhaps  20  to  1.  The  new  electrostatic  "beam  pen"  being  developed  by  the 
group  IS  being  used  in  this  work,  although  the  technique  could  also  be 
applied  to  light  pens. 

The  possibilities  of  infrared  (IR)  data  links  for  remote  operation  of  dis¬ 
plays  is  being  investigated  in  a  thesis  project  by  Mr.  George  Ling.  We 
have  obtained  an  e.xpe rimental  IR  link,  on  loan  from  the  Lincoln  Laboratory, 
that  is  capable  of  10  bits  per  second  over  a  two-mile  range,  and  have 
evaluated  its  performance  in  transmitting  digitally-coded  messages.  The 
link  is  now  being  set  up  to  transmit  display  information  from  the  ESL 
Console  in  Technology  Square  to  a  storage  scope  in  the  Electronic  Systems 
Laboratory,  a  distance  of  about  600  yards. 

As  part  of  the  continuing  program  to  improve  the  ESL  Display  Console.  Mr. 
Endre  Guttman  completed  his  thesis  study  of  a  digital  differential-analyzer 
(DDA)  rotation  matrix  to  replace  the  present  binary  rate -multiplier  (BRM) 
rotation  matrix.  In  the  BRM  matrix,  shown  in  Figure  4.  the  10-bit  signed 
numbers  ij^,  i^,  etc.,  are  set  by  the  computer  and  are  normally  the  vector 
cosines  specifying  rotation  between  the  x,  y,  z  and  h.  v.  d  coordinate  sys¬ 
tems  (depth  coordinate,  d,  is  not  computed).  Scaling  of  these  numbers 
controls  the  size  of  the  displayed  image.  The  boxes  labeled  "increment 
logic"  combine  and  buffer  the  BRM  outputs,  releasing  them  to  the  scope  de¬ 
flection  registers  in  binary  weights  of  1.  2.  4.  or  8  steps,  which  the  regis¬ 
ters  (up-down  counters)  can  accept. 


Although  the  BRM  matrix  is  generally  satisfactory,  it  has  significant 
round-off  errors  which  create  picture  distortions.  Also,  the  buffering 
action  of  the  increment  logic  adds  some  roughness  to  displayed  lines. 
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Figure  5  shows  the  proposed  matrix,  which  would  have  two  identical  three- 
input  parallel  DDA's,  one  for  h.  and  one  for  v.  The  advantage  of  the  DDA 
over  the  BRM  is  that  it  has  a  remainder  register  to  "remember"  round-off 
in  a  particular  step  of  the  calculation,  and  apply  it  to  the  subsequent  step. 
Also,  the  h  and  v  outputs  would  be  obtained  directly,  eliminating  part  of 
the  function  of  the  increment  logic.  The  increment  logic  could  then  be 
eliminated  completely  by  converting  the  first  few  stages  of  the  deflection- 
registers  from  counters  to  adders. 

A  DDA  for  this  purpose,  however,  must  be  able  to  add  four,  ten-bit.  signed 
numbers  in  a  total  elapsed  time  of  only  1.5  microsecond.  It  was  s^hown  that 
multi-input  threshold-logic  elements  offer  the  only  feasible  way  of  meeting 
such  a  severe  speed  requirement  at  reasonable  cost.  A  potentially  suitable 
six-input  threshold-adder  stage  (four  register  inputs,  plus  two  carries) 
requiring  only  three  transistors  was  designed  and  breadboarded.  It  was 
found  to  have  a  transition  time  of  80  nanoseconds;  fast  enough  for  this 
purpose.  Further  development  and  testing  of  threshold  circuits  will  be 
necessary,  however,  before  the  DDA  rotation  matrix  can  be  constructed. 

F.  FUTURE  PLANS 

All  reported  activities  will  continue  during  the  coming  year,  with  particular 
emphasis  on  developing  the  required  display  technology  for  future  graphical 
man-machine  communication  and  studying  the  hardware/software  problems 
in  interfacing  large  numbers  of  displays  with  the  next  generation  time¬ 
sharing  computers. 
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Computer-Aided  Design  -  Douglas  T.  Ross  and  Clarence  G.  Feldmann 

The  Computer  Applications  Group,  along  with  faculty  and  staff  of  the 
Mechanical  Engineering  Department,  has  continued  its  Computer-Aided 
Design  work,  with  emphasis  on  problems  of  mechanical  design.  This 
program  is  being  sponsored  by  the  Air  Force  Materials  Laboratory, 
Wright-Patterson  Air  Force  Base,  under  Contract  AF-33(657)- 10954.  The 
work  IS  an  evolution  of  earlier  work,  conducted  at  ESL,  to  develop  the 
Automatically  Programmed  Tool  (APT)  System  now  being  used  by  industry 
to  program  numerically-controlled  machine  tools.  The  approach  being 
taken  to  computer-aided  design  is  sufficiently  general  to  make  it  applicable 
to  a  wide  variety  of  design  problems  in  engineering  and  science. 

The  current  research  effort  has  four  phases.  The  first  phase  involves 
research  in  fundamentals  of  man-machine  problem-solving,  the  structure 
of  language  and  language  processing,  the  modeling  of  problems,  and  the 
development  of  generalized  computing  techniques.  The  second  phase 
concerns  application  of  the  theoretical  foundations  resulting  from  this 
research  to  the  design  and  construction  of  a  generalized  man-machine 
problem-solving  system,  including  development  of  a  powerful  language  for 
writing,  operating,  and  debugging  programs,  and  facilities  for  efficient 
construction  of  specialized  problem-oriented  systems.  The  third  phase 
is  actual  development  of  such  systems  for  specific  design  areas,  and  the 
remaining  phase  is  adaptation  and  application  of  special  systems  by 
designers  to  meet  their  individual  needs. 


In  the  realm  of  generalized  languages  and  systems,  the  AED-0  system 
(standing  for  Automated  Engineering  Design)  was  the  first  system  completed 
by  the  group.  It  has  been  in  operation  for  over  a  year  and  is  being  employed 
extensively  by  members  of  the  Design  Group  and  by  other  groups  at  Project 
MAC.  More  recently,  through  the  efforts  of  Mr.  Feldmann  and  his  associates, 
AED-0  has  been  made  available  as  a  batch-processing  system  and  distributed 
to  interested  industrial  users.  Mr.  Ross  presented  a  series  of  lectures  on 
the  salient  features  of  AED-0  language,  and  his  notes  now  constitute  the 
basic  AED-0  Programmers'  Guide, 


1 


ELECTRONIC  SYSTEMS  LABORATORY 


81 


During  the  year,  many  .substantial  additions  were  made  to  AED-0  in  order 
to  increase  its  power  for  writing  advanced  programming  systems.  The 
language  is  based  on  Algol-60,  with  many  additions  and  modilications  to 
perrhit  construction  and  manipulation  of  complex  data  structures,  includ¬ 
ing  the  ability  to  reference  subparts  of  computer  words  for  efficient  use  of 
storage.  The  system  also  includes  an  elaborate  Macro  Preprocessor  which 
permits  extensive  manipulation  of  the  input  string  of  characters  before 
compilation,  thus  providing  the  user  with  greatly  enhanced  expressive 
capabilities.  A  further  feature  is  that  the  Preprocessor  includes  facilities 
for  regenerating  AED-0  language  programs  in  a  meaningful  format,  which 
exhibits  the  logical  structure  of  the  program.  Combined  with  flexible 
"comment"  and  "remark"  features,  this  capability  makes  AED-0  programs 
almost  self-documenting,  minimizing  the  need  for  elaborate  flow  diagrams 
to  describe  programs.  A  source-language  debugging  system  has  been 
written,  so  that  the  user  may  trouble-shoot  programs  at  the  AED-0  lan¬ 
guage  level,  with  little  concern  for  the  actual  machine  code  generated  by  the 
compiler.  AED-0  development  is  now  complete,  and  primary  attention  is 
now  focused  on  its  use  as  a  tool  to  achieve  even  more  sophisticated  gen¬ 
eralized  programming  and  compilation  systems. 

The  immediate  successor  to  AED-0  will  be  the  AED-1  System  now  under 
development.  A  salient  feature  of  AED-1  is  that  it  will  be  highly  machine- 
independent.  That  is,  it  is  being  designed  so  that  both  the  language  and 
techniques  used  to  compile  machine-code  programs  will  be  compatible 
with  virtually  any  type  of  large-scale  computer.  Several  valuable  charac¬ 
teristics  of  AED-0  will  be  contained  in  its  successor,  as  well  as  new 
features  to  give  added  flexibility  and  capability.  The  AED-1  System  will 
be  able  to  process  AED-0  as  well  as  AED-1  language  programs.  Several 
building -block  (subsystem)  programs  are  being  written  to  achieve  this  goal, 
including  a  subsystem  (AED  Jr.)  which  permits  the  user  to  define  and  check 
out  descriptions  of  arbitrary  programming  languages. 

With  respect  to  specialized,  problem-oriented  programming  systems,  the 
group's  first  major  system  in  this  category  is  now  being  prepared.  Called 
CADET  (Computer-Aided  Design  Experimental  Translator),  it  permits 
users  to  carry  out  computer-aided  design  applications  using  verbal  or 
graphical  language  or  both.  That  is,  its  purpose  is  to  enable  construction 
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and  manipulation  of  graphical  and  analytic  structures.  CADET  will  also 
permit  analytic  functions  and  physical  properties  to  be  associated  with 
graphical  representations  for  things.  Since  programs  defining  behavioral 
properties  of  graphical  elements  may  take  arbitrary  form,  they  are  written 
in  standard  AED  language. 

To  assist  the  initial  development,  of  CADET,  which  will  be  a  general- 
purpose  problem-oriented  system,  several  specific  design  problems  are 
being  explored  through  its  use.  One  of  these  is  concerned  with  the  general  ’ 
problem  of  three-dimensional  shape  description,  and  will  use  the  generalized 
parametric  surface  patches  developed  by  Professor  Coons  of  the  Mechani¬ 
cal  r-nginee  ring  Department,  as  well  as  the  three-dimensional  pseudo -pen 
program  written  by  Mr.  Polansky,  and  the  routines  prepared  by  Mr.  Lang 
for  ESL's  Graphical  Display  Console.  Another  application  is  in  the  area 
of  linear  and  nonlinear  electronic  circuits,  described  in  the  next  section. 

A  further  aspect  of  the  entire  program  in  Computer-Aided  Design  is  the 
woik  being  performed  on  cathode -ray-tube  displays  for  graphical  mani¬ 
pulation  of  design  information  and  symbols.  This  work  was  described  in 
the  previous  section. 


An  important  goal  of  the  Computer-Aided  Design  Project  is  the  dissemina¬ 
tion  of  information  generated  by  the  group  to  interested  users  outside  M.I.  T. 
It  was  our  pleasure  to  have  seven  system  programmers  from  industry  as 
guests  of  ESL  and  Project  MAC  for  the  past  year.  Our  guests  have  parti¬ 
cipated  in  developing  the  programming  systems  described  previously,  and 
they  plan  to  apply  and  further  expand  these  basic  concepts  upon  return  to 
their  respective  organizations.  This  program  of  cooperative  research  and 
development  with  visiting  staff  from  industry  will  continue  as  an  important 
aspect  of  ESL's  work  in  computer-aided  design. 
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Computer-Aided  Electronic  Circuit  Design  -  Michael  L.  Dertouzos  and 

J.  Francis  Reintjes 


A.  SIMULATION  OF  ELECTRICAL  NETWORKS 

A  method  is  being  investigated  for  on-line  simulation  of  electrical  net¬ 
works.  This  method  treats  networks  consisting  of  linear  and  nonlinear 
resistors,  as  well  as  inductors  and  capacitors.  A  given  network  can  be 
excited  with  voltage  and  current  sources  which  may  be  arbitrary  functions 
of  time.  The  voltage  response  across  any  node  in  the  network  may  then 
be  called  for  as  output. 

The  computer  model  for  the  network  consists  of  lists  of  storage  registers 
representing  branches  connected  to  other  lists  representing  nodes.  The 
method  for  solution  at  each  instant  of  time  is  as  follows:  The  nodes  are 
assigned  fixed  node  potentials.  The  voltage  across  each  branch,  and  hence 
the  current  through  it,  is  determined  by  the  two  nodes  to  which  it  is  con¬ 
nected.  Equilibrium  requires  the  sum  of  the  branch  currents  at  each  node 
to  be  zero.  If  the  sum  is  not  zero,  new  node  potentials  closer  to  the 
equilibrium  values  are  computed  by  perturbing  the  old  ones  in  a  prescribed 
manner.  This  process  is  repeated  until  the  equilibrium  values  are  obtained 
within  a  specified  tolerance. 

Thus  far,  a  computer  program  for  use  on  the  Project  MAC  time-sharing 
system  has  been  written  to  apply  the  foregoing  method  to  linear  RLC  net¬ 
works  excited  with  step  and  sinusoidal  current  and  voltage  sources.  The 
program  computes  an  accurate  solution  to  any  given  linear  network.  Com¬ 
puter  time  for  solution  of  networks  containing  up  to  about  8  nodes  and  12 
branches  is  on  the  order  of  a  few  seconds.  An  example  of  a  circuit 
analyzed  by  the  foregoing  program  is  shown  in  Figure  6. 

B.  DISPLAY  PROGRAM  FOR  ELECTRONIC  CIRCUITS 

A  computer  program  has  been  written  which  allows  an  operator  to  quickly 
and  easily  portray  planar  electronic  circuits  on  the  Electronic  Systems 
Laboratory  Console  Display  at  Project  MAC.  Circuit  elements  are  formed 
by  moving  a  tracking  light  pen  to  the  desired  position  on  the  screen  and 
actuating  an  "element  button."  Circuit  characteristics  are  described  in  a 
list  structure  with  pointers.  Used  in  conjunction  with  the  analysis  program 
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All  values  in  ohms, 
henrys,  and  farads 

(  a  )  Twin-T  Filter 


VOLTS 


SECONDS 

(  b  )  voltage  v(t)  in  Response  to  Unit  Step 


Figure  6.  Response  of  a  Twin-T  Filter  as  Computed  and 
Plotted  by  the  Circuit-Analysis  Program 
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described  previously,  circuits  with  multiple  elements  can  be  displayed  and 
leir  response  to  various  types  of  excitation  readily  determined. 

C.  ON-LINE  SWITCHING-FUNCTION  SYNTHESIS 

The  central  ain,  logical  design  is  the  synthesis  of  any  given  switching 
Uhctton,  t„  terms  of  given  sets  of  elementary  building  blocks,  for  the 
optimtaation  of  some  performance  index  in  the  presence  of  constraints. 
Although  the  present  state-of-the-art  yield,  algorithmic  method,  for  the 
so  ution  of  certain  specific  instance,  of  the  problem  (such  as  minimisation 

U1  “e-Whsk  input,  with  a  two-level,  AND-OR  realisation),  no  fully- 

a  gortthm.c  method  exist,  for  the  solution  of  the  more  general  problem. 

Solutions  to  the  general  problem  are  being  developed  through  use  of  an  on- 
ine  process,  where  the  machine  accomplishes  those  computational  tasks 
witch  can  be  algorithmically  specified,  and  the  user  provide,  those  de¬ 
cisions  which  he  is  belter  gualified  to  make.  The  machine  portion  of  the 
.ystem  IS  based  on  a  set  of  heuristic  procedures  which,  subject  to  certain 
onditions,  guarantee  convergence  of  the  process.  Result,  obtained  to 
date  indicate  that  this  method  o,  synthesis  yields  more  economical  struc- 
ure.  than  those  uiing  succe.sive  loc.l-optimiaation  procedures,  and  does 
not  depend  on  impkctical  (and  usifally  impossible)  exhaustive  searches 
through  all  possible  solution.,  I,  i.  expected  that  a  report  on  this  work 
Will  be  issued  during  the  Fall,  1965. 

D.  THRESHOLD  ELEMENT  REALIZABILITY 

Previous  work  has  shown  that  the  problem  of  determining  whether  an 
arbitrary  Boolean  function  of  N  variable,  is  realisable  with  a  single  threshold 
e  ement  can  be  reduced  to  the  problem  of  minimising  a  quantity  which  is  a 
unction  of  N+1  variables.  Sevefal  properties  of  this  quantity  are  being 
derived  and  studied.  It  can  be  shown,  for  example,  that  the  quantity  may 
be  interpreted  geometrically  a.  a  .tructure  of  intersecting  N+1  dimen¬ 
sional  hyperplanos.  Each  of  these  hyperplanes  corresponds  to  some  N- 
variable  Boolean  function,  and  if  the  given  Boolean  function  is  realisable 
tt.  associated  hyperplane  will  be  'Tlat"  with  sero  height.  The  hyperplane, 
are  so  situated  that  they  form  the  boundary  of  a  convex  body  and  hence  the 
quantity  contain,  no  "pocket."  of  local  minima.  This  property  enables  a 
computer  to  employ  hill -descending  techniques  in  N+1  space  which  seek  a 
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local  minimum  point.  Since,  by  virtue  of  the  convex  nature  of  the  quantity, 
any  point  satisfying  local  minimum  conditions  must  also  satisfy  absolute 
minimum  conditions,  hill-descending  techniques  will  lead  to  the  required 
solution. 

Several  hill-descending  techniques  for  performing  a  minimization  process 
are  being  studied  and  are  yielding  dependable  results.  The  question  of  the 
complete  generality  of  these  methods,  however,  is  still  an  open  one. 

E.  A  LOW-COST  PLOTTING  BOARD 

A  hard-copy  graphical-display  device  that  can  be  operated  from  low- 
speed,  teletype  data  sources  such  as  those  found  in  time-shared  computer 
installations,  is  being  investigated.  A  novel  feature  of  this  device,  shown 
in  Figure  7,  is  the  conversion  of  incremental  binary  data  into  graphical 
form.  This  conversion  is  primarily  accomplished  by  the  use  of  tuned 
dynamic  networks  at  each  axis.  Inclusion  of  mechanical  networks,  use  of 
certain  specific  torquers  for  actuators,  and  over-all  system  organization 
were  motivated  by  a  desire  to  keep  cost  at  an  absolute  minimum.  This 
objective  has  not  adversely  influenced  performance  or  speed  of  operation. 

The  latter  is  fundamentally  determined  by  low  teletype  information  rates 
(100  bits  per  second)  and  by  upper  bounds  on  the  torque-to-inertia  ratios 
of  reasonably-sized  actuators.  Theoretical  e.xtensions  of  the  foregoing 
basic  principle  have  shown  that  it  will  be  possible  to  implement  curvature 
control,  data  smoothing,  high-bandwidth  plotting  and  other  desirable 
prope  rties. 

F.  SIMULATION  OF  NONLINEAR  CIRCUITS 

A  different  attack  on  comfiuter  simulation  of  electrical  networks,  from  that 
described  in  Section  A.,  previously,  is  being  taken  in  this  investigation. 
Primary  emphasis  is  on  nonlinear-network  simulation,  and  network  response 
is  obtained  through  successive  solutions  and  modifications  of  a  matrix  which 
is  a  mathematical  description  of  the  network  being  simulated.  A  computer 
program  is  now  being  written  in  the  AED-0  lairguage  to  establish  the  net¬ 
work  equations  and  solve  them.  Later  work  will  include  programming  of 
the  ESL  Display  Console  for  graphically  portraying  network  diagrams  and 
responses  to  various  input  signals.  The  system  will  be  able  to  handle  time- 
dependent  sources  and  RLC  elements  whose  characteristics  are  piecewise 
linear  and  can  be  either  monotonic  or  nonmonotonic. 
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Figure  7.  A  Low-Cost  Teletype-Operated  Plotting  Board 
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Accelerometer  System  Studies  -  Alfred  K.  Susskind  and  Donald  R.  Haring 

A  novel  approach  to  obtaining  digital  data  from  force-rebalance  accelero¬ 
meters  is  under  investigation  for  the  U.S.  Air  Force  under  Contract 
AF-33{657)-l  131 1.  Because  the  proposed  system  contains  both  sampling 
and  quantization,  it  cannot  be  treated  adequately  by  analytical  techniques. 
Therefore,  experinientation  is  required,  which  can  be  economically  done 
only  by  means  of  simulation.  The  Project  MAC  Time-Sharing  System  has 
been  found  particularly  well  suited  for  this  purpose,  because  of  terminal 
equipment  such  as  the  ESL  Display  Console. 

By  interplay  between  theory  and  simulation,  good  progress  has  been  made 
in  understanding  and  evaluating  the  proposed  scheme.  The  observed 
system  performance  is  encouraging  and  the  scheme  will  therefore  be 
investigated  in  greater  depth,  gradually  removing  more  and  more  of  the 
idealizations  now  being  made  and  also  observing  the  effect  of  noise. 
Eventually,  a  rather  complex  dynamic  system  will  be  simulated,  and  it 
is  expected  that  sufficient  information  will  be  obtained  to  act  as  sound  guide 
lines  in  the  design  of  an  actual  accelerometer  system. 
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Simulation  of  Strapped-Down  Naviation  Systems  -  F,  B.  Hills  and 

J.  A.  Parisot 

In  a  research  program  for  the  Air  Force  under  Contract  AF-3  3(657)- 1  1  3 1 1 , 
we  are  studying  the  computations  peculiar  to  "strapped  down"  navigation 
systems,  in  which  inertial  sensors  (gyros  and  accelerometers)  are 
rigidly  mounted  to  the  airframe  rather  than  being  mounted  on  a  stabilized 
platform.  The  accelerometer  outputs  must  therefore  be  transformed  from 
body  axes  to  the  coordinates  in  which  the  navigation  computations  are 
made.  We  are  using  the  direction  cosine  transformation.  Because  of 
vehicle  rotations,  the  direction  cosines  are  not  fixed  and  hence  must  also 
be  computed.  This  is  done  by  numerically  solving  nine  simultaneous 
differential  equations  in  which  gyro  data  are  the  independent  variables. 

Analytical  studies  made  by  the  group  have  shown  that  the  major  sources 
of  error,  insofar  as  the  computations  are  concerned,  are  due  to  the  un¬ 
certainties  that  arise  in  vehicle  position  and  attitude  because  the  input  data 
must  be  sampled  and  quantized.  The  uncertainties  arise  primarily  be¬ 
cause  acceleration  and  rotation  are  not  communtative  and  the  components 
of  rotation  are  not  commutative.  Hence,  not  only  is  the  intersample 
behavior  of  the  input  data  important,  but  even  more  important  is  the 
relative  intersample  behavior  between  the  various  input  variables.  Some 
analytical  studies  into  the  size  of  the  uncertainty  have  been  made;  however, 
the  investigation  becomes  intractable  for  all  but  simple  flight  paths.  There¬ 
fore,  analysis  must  be  augmented  by  simulation. 

Project  mac's  PDP-6  was  chosen  for  these  simulation  studies  because  of 
its  ability  to  perform  double-precision  arithmetic  at  high  speed  (46-bit 
accuracy  is  required),  and  its  suitability  for  on-line  operational  control 
(e.g.,  possibility  of  connecting  external  equipment,  such  as  a  stick  to 
control  the  flight  of  a  simulated  vehicle).  Also,  macro  definitions  allowed 
by  the  PDP-6  assembly  program  make  it  possible  to  set  up  desired  tests 
with  a  minimum  of  debugging. 

A  diagram  of  the  simulation  is  shown  in  Figure  8.  In  the  box  labeled  vehicle 
simulation,  the  input  data,  acceleration  (a)  and  rotation  rate  (to),  are  gen¬ 
erated.  At  present,  these  quantities  are  generated  in  accordance  with  the 
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geometry  of  simple  flight  paths.  In  the  future,  a  simple  dynamic  model 
of  the  vehicle  with  manual  control  may  be  used  to  allow  a  search  for 
maximum  uncertainty. 

Two  computations  of  the  coordinate  conversion  computation  must  be  made: 
one  for  the  system  under  test,  the  other  for  a  reference  computation  which 
is  considered  to  be  the  true  solution;  The  reference  computation,  like 
that  of  the  system  computation,  requires  sampling  and  quantizing  of  the 
input  data.  Analysis  has  shown  that  all  numerical  formulas  compute  those 
solutions  to  the  direction  cosines  that  would  result  if  the  components  of  a 
were  proportional  between  samples.  Therefore,  to  make  the  reference 
computation  more  accurate  than  system  computation,  the  samples  for  the 
reference  must  be  taken  more  often  than  those  taken  for  the  system.  This 
produces  a  reference  solution  based  on  information  in  the  input  data  not 
available  to  the  system  computation,  and  hence  a  meaningful  uncertainty 
can  be  computed.  The  high  sample  rate  of  the  reference  computation 
requires  high-speed  computation  to  obtain  a  close  to  real-time  simulation 
and  long  word  lengths  to  limit  the  accumulation  of  round-off  errors  to 
acceptable  values.  These  two  aspects  are  the  major  sources  of  difficulty 
in  programming  the  simulation. 

The  block  entitled  accelerometer  and  gyro  simulation  performs  the  sampling 
and  quantizing  done  by  the  inertial  sensors.  Observe  that  the  outputs  are 
not  acceleration  and  rotation  rate,  but  incremental  changes  in  velocity  and 
angle,  respectively. 

A  set  of  macros  have  been  written,  so  that  tests  can  be  set  up  and  para¬ 
meters  changed  with  a  minimum  of  effort.  Where  several  tests  are  to  be 
run  on  the  same  input  data  (e.g.,  changes  in  word  length),  provision  has 
been  made  to  store  pertinent  results  of  the  reference  computation.  Hence, 
the  time-consuming  reference  computation  need  be  performed  only  once. 
Also,  an  output  routine  has  been  written  which  produces  pages  of  columnar 
data,  each  column  with  a  descriptive  title.  Numbers  can  be  printed  in 
decimal  or  octal  radices,  and  in  fixed-point  or  floating-point  forms. 
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T-i  •  f  /TTPI  Mver  M  Kessler,  William  D.  Mathews 
Technical  Information  Pro.iect  (TIP)  -  Myer  JYi.  tvessie  , 

The  M  I  T.  Libraries'  Technical  Information  Project  is  attempting  to 
ntillse  the  time-sharing  facilities  of  Project  MAC  for  literature  search 
and  retrieval.  The  system  is  a  prototype,  based  on  twenty-five  phys.cs 
journals  chosen  Iron,  a  reference  matrlst  started  by  the  Physical  Revrew. 

This  literature  extends  back  five  years,  the  average. 


For  each  journal  article,  we  record:  location  of  the  article  (journal, 

volume,  page),  title,  authors,  institutional  affiliation  of  the  authors ,  cita 

tions  (journal,  volume,  page),  location  of  the  article  in  Physics  Abstract: 
and  subject-index  information  if  available  from  a  published  source. 


This  information  is  edited,  compressed,  and  put  on  the  MAC  system  as 
data  files.  The  user  may  perform  search  and  retrieval  operations  by 
means  of  a  special  language  that  consists  of  several  dozen  words  to  control 
search  programs.  The  simplest  form  of  engagement  may  be  initiated  by 
specifying  three  commands: 

defines  the  range  of  literature  to  be  searched, 
states  the  items  to  be  found, 

defines  the  nature  and  content  of  the  computer 
output. 


SEARCH 

FIND 

OUTPUT 


several  forms: 


The  SEARCH  command  may  take 
SEARCH  ALL 

SEARCH  ALL  NEW 

SEARCH  PHYREV  ALL 


will  search  the  entire  literature 
in  store, 

will  search  the  last  volume  of 
each  journal, 

will  search  everything  in  store  of 
a  given  journal  (in  this  case  the 
Physical  Review), 

only  volumes  120  through 

135, 

only  volume  135. 

item  described 
sibilities. 


SEARCH  PHYREV  V.  120  to  V.  135 

SEARCH  PHYREV  V.  135 

The  program  will  search  as  specified  above  and  detect  any 
under  the  FIND  command,  which  itself  has  a  variety  of  pos 


96 


LIBRARY  RESEARCH 


FIND  AUTHOR  SMITH  -  will  find  all  papers  that  include 

Smith  among  its  authors. 

FIND  TITLE  CRYOGENICS  -  will  find  all  papers  that  contain 

the  word  "cryogenics"  in  the 
title. 

One  may  similarly  define  a  location, 

FIND  LOCATION  IOWA  STATE  - 
or  a  citation, 

FIND  CITATION  1  131  1165  -  to  find  all  papers  that  cite 

the  article  in  Physical 
Review  (code  =  1),  volume 
131,  page  1165.  The  find 
instruction  includes  the 
logical  manipulations  of  "and/ 
or/ but  not." 

Several  other  sophisticated  FIND  commands  are  available,  the  most  in¬ 
teresting  of  which  is: 

FIND  SHARE  B  journal,  volume,  page. 

In  this  command  we  are  looking  for  papers  in  the  search  range  that  share 
a  ciitation  item  with  a  given  paper.  This  mode  of  literature  search,  called 
bibliographic  coupling,  is  a  powerful  tool  for  information  retrieval.  The 
final  command,  OUTPUT,  selects  one  or  more  of  the  options  available  as 
output. 

/ 

The  literature  and  program  facilities  of  TIP  are  available  to  all  MAC  system 
users,  and  a  monitor  program  has  been  written  to  collect  statistics  on  the 
type  and  extent  of  MAC-TIP  utilization.  The  most  extensive  use  of  MAC- 
TIP  to  date  has  been  in  connection  with  a  book  on  Plasma  Physics,  written 
by  Prof.  Sanborn  Brown,  as  noted  later  in  this  section. 


I 

I 


/ 
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Process  Control:  Serials  and  Journals  -  Patricia  M.  Sheehan,  and 

Marianne  C.  Henneman 

As  libraries  grow  in  size,  the  problems  associated  with  maintaining  con¬ 
trols  grow  in  complexity.  Difficulties  are  most  acute  in  university 
libraries  with  scientific  commitments,  where  periodicals  represent  a 
large  proportion  of  the  holdings. 

During  1965,  TIP  has  made  a  systems  study  of  the  traditional  manual  and 
punch-card  methods  currently  used  to  control  serials  and  journals  in  the 
M.  I.  T.  Libraries.  A  feasibility  report  has  proposed  using  serials  and 
journals  as  a  microcosm  for  a  computer  approach  capable  of: 

1)  e.xtending  the  features  of  the  MACTIP  information 
retrieval  system  to  M.  I.  T.  and  other  regional  users, 

2)  offering  alternates  to  supplement  the  present  card 
catalogues, 

3)  providing  a  two-way  communication  with  the  computer, 
interjecting  reports  and  triggering  actions,  and 

4)  performing  the  usual  data-processing  functions  for 
administration,  accounting,  and  statistics. 

The  proposed  system  will  be  based  on  a  permanent  disk  file  and  a  series 
of  programs  combining  remote  console  and  direct  input  means  for  updating 
and  retrieving  data.  Initial  conversion  of  punch -card  data,  and  a 
limited  retrieval  program  have  been  completed.  Additional  programs  are 
now  being  written  to  purify  and  expand  files  to  meet  defined  objectives. 

As  one  of  the  largest  files  in  the  time-shared  system,  the  library  files 
need  careful  program  logic  to  minimize  the  effects  of  continually  expanding 
and  revising  records  and  items  of  variable  length.  A  linking  technique 
permits  simultaneous  use  of  sorted  and  unsorted  data.  Additionally,  a 
multi-level  sorter  tries  to  reconcile  file  data  having  key  breaks,  after  a 
comparatively  small  number  of  computer  characters,  with  data  having  a 
break  delayed  until  the  end  of  a  lengthy  field.  Finally,  sorting  and  re¬ 
trieval  programs  expand  a  variable  number  of  variable -length  fields  within 

.1 

a  record  and  align  each  field  with  the  corresponding  data  in  other  records. 
The  pragmatic  solution  at  this  point  is  to  have  the  programs  determine 
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rZ^lZirtoTlZr  """  necessary, 

placed  bv  current  sources  of  information  to  be  re- 

Pla  ed  by  computer -produced  catalogues  and  reports  will  be  governed,  in 
.  y  aesthetics.  For  this  reason,  we  hope  to  ultimately  have  an  off 

that  wilTl?  or  Similar  ec.uipment, 

will  be  more  versatile  than  standard  computer  printers. 

Search  Procedures  -  Evan  L.  Ivie 

‘=■‘■>■1“'  niP)  ia  basad  o„  a  „,ea.a,.  or  „.o 
^  e  ness  etween  documents  that  is  a  correlation  coefficient  based 
n  information  theory.  This  measure  was  described  in  the  Project  MAC 
Progress  Report  for  1 964  (AD-465-0881  a  ■  • 

.o  .h.  aeoa.ao.a-  or  '  I  ’•  ''  aasigood 

way  .Ko  libra,,,  la  oaod.  »»  ' 

A  paogaao,  aal  l.aa  boa,,  doviaad  ivhich  uaaa  lb, a  „,a.a„,a  to  oo„v„.  a 

pa;:;,\:::rro?a:ro“::,.r‘; 

yen  apooiria,  a„d  ,000.::?;:  ri:™ "raro: 

bo  ™  "  "’>‘'S‘»a“y-"lat.d  documants,  or  ho  ,oay 

sana.al,  and  obtab,  a  corraspoadlogly  l.,g„  o„.„„  . 

-a 

a\?a::::a:rbr::r  --  - 

English.  Some  examples  of  this  language  are; 

'■'yblch  papora  ara  rola.od  .0  papara  by  J.  R.  jonas  but  „o. 
papers  using  acoustic." 

"Pribt  ror  dociaion  ,|.a  tiilaa  ol  ariidoa  rolaiad  to  articloa 
Citing  Phys.  Rev.,  Vol.  136,  page  22." 

b  l"a'.raT‘'T‘  S'-y-D  by 

■lava  :  r: “*•  --  - 

raaldoloordaro d  ata  hay, 

l.ata  hayiog  coanootiyoa  to  othar  paramotara. 
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Measures  of  Relatediicss  -  John  P.  Casey 

The  volume  of  journal  literature  is  increasing  -rapidly.  An  understanding 
of  the  structure  of  this  literature  may  provide  some  interesting  informa¬ 
tion  about  an  e.xpanding  communication  network,  and  may  help  in  the 
design  of  systems  to  manipulate  and  order  these  articles.  It  is  not  clear 
what  type  of  structure  will  fit  journal  literature,  so  tentative  beginnings 
have  been  made  in  three  different  directions. 

Eirst,  a  collection  of  articles  may  be  considered  as  a  graph,  with  the 
articles  as  points  on  the  graph.  An  arc  connects  two  articles  if  they  both 
cite  the  same  article.  Questions  can  now  be  asked  how  tightly  a  subgraph 
is  connected  together  by  its  arcs,  or  how  closely  connected  an  article  is 
to  a  subgraph.  These  questions  are  closely  related  to  the  difficulties  of 
designing  a  mechanical  system  to  find  a  group  of  articles  all  related  to 
a  given  reference. 

Secondly,  this  same  network  may  be  considered  as  a  group  of  nodes  con¬ 
nected  by  bars,  with  the  subject  discussed  in  a  given  article  diffusing  like 
heat  along  the  bars.  Coefficients  of  permeability  may  be  assigned  to  the 
bars,  and  rather  complex  formulas  derived  for  the  amount  of  "heat"  e.x- 
changed  between  an  article  and  others  scattered  through  various  journals, 
to  determine  which  are  most  closely  related  to  it. 

Thirdly,  an  analysis  using  the  methods  of  Rashevsky's  mathematical 
biology  has  been  started.  From  a  number  of  simplifying  assumptions 
about  the  distribution  of  useful  ideas  in  articles  and  related  factors,  a 
simple  model  may  be  constructed  which  predicts  the  likelihood  of  an 
article  being  cited  and  the  time  span  until  an  article  is  no  longer  men¬ 
tioned  in  the  literature.  This  may  then  be  checked  against  statistical 
data  available  through  the  TIP  system.  It  is  hoped  that  a  model  which 
makes  reasonably  accurate  predictions  will  clarify  the  general  charac¬ 
teristics  desirable  in  a  literature-searching  system. 

We  are  investigating,  for  a  short  distance  at  least,  each  of  these  ap¬ 
proaches.  Thus  far,  there  have  been  no  results  that  would  show  the 
feasibility,  or  lack  of  such,  for  any  of  these. 
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Statistics  of  Words  in  Titles  -  Elizabeth  A.  Dole 


Present  studies  attempting  to  describe  several  characteristics  of  TIP's 
vocabulary  are  based  on  word-count  lists  compiled  from  six  years  of 

Pl^sical  Review,  Volumes  113-132  (1958-1963).  and  Volumes  133-136 
(1964)  Series  A. 


Several  criteria  are  used  to  assessed  total  vocabulary,  gross  word  count, 
and  average  words  per  title.  The  basis  of  each  criterion  is  an  exact 
definition  by  formula  or  by  known  arbitrary  rule.  The  first  criteria 
developed  are  those  describing  what  is  a  "word,"  and  the  groups  defined 
are  successively  rejected  from  the  original  unrefined  vocabulary.  The 
following  steps  are  planned: 


1.  A  "word"  is  any  string  of  characters  between  two  spaces 
Sample  Results; 


number  of  individual  words 
individual  words  times  number  of  uses 
number  of  articles  in  sample 
average  number  of  words  per  title 


4, 826. 
56,  413. 

7, 583. 

7.  43 


2.  Refinement  A.  A  word  is  any  string  of  characters  between 
two  spaces  that  does  not  start  with: 

A.  a  number 

B.  +or-or= 

C.  any  punctuation  mark 


3.  Refinement  B_.  (in  progress)  A  word  is  any  string  of 
characters  between  two  spaces  not  ruled  out  by  step  number  two.  above. 
Rat  is  not  one  of  an  arbitrarily-specified  list  of  words  selected  for  their 
mherent  lack  of  information  context,  such  as  connectives. 

In  this  step,  the  list  of  words  is  selected  and  tabulated.  The 
impact  of  removing  these  items  from  the  word  list  GROSS  (A)  will  be 
calculated  to  give  a  new  GROSS  (B). 

4.  Refinement  C.  A  word  is  any  string  of  characters  be¬ 
tween  two  spaces  not  ruled  out  in  steps  two  and  three,  above,  that  bears 
basic  similarity  to  any  other  word. 
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In  this  case,  words  with  notably  similar  bases  or  roots  are 
grouped,  and  these  counts'  are  compounded  in  the  total  figures.  The 
result  should  be  a  reduction  in  the  number  of  individual  words,  but  not  in 
the  gross  number  used  for  this  step. 

5.  Refinement  D.  Remove  any  remaining  combination  of 
letters  and  numbers  that  is  not: 

a)  an  obvious  symbol  for  a  chemical  compound, 

b)  an  obvious  isotope  symbol. 

6.  Refinement  E_.  Remove  all  chemical  compound  symbols. 

7.  Refinement  F.  Remove  all  isotope  symbols. 

Refinement  G.  Remove  all  unclumped  words  occurring 
five  times  or  less.  Evaluate  the  remaining  vocabulary.  List  all  words 
so  removed. 

9.  Comparison  I.  Compare  a  portion  (or  more)  of  the  results 
above  with  those  gained  from  an  exactly  similar  list  and  count  from  TIP's 
present  full  file  of  25  journals* 

10.  Comparison  II.  Compare  search  products  from  equal 
periods  of  TIP's  library  with  the  library  of  physics  references  available 
through  equivalent  editions  of  Physics  Abstracts,  and  through  the  Physics 
Abstracts  Subject  Index  Headings  List  effective  at  the  same  time. 

In  summary,  these  vocabulary  and  retrieval  studies  should  show  the 
breakdown,  from  an  information  standpoint,  of  the  vocabulary  of  title  words 
of  physics  articles  as  represented  in  the  data  base.  The  studies  are  dis¬ 
tinguished  by  clearly  defined  stages.  These  stages  are  defined  explicitly 
enough  that,  once  accepted,  they  could  easily  be  expressed  in  logical  pro¬ 
gramming  for  computer  processing.  The  studies  so  outlined  will  afford 
needed  information  to  designers  of  information  files,  particularly  of  in¬ 
verted  title  word  files,  in  the  field  of  scientific  literature. 
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Educational  Use  of  TIP  -  Myer  M.  Kessler,  Sanborn  C.  Brown,  and 

Irma  Y.  Johnson 

A  freshman  seminar  based  on  TIP  is  being  offered  in  the  fall  of  1965, 

This  seminar  is  conducted  with  the  cooperation  of  Prof.  Sanborn  C.  Brown, 
of  the  Physics  Department,  and  Mrs.  Irma  Johnson,  Reference  Librarian 
at  the  M.  I.  T.  Libraries. 

Students  will  investigate  the  literature  of  plasma  physics  using  TIP  and 
compare  it  with  the  results  of  a  search  based  on  traditional  library 
methods.  The  course  is  experimental  and  is  being  observed  from  two 
points  of  view: 

1.  Can  beginning  students  more  profitably  study  a  field  of 
physics  from  current  literature  references  on  a  computer,  instead  of 
traditional  text  books? 

2.  Can  beginning  students  compile  useful  bibliographies  and 
compilations  of  physics  literature  from  a  computer-based  system? 
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Use  of  TIP  to  Update  a  Data  Compilation  -  Sanborn  C.  Brown 

I  have  been  investigating  the  possibility  of  writing  a  report  which  will 
demonstrate  using  a  computer  to  update  a  compilation  of  scientific  data. 
Six  years  ago  I  published  a  reference  book,*  and  the  present  revision  is 
based  on  bibliographic  material  contained  in  the  physics  literature  pro¬ 
grammed  in  the  M.  I.  T.  Library's  Technical  Information  Project  (TIP). 
This  program  is  available  to  all  Project  MAC  users. 

Using  this  system,  the  attempt  is  to  create  material  in  an  open-ended 
form,  so  that  anyone  with  access  to  the  computer  program  can  search  the 
literature  for  material  which  will  appear  after  the  text  of  my  book  is 
printed.  As  an  aid  for  doing  this,  references  to  individual  displays  of 
data  are  given  in  com.puter-language  form,  rather  than  the  usual  biblio¬ 
graphic  form. 

To  illustrate  the  method  of  keeping  up-to-date,  let  me  give  an  example. 
Suppose  that  the  reader  is  interested  in  the  elastic  collision  cross  section 
between  electronics  and  drypton  atoms.  He  finds  a  collection  of  data  from 
L.  S.  Frost  and  A.  V.  Phelps  (1964)  PHYREV  V00136  1538.  To  find  out 
what  further  work  has  been  published,  assume  that  new  work  would  cite 
this  paper  in  its  bibliography.  The  computer  program  is,  therefore, 
searched  to  see  if  this  has  occurred. 

First,  the  reader  must  log  in  to  the  computer  in  a  manner  appropriate  to 
his  connection  with  the  computer.  Once  logged  in,  he  calls  for  the  entire 
library  file  with  the  command  TIP  and  then  asks  the  computer  to  search 
for  the  citation  and  deliver  the  information  with  the  command  to  search 
ALL  and  find  the  appropriate  citation,  by  typing  as  output  the  author,  title, 
and  identification.  The  command  actually  typed  into  the  console  is; 

S  ALL  from  1964 
F  CITE  PHYREV  V00136  1538 
O  PATI 


'"Basic  Data  of  Plasma  Physics,  John  Wiley  &  Sons,  Inc.,  New  York,  1959. 
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On  the  assumption  that  some  new  reference  was  found,  the  reader  then 
may  want  to  browse  further  in  the  associated  literature.  This  is  most 
efficiently  done  by  having  the  computer  search  for  all  papers  which  share 
^ast  one  bibliographic  reference  with  the  new  paper  just  found.  This 
can,  of  course,  be  done  against  the  entire  library  file,  but  since  the 
material  of  "Basic  Data"  was  compiled  using  this  technique  of  literature 
search,  computer  time  can  be  saved  by  excluding  from  the  search  all  those 
papers  which  already  appear  in  the  text.  This  can  be  accomplished  by  a 
subroutine  called  "BROWN"  created  in  the  process  of  putting  this  present 
VO  ume  together.  The  command  on  the  console  to  get  a  shared  bibliography 
based  on  the  program  for  this  volume  would  be: 

S  ALL  BUT  NOT  BROWN 

F  SHARE  B  PHYREV  VOOI6O  1234 

O  PATI 


Dr.  MyerM.  Kessler.  Associate  Director  of  Libraries  at  M.  I.  T  and 

Director  of  Project  TIP.  studied  statistically  where  most  of  the  physics 
iterature  was  published  and  developed  a  list  of  25  journals  which  contained 
percent  of  physics  papers  (Physics  Today.  18.  28.  1965).  In  develop- 
mg  the  "BROWN"  program  for  this  report,  certain  journals  were  further 
excluded  as  not  being  appropriate  to  experimental  gaseous  electronics  and 
p  asma  physics.  The  list  of  journals  in  the  "BROWN"  program  and  their 
program  code  designation  is  given  in  the  following  list: 


Applied  Physics  Letters  a  ppr  p.-- 

Canadian  Journal  of  Physics  PHYPAM 

Helvetica  Physica  Acta  'SyS^? 

Indian  Journal  of  Physics  INDtph 

Japanese  Journal  of  Applied  Physics  '  ‘  PHAPJA 

Journal  of  Applied  Physics  PHyapp 

Journal  of  Chemical  Physics  JCHPPP 

PhysS  Society  of  Japan  PHYSOJ 

Physics  Letters  PHYLET 

The  Physics  of  Fluids  ■ 

Physical  Review  A  PHYRPV 

Physical  Review  Letters  PHYRPT 

Proceedings  of  the  Physical  Society  (London)  PHYPRO 
Soviet  Physics  -  JETP  '  cp/pl^ 

Soviet  Physics  -  Technical  Physics  SPTPHY 
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In  addition  the  total  TIP  program  will  also  search: 

Annals  of  Physics 

II  Nuovo  Cimento 

Nuclear  Physics 

Physical  Review,  Series  B 

Progress  of  Theoretical  Physics  (Kyoto) 

Soviet  Physics  -  Solid  State 


ANNPHY 

NUOCIM 

NUCPHY 

PHYREB 

PROPHJ 

SPSOLS 


Another  source  of  data  which  are  not  always  to  be  found  in  the  periodical 
literature  are  papers  published  in  the  proceedings  of  conferences.  A 
number  of  th^se  were  of  importance  to  the  material  of  this  report  and 
these  have  been  especially  incorporated  into  the  TIP  program.  They  are 
Proceedings  of  the  Fifth  International  Conference  on  Ionization  Phenomena 
in  Gases  held  in  Munich  in  1961  (MUNICH),  Comptes  Rendus  de  la  VI— 
Conference  Internationale  sur  les  Phenomenes  d'lonisation  dans  les  Gaz 
held  in  Paris  in  1963  (OPARIS),  and  Proceedings  of  the  Seventh  Inter¬ 
national  Conference  on  Ionization  Phenomena  in  Gases  held  in  Belgrade  in 
1965  (BELGRA). 


If  the  reader  does  not  wish  to  exclude  the  already  search-shared  biblio¬ 
graphies  in  the  "BROWN"  program,  both  lists  of  journals  will  be  searched 
when  "BUT  NOT  BROWN"  is  omitted  from  the  console  command. 

This  study  is  aimed  at  serving  two  purposes:  one  is  to  bring  the  "Basic 
Data"  book  up-to-date;  the  other  is  to  provide  a  computer  method  of  find- 
ing  new  data  which  appear  after  the  report  is  issued.  The  second  purpose 
is  accomplished  through  the  citation  process,  and,  therefore,  will  be 
effective  only  if  there  is  some  article  to  CITE.  Thus,  if  no  published  work 
has  been  found  which  supercedes  that  published  in  the  original  "Basic  Data" 
book,  the  previous  volume's  citation  will  be  repeated  with  its  proper 
reference  coding,  so  that  subsequent  material  may  be  found  by  the  com¬ 
puter-search  method. 

This  is  possible,  because,  although  the  TIP  program  does  not  contain 
literature  published  prior  to  1959,  the  bibliographic  coding  is  complete.  If 
the  reader  so  desires,  he  may  call  for  the  computer  to  present  him  with 
any  complete  bibliography  he  desires  by  calling  for  "B"  as  an  output 
command. 
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In  searching  the  literature,  it  is  sometimes  desirable  to  determine  what 
the  common  bibliographic  references  are  between  two  articles  found  by 
the  'SHARE  B"  technique.  This  is  provided  for  in  the  TIP  program  by 
asking  the  computer  to  print  the  linkage  between  any  two  papers. 

The  program  has  been  written  so  that  it  is  not  necessary  for  a  cited 
reference  to  be  for  articles  actually  available  in  the  program  by  title, 
author,  etc.,  but  only  that  it  be-in  the  references.  Thus,  the  classic 
paper  on  elastic  collision  was  published  by  Erode  in  the  Reviews  of 
Modern  Physics  in  1932,  long  before  the  time- interval  covered  by  TIP, 
but  use  of  the  LINKAGE  request  will  uncover  the  proper  volume  and  page 

location  of  Brode's  paper  if  it  has  been  used  as  a  reference  in  a  recent 
article. 
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Baseband  Design  of  a  Unified  Carrier  -  Homer  C.  Peterson 

Work  on  a  design-optimization  prograrn  was  started  using  FORTRAN  and 
batch-processing  compiling.  The  very  nature  of  the  equation  and  the 
technique  of  handling  were  not  clearly  defined  and  the  infinite  limits  on 
the  summation  made  it  hard  to  determine  what  intervals  to  investigate. 

The  first  runs  were  so  long  that  they  were  dumped  from  the  Lincoln 
computer  without  progressing  far  enough  to  even  show  a  trend  to  the 
results.  With  usually  20  hours  elapsing  between  submitting  a  run  and 
seeing  the  results,  it  became  apparent  that  the  problem  would  either 
have  to  be  over-simplified  or  dropped. 

At  this  time  CTSS  became  available  to  use.  The  first  big  step  was  to 
recode  the  program  in  the  MAD  language.  The  program  was  then  modi¬ 
fied  to  be  handled  a  section  at  a  time;  where  hand  comparisons  could  be 
made,  a  rough  plot  obtained  from  the  teletype,  and  programming  errors 
readily  corrected.  Within  a  week  the  first  outputs  were,  available,  while 
the  original  had  dragged  out  for  nearly  a  year.  Considerable  rewrite  and 
modification  was  necessary,  as  it  became  clear  that  certain  sections  of 
the  program  did  not  agree  either  with  the  model  or  with  the  actual  world. 

A  new  Bessel  function  subroutine  was  written,  to  obtain  the  desired 
accuracy,  only  to  find  that  it  worked  about  half  the  time.  The  CTSS 
loader  didn't  always  properly  load  a  double-precision  hardware  program, 
so  the  program  had  to  be  rewritten  to  be  self-shifting  if  it  was  started  in 
an  odd  register. 

For  nearly  a  calendar  year  (but  probably  less  than  two  man-months  of 
work)  the  program  was  in  constant  change  with  as  many  as  ten  recompila¬ 
tions  occurring  in  one  day.  During  this  time,  the  program  was  lost  only 
three  times  due  to  computer  failure.  The  difference  between  using  a  CTSS 
system  and  batch-processing  was  a  solution  with  about  two  man-months 
of  work  as  opposed  to  abandoning  the  project. 
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On-Line  Data  Storage  and  Retrieval  -  John  F.  Nolan 

During  the  past  year,  the  on-line  data  storage  and  retrieval  systeni  (as 
described  in  last  year's  progress  report)  was  coded,  checked,  and 
operated  on  a  variety  of  sample  files.  Although  limited  in  data  capacity 
by  the  disk  space  allocated  to  the  research  problem,  the  experimental 
system  demonstrated  successful  implementation  of  advanced  techniques 
for  on-line  access  and  manipulation  of  formatted  data  files.  The  most 
significant  features  of  the  design  are; 

1.  porage  of  coded  information  with  the  data  files  describing 
their  organization  and  content. 

2.  Representation  of  this  descriptive  information  by  the 
identical  conventions  used  to  store  the  datafiles  themselves:  thus  all 
system  commands  for  searching,  ’  cross-associating,  and  manipulating 
data  files  can  be  used  to  operate  upon  the  file  descriptions  as  well. 

3.  The  internal  use  of  list-structures  to  allow  unconstrained 
cross-association  of  filed  data  (including,  as  above,  the  file  descriptions) 
according  to  the  user's  needs. 

4.  Simple  one-word  commands  and  dynamic  instructions  which 
allow  immediate  use  without  any  indoctrination  or  programming  knowledge. 
(See  Nolan,  Appendix  C) 

Compilation  for  a  Digital  Differential  Analyzer  -  Harold  K.  Knudsen 

Lincoln  Laboratory  is  constructing  a  Digital  Differential  Analyzer  (DDA) 
as  part  of  a  hybrid  computer  facility.  A  compiler  is  being  written  on  the 
MAC  system  to  make  it  possible  for  the  nonspecialist  to  generate  solutions 
to  ordinary  differential  equations  on  the  DDA,  The  inputs  to  the  compiler 
will  be: 

1.  The  diffeJ'ential  equation  written  in  a  FORTRAN-like 

language, 

2.  The  initial  conditions  of  the  variables,  and 

3.  Estimated  maximum  absolute  values  of  the  variables. 
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Tlie  output  of  the  compiler  will  be; 

1.  An  interconnection  matrix  which  specifies  the  internal 
connections  of  the  integrators  and  servos  in  the  DBA,  and 

2.  The  scaled  initial  conditions  of  the  problem  variables. 

The  compiler  outputs  will  be  loaded  directly  into  the  DBA  from  the  MAC 
computer  through  a  channel  consisting  of  a  Model  35  teletype  and  the 
LINC  computer. 

On-Line  Experimentation  -  Joseph  W.  Bex 

Buring  the  past  year,  CTSS  was  used  for  the  following  activities: 

1.  Searching  for  sequences  with  desirable  autocorrelation 
properties,  for  use  in  coding  patterns  for  the  Arecibo  radar  antenna  in 
Puerto  Rico, 

2.  Testing  the  stability  of  a  new  circuit  design  to  be  used  to 
control  the  carbon  arc  for  Lincoln  Laboratory's  new  environmental 
facility, 

3.  Simulating  designs  for  magnetic  orientation  systems  used 
on  Lincoln's  experimental  communications  satellite,  LES. 

Exploratory  on-line  computation,  done  with  the  responsible  scientist  or 
engineer  present,  was  performed  for: 

1.  Besigning  antenna  arrays  for  use  with  the  LES  satellite, 
involving  the  adjustment  and  testing  of  various  parameter  values, 

2.  Betermining  the  stability  or  instability  of  the  theoretical 
model  for  a  new  solid-state  device  which  would  have  unusual  high-gain 
prope  rties . 
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In  addition,  th,  Projoot  MAC  computer  sy.tem  began  to  be  used  by  staff 
o  e  Bell  Telephone  Laboratories  and  of  the  General  Electric  Company 
participating  in  the  development  of  the  MULTICS  system. 
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Model  Testing  -  Merrill  M.  Flood  (University  of  Michigan) 

The  purpose  of  the  model  testing  project  is  to  develop  programs  and  on¬ 
line  techniques  adequate  for  parameter  estimation  and  goodness -of-fit 
testing  on  several  stochastic  process  models.  These  models  derive  from 
the  data  of  various  learning  and  decision  experiments  with  human  subjects 
and  animals.  We  have  adapted  programs  used  previously  on  IBM  7094 
computers  at  the  University  of  Michigan,  at  the  University  of  California, 
and  at  the  System  Development  Corporation,  for  CTSS  use.  The  pro¬ 
grams  are  nov'  suitable  for  one  specific  stochastic  model  and  have  been 
used  successfully  on  a  few  small  test  runs. 

One  set  of  FORTRAN  Subroutines,  called  LOOK,  are  part  of  the  general 
program,  STOl.  They  are  an  adaptation  of  a  direct  search  code  developed 
at  the  Wcstinghouse  Electric  Company  for  the  iterative  solution  of  maxi¬ 
mization  problems  with  constraints.  The  LOOK  sub-routines  will  be  use¬ 
ful  as  part  of  other  computer  codes  and  include  optimization  calculations. 
Eventually,  we  hope  to  make  LOOK  and  other  similar  optimization  codes 
available  to  all  interested  CTSS  users.  One  FAP  Subroutine,  called 
PRODUC,  within  STOl  was  written  to  calculate  rapidly  and  accurately  a 
lengthy  product  of  probabilities.  Some  of  our  current  programming  is 
devoted  to  the  development  of  improved  codes  for  this  calculation  and 
CTSS  is  useful  in  testing  such  alternative  codes  for  speed  and  accuracy. 

Wa  are  using  these  optimization  codes  to  investigate,  with  experimental 
data,  a  few  stochastic  process  models  of  rat  behavior. 
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Discovery  and  Learning  Techniques  and  Pro£rrammed  Instructions  - 
Leonard  Uhr  (University  of  Michigan) 


Wc  are  working  on  a  program  that  learns  how  to  handle  one  or  two-dimen¬ 
sional  inputs,  such  as  language  strings  and  sensory  patterns.  With  the 
exception  of  those  parts  that  build  up  class  structures,  sections  of  the 
program  that  use  previously-learned  information  to  generate  the  appropriate 
response  to  an  input  have  been  largely  completed  and  debugged.  The  sec¬ 
tions  that  learn  have  been  coded,  and  are  now  being  debugged.  We  are  also 
continuing  tests  of  the  earlier  program  that  handles  one-dimensional  strings 
and  we  are  looking  into  methods  for  increasing  its  abilities  to  handle 
different  types  of  problems. 


In  addition,  we  have  coded  a  program  that  turns  written  materials,  such  as 
textbooks  and  speck  lly-prepared  instructional  materials,  into  rough  drafts 
of  teaching-machine  programs,  conveniently  edits  such  programs  on-line, 
and  presents  the  finished  instructional  program  to  students.  In  this  pro-  ' 
gram,  content  material  is  interspersed  with  short  fill-in  questions  that 
lead  to  forward  or  backward  branches  as  a  function  of  the  student's  success 
or  failure  in  giving  a  correct  or  synonymo'us  answer.  We  can  now  generate 
teaching-machine  programs  in  three  ways:  by  compiling  from  materials 
prepared  according  to  a  simple  standard  format;  by  generating,  during - 
interactions  with  the  student,  from^a  standard  description  of  a  problem 
domain;  and  by  converting  ordinary  books  into  programs.  We  should  now 

like  to  combine  these  methods  into  a  single  prdgram  that  would  offer  all 
these  options. 
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Generalized  Desk  Calculator  -  Anthony  G.  Oettinger  and 

Adrian  Ruyle  (Harvard  University) 

Project  TACT  (Technological  Aids  to  Creative  Thought)  has  been  develop¬ 
ing  the  on-line  system  of  Glen  Culler,  which  functions  roughly  as  a  gen¬ 
eralized  desk  calculator.  With  this  system,  not  only  can  arithmetic  opera¬ 
tions  be  carried  out  on  single  numbers,  as  they  would  on  a  typical  desk 
calculator,  but  operations  of  more  generality  —  such  as  exponentiation 
and  forward  differencing  —  can  be  carried  out  on  functions  considered  as 
such. 

Execution  is  as  quick  and  easy  to  initiate  as  using  a  desk  calculator. 
Operations  on  complex  functions  are  also  possible,  as  are  plots  of  functions 
on  a  cathode-ray-tube  screen,  and  display  of  text.  There  is  also  a  facility 
for  on-line  storing  of  "programs"  of  operator  sequences,  which  can  be 
executed  en  bloc  as  single  calculator  operations. 

The  calculator  is  now  developed  to  the  stage  where  it  is  a  useful  tool,  for 
about  a  half-dozen  users  who  will  be  variously; 

1.  Preparing  demonstrations,  such  as  prototypes  of  course 
lectures  in  mathematics,  and  giving  calculator  exhibitions, 

2.  Conducting  research,  such  as  kidney  physiology  and 
frog  demography, 

3.  Exploring  the  imaginative  possibilities  of  the  calculator 
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Aspects  of  Integer  Linear  Programming  -  R6mulo  H.  Gonzalez-Zubieta 

Under  U.  S,  Army  Research  Office  (Durham)  Contract  No.  DA-31 -124- 
ARO-D-209,  a  primal  feasible  (all -integer)  integer  linear-programming 
algorithm  has  been  developed  and  programmed,  together  with  a  related 
proceduie  for  obtaining  a  first  feasible  solution.  Once  a  feasible  solution 
is  found,  the  algorithm  maintains  feasibility  at  each  stage,  in  contrast  to 
other  algorithms  that  have  been  programmed  and  are  currently  available. 
These  other  algorithms  do  not  achieve  feasibility  until  the  optimal  solution 
is  reached.  The  primal  feasible  algorithm  is  based  on  a  particular  way  of 
applying  the  cutting  planes  previously  developed  by  R.  E.  Gomory,  and  on 
a  specific  interpretation  of  their  role. 

The  finiteness  of  convergence  has  been  established  for  two-dimensional 
problems  but  not  for  the  general  case;  however,  there  appears  to  be  at 
least  computational  convergence  in  a  considerable  fraction  of  the  cases. 

In  addition,  a  Generalized  Euclidean  Algorithm  for  finding  the  greatest 
common  divisor  for  more  than  two  numbers  has  been  defined.  The  solution 
of  systems  of  linear  diophantine  equations  is  presented  in  terms  of  integer 
linear  programming. 
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Introduction 

This  interdepartmental  laboratory  provides  facilities  for  academic  research 
in  three  categories,  designated  as  general  physics,  plasma  dynamics  and 
communication  sciences. 

Major  support  for  the  research  is  provided  by  the  Joint  Services  Electronics 
Program  of  the  Army,  Navy  and  Air  Force  as  well  as  the  Atomic  E^nergy 
Commission,  the  National  Science  Foundation,  the  National  Institutes  of 
Health  and  the  National  Aeronautics  and  Space  Administration. 

As  indicated  by  the  following  section,  a  substantial  number  of  RLE  research 
projects  have  received  invaluable  assistance  from  Project  MAC  facilities. 
These  faculty  and  student  research  activities  span  a  wide  range  of  scien¬ 
tific  and  engineering  subjects. 


Dynamics  of  Active  Plasma  Systems  -  Abraham  Bers 

During  the  past  year  we  have  made  extensive  use  of  the  MAC  computer 
and  the  ESL  display  system  to  solve  a  large  number  of  difficult  and  impor¬ 
tant  problems  in  plasma  stability.  The  complexity  of  these  problems  is 
such  that  without  direct  man-computer  interaction  many  of  these  problems 
would  not  have  been  attacked.  The  basic  problem  of  plasma  stability  to 
small-signal  perturbations  was  described  in  our  progress  report  last  year. 
The  most  significant  programs  developed  during  the  last  year  for  our  prob¬ 
lems  are  summarized  in  three  reports  that  follow  (M.  A.  Lieberman;  C. 
Speck  and  M.  Pennell;  and  J .  D.  Mills  and  A.  Bers). 


Our  project  has  brought  a  large  number  of  our  students  in  contact  with  the 
Project  MAC  facilities.  The  thesis  research  of  these  students  has  ranged 
from  computer-display  system  programming  to  computer  stability  analyses 
of  hot,  gaseous  plasmas,  and  plasmas  in  solids.  From  our  group,  a  total 
of  about  a  dozen  graduate  students,  three  faculty  members,  and  two  post¬ 
doctoral  visitors  have  made  use  of  the  Project  MAC  facilities  during  this 

past  year. 
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JJispersion  Relation  for  Hot  Plasmas  -  Michael  A.  Lieberman 

We  have  written  a  computer  program  that  finds  the  zeros  of  a  trans¬ 
cendental  dispersion  function  D(^,  k,  .  .  .  )  in  the  complex  uj  plane.  Given 
the  complex  frequency^,  the  wave  number  k.  and  any  other  parameters 
that  a  user  may  desire,  he  must  provide  a  subroutine  within  CTSS  that 
computes  the  value  of  D.  This  subroutine  may  be  written  in  MAD, 
FORTRAN,  or  FAP. 

In  operation,  the  program  continuously  steps  up  k  by  an  increment  Alt  and, 
at  each  step,  finds  a  zero,  of  the  dis,persion  function  D; 


D{co  ,  k+  nAlv,  .  .  .  )  =  0. 


The  zero,  is  found  through  the  construction  of  a  grid  of  values  in  the 
complex  u  plane  around  an  initial  guess  to  as  follows- 

gn 


=  provided  by  the  user 


to  ,  =  to 
-gl  -O 

“gn=^n-l  +  ^<^n-l  >  2. 


The  program  evaluates  the  real  function  F  =  |  D{u,  k+  itAlt,  .  .  .  )  for  every 

point  on  the  grid.  Whenever  a  minimum  of  F  is  found  at  to  ,  the  grid  is 

n  “ 

refined  several  times  until  to is  given  to  three  significant  figures;  then,  D 
is  checked  to  verify  whether  its  real  and  imaginary  parts  change  sign  in 
the  neighborhood  and,  if  they  do,  the  zero  is  printed  and  k  is  stepped 

up.  Thus,  the  zeros  of  D  in  the  complex^  plane  are  computed  as  a  func¬ 
tion  of  the  wave  number  k. 


Without  time  sharing,  a  solution  of  this  problem  would  be  prohibitive;  with 
it,  a  user  may  specify  and  alter  at  will  the  grid  size  and  spacing  in  the  com 
plex  u  plane,  the  wave  number  k  and  its  increment  Ak,  and  all  other  param 
eters.  If  at  any  step  a  zero  of  D  is  not  found,  the  program  requests  the 
user  to  change  the  grid  size,  spacing,  and  location  in  the  complex  uplane. 

To  help  the  user,  the  program  will  print  the  values  of  D  and  F  at  the  grid 
points  if  they  are  requested. 
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Plasma  Dispersion  Relations  with  Infinite  Roots  -  Carlton  E.  Speck  and 

Martha  M.  Pennell 


In  the  study  of  a  particular  class  of  plasma  waves,  the  following  equation 
must  be  solved  for  its  roots,  VC  =  RVC  +  jCVC: 


1  -  (VPC)^ 


(VC  -  n) 


’^^Jn<P>  i 

(VC-n)^  j 


=  0 


We  had  to  determine  those  values  of  px  and  VPC  which  had  ranges  of  p 
that  resulted  in  roots  with  a  negative  imaginary  part.  Since  the  equation 
has  an  infinite  number  of  roots  and  is  not  expressible  as  a  polynomial, 
we  used  Newton's  method  to  locate  the  roots.  However,  Newton's  method 
requires  an  accurate  guess  of  the  root  location.  To  obtain  these  guesses, 
we  wrote  subsidiary  programs  for  the  approximate  portion  of  the  root. 
These  approximations  derive  from  an  approximate  equation  in  only  those 
terms  of  the  infinite  sum  which  are  large  when  |VC|  ~n.  With  such 
"guessing"  programs  an  operator  can  find  the  required  approximation  for 
the  more  exact  Newton's  method  and,  by  varying  VPC  and  px,  can  rapidly 
search  the  VC  -  p  plane  until  he  finds  the  approximate  regions  of  the  nega¬ 
tive  imaginary  part  of  VC. 
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Stability  Analysis  of  Dispersion  Relations  -  James  D.  Mills  and 

Abraham  Bers 

A  program  has  been  written,  checked  out,  and  is  now  being  used  by  mem¬ 
bers  of  the  Beam-Plasma  Group  of  the  Research  Laboratory  of  Electronics 
in  the  stability  analysis  of  plasma  dispersion  relations.  The  program 
makes  use  of  the  ESL  Display  Console  for  graphical  outputs.  The  time¬ 
sharing  system  facilitates  man -machine  communication  for  changing  pa¬ 
rameters  and  equations  in  the  analysis. 

The  motivation  for  the  prog’ am  came  through  the  formulation  of  the  Bers- 
Briggs  stability  criteria.  These  criteria  enable  one  to  analyze  the  disper¬ 
sion  relation  of  a  plasma  (or  other)  system,  to  determine  and  classify  the 
kinds  of  stable  and  unstable  wave  responses  which  are  inherent  to  the  sys¬ 
tem.  The  technique  involves  the  general  problem  of  mapping  a  function 
of  a  complex  variable.  The  dispersion  equation  of  the  physical  system  to 
be  analyzed,  D(co,  k)  =  0,  relates  the  two  complex  variables  which  are 

Gj  =  oj  +  joj.  and  k  =  k  +  jk..  One  then  maps  lines  of  constant  oj  or  tj.  into 
r  i  r  i  ^  ri 

the  complex  k-plane  or  lines  of  constant  k^  or  k.  into  the  complex  oj -plane. 

The  program  restricts  the  user  to  polynomial-type  equations  of  no  more 
than  degree  twenty.  The  coefficients  of  the  terms  of  the  equation  are  speci¬ 
fied  simply  in  a  subroutine  provided  by  the  user.  He  then  translates  this 
subroutine  using  the  MAD  command.  It  is  then  loaded  together  with  the 
display  program  which  queries  the  user  concerning  various  parameters. 
When  these  are  specified,  the  program  makes  all  of  the  calculations  neces¬ 
sary  to  the  mapping  and  displays  the  results  on  the  Display  Console.  In 
Figure  9  for  example,  the  user  specifies  the  lines  indicated  on  the  complex 
to -plane  (a)  and  he  sees  the  complex  k-plane  (b)  plotted  on  the  display. 

The  program  has  been  used  in  various  stages  of  development  since  April, 
1964.  A  joint  RLE  -  Project  MAC  technical  report  is  being  prepared  for 
publication.  (See  J.  D.  Mills  thesis  entry  in  Appendix  B.  ) 
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Figure  9.  Mapping  of  Dispersion  Relation  D(oj,  k) 
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Use  of  CTSS  in  a  Plasma  Physics  Experiment  -  David  T.  Llewellyn -Jones 

This  work  combines  automatic  data  acquisition  and  processing  with  a  labor¬ 
atory  experiment  in  plasma  physics  which  requires  rapid  feedback  between 
experimenter  and  his  apparatus.  The  last  part  of  this  work  will  establish 
an  interface  between  a  digital  voltmeter,  connected  to  the  laboratory  ex¬ 
periment,  and  a  teletype,  connected  to  the  MAC  computer.  Since  a  tele¬ 
type  has  been  installed  beside  the  experiment,  the  experimenter  has  been 
able  to  simulate  actual  runs  and  at  the  same  time,  to  act  and  manually 
type  the  data  into  the  teletype  as  it  is  displayed  on  the  voltmeter.  In  this 
way,  speeds  of  approximately  one,  three-digit  data-point  in  every  2.  5 
seconds  are  possible.  Since  the  experiment  can  run  at  speeds  up  to  about 
one  data-point-per-second,  realistic  simulation  is  only  possible  for  slow 
runs . 

In  simulation  runs,  the  .  TAPE,  feature  of  CTSS  is  very  useful.  This  fea¬ 
ture  immediately  creates  a  permanent  data  record  before  the  data  reaches 
the  program,  and  it  has  all  the  useful  correction  and  substitution  facilities 
afforded  by  the  INPUT  mode.  When  the  data  was  read  directly  into  the 
program  and  the  necessary  correction  and  substitution  facilities  had  to  be 
written  into  the  program  itself,  the  program  became  cumbersome  and  slow. 
After  installation  of  a  reliable  automatic  recording  system,  the  data  will 
be  read  directly  into  the  program  and  the  program  itself  will  automatically 
create  a  .  TAPE,  of  the  data  as  a  record. 


RESEARCH  LABORATORY  OF  ELECTRONICS 


131 


Analysis  of  Speech  -  John  M.  Heinz  and  Eleanor  C.  River 

As  part  of  a  study  on  the  constraints  which  the  vocal  organs  impose  on  a 
speech  signal  as  a  result  of  their  physical  configuration,  we  are  trying 
to  model  the  process  of  speech  production  from  its  description  by  lateral 
cineradiographs  up  to  the  resulting  acoustic  signal.  Four  stages  exist  in 
analyzing  these  lateral  cineradiographs  and  relating  them  to  the  corres¬ 
ponding  quasi-static  acoustic  spectra:  1)  obtaining  from  the  film  a  frame- 
by  frame  tracing  which  includes  all  of  the  relevant  detail,  2)  devising  a 
reference  coordinate  system  for  quantitative  measurements  which  can  be 
precisely  defined  in  terms  of  stable  landmarks  readily  observable  on  the 
film,  3)  devising,  from  a  knowledge  of  the  structural  anatomy,  a  set  of 
transformations  by  which  the  acoustically-relevant  cross-sectional  area 
may  be  inferred  from  the  lateral  measurements,  and  4)  calculating  an 
acoustic  spectrum  from  a  specification  of  the  cross-sectional  area  as  a 
function  of  distance  along  the  vocal  tract  midline  .  The  programs  that  have 
been  written  to  aid  in  the  construction  and  testing  of  the  models  of  speech 
production  correspond  to  parts  3)  and  4)  above. 

A  preliminary  version  of  the  midsagittal-to -area  function  transformation 
program  is  now  in  operation.  With  this  program,  we  can  preset  the  cross - 
dimension-to-area  function  transformation  for  a  particular  talker,  either 
by  reading  in  data  tables  prepared  from  anatomical  measurements  of  him 
or  by  specifying  certain  pararneters  which  determine  an  analytic  descrip¬ 
tion  of  the  transformation.  One  of  the  objectives  of  the  present  study  is 
the  investigation  of  various  possible  types  of  analytic  descriptions.  Once 
the  progra.m  has  been  preset  for  the  talker,  cross -dimension  data  from 
the  films  are  read  in.  The  program  then  calculates  both  the  cross- 
sectional  area  as  a  function  of  distance  along  the  vocal  tract  midline  and 
the  effective  vocal-tract  length. 

We  have  written  a  second  program  which  accepts  cross-sectional  area  func¬ 
tions  derived  by  the  first  program  as  input  and,  with  a  Webster's  equation 
approximation  to  an  acoustic  description  of  the  vocal  tract,  determines  the 
corresponding  acoustic  output  spectra.  We  have  then  compared  a  computed 
spectrum  with  the  spectrum  of  the  actual  sound  produced  at  the  time  of  the 
x-ray  filming.  Such  a  comparison  tests  the  validity  of  an  investigated  model. 
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Simulation  of  the  Human  Laryiix  -  Tnomas  H.  Crystal 

During  the  past  six  months,  we  have  developed  a  set  of  programs  for  sim¬ 
ulation  of  the  behavior  of  a  human  larynx  during  the  productron  of  speech^ 
sounds  and  for  control  and  display  of  the  results  of  the  simu  ation  on 
:::  typewriter  console.  The  programs  still  need  some  addUmnal  modi¬ 
fications ,  but  have  produced  some  interesting  results. 

The  larynx  itself  is  represented  by  two  systems.  A  mechanical  system. 
Lch  n  odels  the  mass  and  tension  of  the  vocal  cords,  is  used  o  de  e  - 
mine,  as  a  function  of  time,  the  opening  between  the  vocal  cor  s  and  ai 
„„„  vocal  coeds.  The  paeac.ce.  ol  ihis  sv=.=m  ^ 

,h.  vocal  cord  opaainp  as  de.otmined  by  th.  n.ecb.mcal  sysl.m.  D 
eace  e,ua.lons  Iha.  describe  these  sys.enis  are  Iterated  successively  to 

1  b-  Tn  addition  the  program  contains  provisions  foi 

produce  the  simulation.  In  addition,  me  p  b 

^  f  of  a  number  of  vocal  tract  and  sub-glottal  models  to 

the  connection  of  any  of  a  numbei 

the  larynit  oiodel.  These  systems  should  together  provide  a  complete  p 

ture  of  larynx  activity. 

A  feature  of  the  simulation  programming  is  the  ease  with  which  the  p^g- 
ress  of  simulation  may  be  controlled  and  its  parameters  varied.  T 
progress,  together  with  various  plotting,  printing,  and  freque  y  y 

routines,  is  controlled  through  simple,  one-word  commands.  Thu  . 

Lay  perform  a  simulation,  examine  the  various  waveforms  generated,  p  i- 
form  Fourier  analyses  on  any  particularly  interesting  waveforms  cha  g 
the  parameters  controlling  the  simulation  dependent  on  the  results  ole 
previous  simulation,  and  simulate  again.  Because  of  extensive  con  rol 
facilities  programmed  into  the  model,  full  advantage  of  time-shared  com¬ 
putation  can  be  obtained  in  this  simulation  problem. 
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Articulatory  Events  of  Speech  Generation  -  William  L.  Henke 

This  project  applies  computer  information-processing  techniques  to  the 
study  of  the  articulatory  events  of  human  speech  generation.  Our  data 
which  consist  of  cineradiograph  films  of  human  speakers,  together  with 
their  corresponding  phonemic  input  and  acoustic  output,  are  studied  and 
modeled  mechanically  for  the  purpose  of  understanding  the  human  speech 
process,  in  particular  the  transformation  between  a  discrete  phonemic  in¬ 
put  and  a  continuous  acoustic  output.  As  we  are  dealing  with  spacially- 
distributed  structures,  we  require  the  input/output  facilities  afforded  by 
the  ESL  Display  Console.  Since  the  information-retrieval  and  modeling 
of  these  data  depend  upon  techniques  of  on-line  interaction  of  a  user  with 
a  computer  and  non-textual  types  of  information,  part  of  our  work  is  the 
development  of  such  techniques.  Three  classes  of  programs  have  been 
developed  to  the  present:  1)  cineradiographic  graphical  data  input,  dis¬ 
play,  and  analysis;  2)  implementation  and  display  of  a  model  of  tongue 
action;  and  3)  a  very  general  solution  for  the  acoustic  transfer  function  of 
the  vocal  tract,  with  associated  displays  and  inputs.  They  are  currently 
disjointed,  but  further  work  will  increase  the  communication  between  them 
so  that  they  may  be  used  singly  or  together,  depending  upon  the  needs  of  the 
aspect  of  speech  being  studied. 

One  part  of  our  work  involves  the  input,  storage,  and  display  of  a  graphical 
representation  of  the  mid-sagittal  plane  of  a  human  vocal  tract  in  action. 

The  MIT  Speech  Group  has  recently  made  frame-by-frame  tracings  of  pro¬ 
jections  of  cineradiographic  films  of  speakers:  These  tracings  are  fortu¬ 
nately  about  the  size  of  the  ESL  Display  Console  scope  face,  so  that  with 
tracing  paper  affixed  to  the  face  of  the  scope,  the  light  pen  can  see  the 
tracking  cross  through  tracing  paper.  We  enter  points  and  connected  line 
segments  for  various  anatomical  structures  and  landmarks,  together  with 
Identifying  information,  into  a  dynamic  data  structure.  At  the  end  of  the  in¬ 
put  for  one  frame,  the  data  are  edited  and  written  out  on  the  disk.  The  first 
generation  display  program  searches  the  file  for  desired  frames  and  si¬ 
multaneously  displays  up  to  three  at  different  intensity  levels.  The  organ¬ 
ization  of  this  file  and  the  specification  and  recording  of  complementary 
acoustic  and  phonemic  data,  which  wiF  be  suited  to  analysis  by  visual  com¬ 
parison,  has  prompted  our  studies  on  graphical  input. 
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We  are  developing  a  model  of  articulatory  action  which  will  distinguish  be¬ 
tween  operators  that  are  associated  with  inputs  representing  linguistic 
categories,  and  states  that  are  associated  with  the  responding  articulatory 
structures.  Considering  the  behavior  of  the  tongue  as  a  responding  struc¬ 
ture,  our  work  develops  general  subprograms  for  the  display  of  graphical 
and  textual  data  and  constructs  a  time-domain  simulation  model  which 
stops  at  any  desired  time  and  displays  the  present  state.  With  this  model, 
a  user  can  then  compare  the  state  with  x-ray  tracings  and  modify  some  of 
the  parameters  of  the  model. 

For  acoustic  comparison,  an  algorithm  which  calculates  the  frequency- 
domain  response  of  the  vocal  tract  is  desirable.  We  have  written  one  for 
general  situations  w'hich  calculates  both  the  complex  spectrum  from  input 
data  of  the  area  function  and  the  loss  as  a  function  of  frequency,  position, 
and  area.  From  the  complex  spectrum,  it  further  calculates  poles  of  the 
transfer  function.  Since  all  this  information  is  displayed  on  a  CP,T,  the 
volume  velocity  distribution  at  any  frequency  is  available.  Programs 
developed  by  other  members  of  the  Speech  Group  to  obtain  area  functions 
from  mid-sagittal  plane  representations  of  the  vocal  tract  will  combine 
with  the  programs  discussed  above, 

Grapheme-t'o-Phoneme  Translation  of  English  -  Francis  F.  Lee, 

It  appears  that  when  an  English-speaking  person  "sounds  out"  a  printed 
English  w'ord  or  an  English-like  non-word,  be  makes  use  of  sev^eral  often 
mutually-supplementary  processes,  namely,  recognition  of  the  pattern  as 
a  whole,  identification  of  certain  affixes,  in  particular  the  suffixes,  iden¬ 
tification  of  compound  words  (e.  g.  Whitestone  as  white-stone,  not  whitest- 
one)  and  a  mapping  of  graphemes  to  phonemes  employing  the  graphemical 
context,  with  the  last  process  occurring  at  the  latest  or  innermost  level. 

The  idea  of  a  mechanical  procedure,  translating  a  printed  English  word 
into  its  pl^onemic  representation  using  oiily  a  single  and  small  exception 
list,  runs  into  immediate  trouble.  By  using  multiple  exception  lists,  con¬ 
sulted  when  necessary,  the  total  number  of  exceptions  and  the  list  of  con¬ 
text  sensitive  rules  can  be  substantially  reduced.  The  extent  of  the  neces¬ 
sity  of  affix  identification  is  being  studied. 
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It  appears  that  the  presence  of  mute  e  in  medial  position  due  to  word  com 
pounding,  such  as  spacec r aft  can  only  be  detected  by  some  sort  of  table- 
lookup  procedure. 

The  18,  000  more  often  used  English  words  (from  Thorndike  and  Lorge)  hav 
been  processed  to  show  the  grapheme-to-phoneme  cluster  relationship  in 
context.  This  body  of  data  has  been  used  for  the  creation  of 
exception  lists.  It  will  be  used  for  the  decomposition  of  complex  clusters 
and  the  formation  of  rules  for  the  pr eviously-mentioned  innermost  level. 

Programming  Support  and  Development  -  Martha  M.  Pennell 

The  Computation  Group  has  provided  both  programming  support  and  devel¬ 
opment  for  users  of  the  time-sharing  system.  Our  program  QUIXOT, 
written  to  alleviate  the  programming  load  of  the  Computation  Group,  has 
been  operational  for  almost  an  academic  year.  A  user  need  only  type  an 
algebraic  expression  for  the  coefficients,  and  QUIXOT  writes,  compiles, 
loads,  and  runs  a  program  to  find  roots  of  a  polynomial  of  degree  twenty 
or  less. 

Encouraged  by  the  success  of  this  program,  the  Computation  Group  modi¬ 
fied  QUIXOT  so  that  the  user  can  typo  in  his  coefficients  as  one  expression, 
rather  than  as  real  and  imaginary  parts.  Because  MAD,  the  source  lan¬ 
guage  for  this  program,  does  not  provide  for  complex  arithmetic,  a  sub¬ 
routine  call  must  perform  each  complex  arithmetic  operation.  Also,  we 
incorporated  a  compiler -writing  program  into  the  QUIXOT  program,  which 
saves  the  user  some  initial  algebra.  An  outgrowth  of  this  program  is 
COMCAL,  a  routine  to  evaluate  an  expression  in  complex  arithmetic. 
Similar  in  operation  to  QUIXOT,  it  first  queries  the  user  for  his  symbols, 
parameter  names,  and  algebraic  expression;  then,  with  the  command¬ 
programming  feature  of  time-sharing,  it  writes,  compiles  and  executes 
a  program  to  evaluate  the  expression. 
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RLE  Statistics  on  CTSS  -  Martha  M.  Pennell 

I  cite  examples  which  might  give  some  idea  of  the  time  difference  involved 
in  a  programmer's  speed  with  CTSS  versus  batch  processing.  They  are 
taken  from  the  work  of  RLE  Computing  Group. 

A  mathematics  niajor  with  one  summer's  FORTRAN  experience  was  hired 
as  a  scientific  programmer.  She  was  given  a  short  demonstration  of  CTSS 
-  compute  A  from  the  formula  A  =  B  +  1  where  B  is  input  from  the  tele¬ 
type  —  and  told  to  duplicate  it  where  B  is  an  array  of  N  numbers  to  be 
given  at  execution  time.  She  did  so  in  about  four  hours.  The  next  day 
this  person  was  given  an  RLE  problem,  and  the  calling  sequences  for  the 
two  necessary  subroutines  (a  Bessel  function  and  roots  of  a  polynomial), 
and  she  was  told  to  code  it  in  MADTRN.  Within  three  days  her  program  pro¬ 
duced  results.  At  first  it  contained  errors,  but  she  had  set  up  the  subroutine 
linkages  correctly  and  had  made  no  errors  in  its  logic.  Her  errors  lay  in 
the  algebra  for  computing  the  sum.  This  initial  coding  was  done  with  a 
minimum  of  supervision.  When  I  later  asked  how  long  she  thought  it  would 
have  taken  her  on  batch  processing,  she  said  a  month,  because  her  errors 
were  due  to  her  unfamiliarity  with  FORTRAN  and  a  series  of  programming 
blunders,  such  as  keypunching  errors  and  omission  of  control  cards  or 
subroutine  decks. 

Another  user  has  five  years  of  experience  in  FORTRAN  and  one  year  with 
MADTRN.  The  majority  of  programs  she  has  had  to  write  on  CTSS  can  be 
coded  and  debugged,  on  the  average,  with  a  two-hour  session  at  the  console. 
These  same  programs,  if  run  and  checked  out  under  batch  processing  here 
at  M.  I.  T.  ,  would  take  at  least  four  to  six  days.  These  programs  involve  a 
main  program  and  possibly  two  or  three  subroutines  of  about  40-50  state¬ 
ments  each.  These  programs  have  found  roots  to  polynomials,  evaluated 
Bessel  functions,  performed  numerical  quadrature,  and  so  on.  She 
finds  that  with  programs  of  this  length,  three  or  four  tries  usually  debugs 
them.  If  you  take  into  account  machine  errors,  bad  tapes,  lost  output,  this 
usually  means  about  a  week  in  a  batch-processing  scheme. 

This  person  also  decided  to  learn  MAD.  The  problem  she  chose  to  code  was 
a  tree  sort.  Given  the  people  in  her  office,  their  sex,  their  marital  status. 
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their  laboratory  affiliation,  and  their  telephone  number,  the  problem  was  to 
to  write  a  program  to  sort  them,  grouped  alphabetically,  according  to  one 
or  more  of  their  characteristics.  The  MAD  program  was  about  40-50 
statements  long  and  took  six  tries  before  it  finally  compiled,  because  the 
programmer  was  not  completely  familiar  with  the  punctuation  and  rules  of 
syntax.  Two  more  tries  debugged  the  program  logic.  This  work  was  done 
in  about  two  hours  of  console  use. 

Sorting  of  Personnel  Records  -  Robert  L.  Rappaport 

The  RLE  personnel  roster  includes  approximately  700  people,  who  are 
classified  by  several  categories;  position,  research  group,  rooms  occu¬ 
pied,  telephone  extension,  and  so  on.  For  example,  an  associate  pro¬ 
fessor  of  the  Electrical  Engineering  Department  who  works  in  the  Plasma 
Physics  Group  must  be  listed  under  each  of  these  files.  At  present,  these 
files  are  brought  up  to  date  and  sorted  manually. 

We  have  written  a  computer  program  to  automate  this  organization.  The 
program  operates  on  a  masterfile  which  contains  all  the  essential  informa¬ 
tion  for  every  member  of  RLE,  so  that  we  need  make  correction  only  once. 
Sorting  can  be  done  on  any  of  the  items  entered,  A  second  program  prints 
an  alphabetical  directory  of  all  RLE  personnel.  We  have  checked  out  both 
of  these  programs  for  a  limited  number  of  names  and  will  shortly  extend 
them  to  the  entire  RLE  roster. 
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COMIT  -  Victor  H.  Yngve 

During  the  year  we  maintained  the  COMIT  system  as  a  standard  program¬ 
ming  language  available  to  users  of  MAC,  and  we  further  developed  the 
language  so  that  a  new  system,  COMIT  II,  will  soon  be  released  for  use 
in  MAC.  COMIT  II  includes  many  new  features  and  conveniences  for  the 
programmer,  while  any  program  that  works  in  COMIT  will  also  work  in 
COMIT  II.  We  also  developed  a  program  VEDIT,  which  is  written  in  COMIT 
and  is  especially  useful  in  the  writing  and  editing  of  COMIT  programs.  It 
introduced  a  number  of  new  features  that  have  found  their  way  into  recent 
edit  programs. 

Jared  Dariington  further  developed  his  set  of  COMIT  programs  that  check 
the  validity  of  arguments  expressed  in  English.  These  programs  are 
unique  in  that  they  incorporate  several  levels  of  analysis  of  the  input  argu¬ 
ment  to  determine  the  lowest  level  of  logical  analysis  sufficient  to  prove 
the  argument  valid.  In  addition,  a  program  in  COMIT  which  partially 
translates  programs  from  FAP  to  GE635  GEM  has  been  developed  by  two 
undergraduate  students. 
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Computer-Aided  Design  -  Charles  Garman,  Richard  Bayles,  and 

Harold  Levin 

During  the  past  year,  Mr.  Bayles  and  Mr.  Levin  were  programming  the 
section  of  the  CTSS  supervisor  responsible  for  maintaining  the  display 
on  the  ESL  Display  Console.  Since  the  design  of  the  software  for  the 
display  system  requires  a  fairly  rigid  and  highly  artificial  structuring 
of  display  data,  operation  of  the  first  version  of  the  supervisory  program 
soon  showed  defects  in  this  design. 

We  began  in  August,  1964,  to  redesign  the  interface  between  the  user 
and  the  display  system  and  succeeded  in  developing  a  new  supervisory 
program.  Among  its  features  were:  the  recognition  of  the  slave  display, 
which  had  been  added  to  the  display  console,  so  that  two  users  can  main¬ 
tain  different  displays  simultaneously;  the  addition  of  a  system  which 
assigns  real-time  tasks  to  various  analog  and  digital  inputs  of  the  display 
console  and  to  different  sections  of  the  image;  and  the  addition  of  sub¬ 
routine  pictures,  which  facilitate  display  of  similar  graphical  objects. 

This  project  was  completed  in  December  ,  and,  with  the  help  of  Michael 
Brescia  of  the  ESL  Display  Group,  we  have  continued  to  modify  and 
improve  the  system. 

At  the  end  of  December  Mr.  Bayles  began  programming  for  the  Disk- 
Drum  Strategy  module  that  was  to  be  incorporated  into  the  revised  File 
Control  System  of  the  CTSS  supervisor.  This  work  was  continued  and 
has  since  been  completed  by  Mr.  Garman.  He  was  also  involved  during 
the  previous  summer  in  programming  a  system  that  generated  and  dis¬ 
played  on  the  ESL  console  that  were  surfaces  derived  from  an  algorithm 
of  Professor  S.  A.  Coons.  These  surfaces  are  completely  general  in 
nature  and  may  model,  with  remarkable  precision,  objects  as  diverse  as 
boat  hulls,  an  octant  of  a  sphere,  adjacent  planes  of  a  regular  tetrahedron, 
or  a  section  of  an  airplane  fuselage.  In  December,  he  began  work  on 
modifications  and  extensions  to  the  CTSS  'ED*  command,  which  is  a  con¬ 
text  editor  for  BCD  card-image  files  and  he  is  reprogramming  the  input- 
output  section  of  'ED'  to  conform  to  the  new  file-control  system. 
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On-Line  Mathematical  Analysis  -  Roy  Kaplow.  Stephen  L.  Strong,  and 

John  W.  Brackett 

With  our  MAP  ,  user  with  little  knowledge  of  computers  or  pro¬ 

gramming  can  solve  a  wide  variety  of  mathematical  problems  expressed 
in  normal  mathematical  terminology,  MAP  contains  many  interrelated 
programs,  only  a  small  fraction  of  which  may  be  executed  in  response  to 
any  one  request.  The  language  consist,  of  functional  equations  and  an 
expandable  group  of  commands.  All  communications  can  be  transmitted 
through  a  standard  Project  MAC  console. 

The  system  evaluates  one-dimensional  functional  equations  that  involve 
constants,  vectors,  library  functions  and  ordinary  arithmetic  operators 
The  form  used  is  identical  to  normal  notation,  except  that  multiplications 
only  implied  in  the  first  form  must  be  specified  in  the  latter  and 
functional  operations,  such  as  sine  and  cosine,  must  be  designated  by  a 
terminating  f.  MAP  automatically  generates  the  q(x)  vector  at  equal 
mtervals  over  any  range  in  x  specified  by  the  user.  If  the  range  and 
interval  in  x,  the  values  of  the  constant  a,  or  the  values  of  the  vectors  y  (x) 
are  not  already  specified,  they  are  requested  by  MAP  and  may  then  be 
typed  in  under  the  control  of  context  editor. 

Commands  in  MAP  correspond  to  requests  for  operations  upon  functions 
such  as  q(x),  or  upon  groups  of  functions,  which  cannot  proceed  in  the 
point-by-point  fashion  of  equations.  With  such  commands  at  his  disposal 
the  user  need  not  remember  the  parameters  required  for  the  operation.  ’ 

As  an  example,  the  typing  of  the  single  word  "integrate"  will  cause  MAP 
to:  ask  the  user  if  he  desires  to  integrate  between  two  constants,  the  full 
range  of  an  independent  variable,  or  some  combination  of  constant  and 
variable  limits;  ask  the  name  of  the  function  to  be  integrated:  the  name  to 
be  given  to  the  constant  or  function  answer;  and.  if  necessary,  the  values 

already  defined,  its  values  will  be  requested  and  may  be  typed  in  with  a 
context  editor.  The  sophisticated  user  may  prescribe  all  of  the  necessary 
in  ormation,  regarding  options  within  the  operation  and  parameter  names 
or  va  ues,  and  avoid  unnecessary  communication.  The  following  operations 
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are  now  available:  integrate,  differentiate,  convolute,  Fourier  trans¬ 
form,  change  of  basis,  and  least-square  fitting. 

In  addition  to  the  equations  and  simple  commands  discussed  above,  other 
commands  are  available  which  facilitate  problem  solving  and  which  extend 
the  range  of  MAP.  These  are  briefly  described  below  under  five  separate 
classifications. 

A.  DATA  HANDLING 

These  commands  provide  for  off-line  input  of  blocks  of  data,  for  the 
deletion  of  data  from  mass  memory,  and  for  the  saving  of  data  and  pro¬ 
grams.  All  information  is  handled  by  its  symbolic  name,  and  the  user 
need  not  be  concerned  with  computer  techniques. 

B.  FUNCTION  HANDLING 

The  "minmax"  command  changes  the  range  and  interval  for  which  an 
independent  variable  or  function  is  defined.  The  "select"  command 
manipulates  a  portion  of  an  individual  function. 

C.  OUTPUT 

The  "print"  command  is  used  for  all  output  other  than  system-generated 
messages  and  questions.  The  format  of  printing  may  be  specified  by  the 
request  of  a  user.  When  he  so  desires,  a  user  may  request  that  only 
part  of  the  tabulated  range  of  a  function  be  printed,  together  with  an 
interval  in  the  independent  variable. 

D.  COMMAND  SEQUENCE  CREATION  AND  EXECUTION 

The  "run",  "create",  and  "edit"  requests  define  a  sequence  of  MAP  state¬ 
ments  and  store  them  in  mass  memory  for  subsequent  use.  "Create" 
allows  typing  in  a  sequence  of  commands  under  control  of  the  context 
editor.  The  sequence  is  assigned  a  name  and  executed  by  "run"  and  the 
assigned  name.  "Edit"  edits  a  sequence  of  commands  to  correct  errors 
or  to  substitute  variable  names  when  the  sequence  is  executed  by  a  "run" 
command.  A  user  who  has  not  learned  a  complex  language,  such  as 
FORTRAN  or  MAD,  can  thus  create  a  programming  structure  for  his 
specific  application  and  call  it  by  its  symbolic  name. 
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E.  EXECUTION  OF  MAD  PROGRAMS 

With  "Execute,  "  a  user  can  write  his  own  command  in  the  MAD  language 
and  use  it  within  the  MAP  system.  The  existence  of  this  command  makes 
the  system  easily  expandable  and  should  encourage  users  to  add  to  the 
system.  After  the  user  has  written  a  MAD  program,  such  as  "root,  "  he 
need  only  type  "execute  root"  to  run  his  program.  As  an  aid  in  the 
creation  of  his  own  commands,  a  large  number  of  system  subroutines 
are  available  to  the  facilities  of  the  MAP  system.  Other  subroutines 
accomplish  such  tasks  as;  loading  any  function,  previously  defined 
within  the  MAP  system,  from  mass  memory  into  any  array  specified  in 
the  user's  program;  obtaining  one  or  more  values  of  a  function  from  mass 
memory,  and  storing  them  as  the  values  of  any  variable  in  the  MAD  pro¬ 
gram;  creating  a  MAP  function  from  an  array  in  the  user's  program:  and 
executing  any  CT.SS  command.  All  of  the  necessary  subroutines  are 
automatically  loaded  along  with  the  user's  program  when  the  "execute" 
command  is  given. 

During  the  past  year  the  MAP  system  was  made  available  to  a  number  of 
users  of  the  MIT  Computation  Center.  Many  of  the  alterations  and 
additions  provided  during  this  year  are  the  result  of  our  experience  with 
those  users.  Plans  for  the  forthcoming  year  include  extension  of  the 
system  to  functions  of  more  than  one  independent  variable,  matrix  opera¬ 
tions,  differential  equations,  and  addition  of  graphical  input  and  output. 

Derivation  of  Preliminary  Ship  Lines  -  Theodore  M.  Pitidis -Poutous 

The  object  of  this  investigation,  supported  by  Bureau  of  Ships  Contract 
No.  NObs  90100,  was  to  develop  a  suitable  procedure  for  the  mathematical 
derivation  of  preliminary  ship  lines,  for  a  new  design,  by  reference  to 
the  lines  of  a  parent  ship.  An  algorithm  has  been  formulated  to  determine 
a  new  curve  or  surface,  in  the  form  of  a  polynomial,  resembling  as 
closely  as  possible  another  graphically- specified  curve  or  surface,  and 
satisfying  a  prescribed  set  of  constraints.  The  computat'ons  involved  are 
facilitated  by  the  use  of  special  polynomials. 


Various  ways  of  applying  this  algorithm  to  the  delineation  of  preliminary 
ship  lines  have  been  investigated.  In  particular,  a  procedure  has  been 
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in  the  region  of  the  bow,  so  that  the  corresponding  parent  section  is 
Ignored.  Further  investigation  concerning  the  use  of  the  constraints 
already  employed,  or  possible  additional  constraints,  will  be  necessary 
to  critically  evaluate  this  procedure. 


These  remarks  apply  to  the  trawler  case  examined  so  far,  as  full  ship 
forms  have  not  been  investigated.  Preliminary  experience  with  very 
full  sections  indicate  that  a  few  more  terms  would  be  required  for  approxi¬ 
mations  comparable  to  those  obtained  for  the  trawler.  Further  investi¬ 
gation  IS  necessary,  however,  before  definite  conclusions  can  be  drawn 
with  respect  to  full  forms.  (See  Pitidis-Poutous.  Appendix  C.) 


Time-Sharing  Reactor  Code  System  -  Kent  F.  Hansen.  Ian  C.  Pyle, 

Billy  V.  Koen  and  Q.  L.  Miller 

The  purpose  of  this  work  is  to  study  the  use  of  time-sharing  in  Nuclear 
Engineering  education  and  research.  In  a  question-and-answer  session 
at  the  console,  a  user  derives  input  data  for  any  nuclear  code,  in  such  a 
way  that  he  need  have  no  knowledge  of  computer  programming  at  all.  In 
the  past  year,  we  have  completed  a  general-purpose  input/edit  program 
that  obtains  input  data  from  the  console  and  numerous  nuclear  codes. 

The  entire  collection  of  programs  is  called  the  TREC  system  (Time- 
Sharing  Reactor  Codes).  ”” 


The  TREC  system  operates  as  follows:  when  the  user  requests  a  particu¬ 
lar  code,  TREC  obtains  the  appropriate  files  and  certain  satellite  routines 
that  read  in  the  input  data.  The  input/edit  program  asks  the  user  for  each 
Item  of  input  in  a  language  natural  to  the  problem.  The  user  may  review, 
change,  or  correct  his  data.  After  all  data  is  properly  given,  he  may 
request  the  code  to  start;  TREC  then  transfers  data  and  control  to  the 
nuclear  code.  At  the  conclusion  of  the  calculation  or  after  an  interrupt,' 
control  is  returned  to  the  input/edit  program  for  possible  changes  in  data. 
At  the  conclusion  of  the  run,  the  user  may  select  another  code. 


The  function  of  the  routine  TREC  is  to  ask  the  user  what  code  he  wants, 
determine  if  the  code  i,s  available  or  not,  and,  if  it  is  available,  locate 
the  symbol  table  and  test  routine  for  the  code.  The  symbol  table  itself 
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is  a  list  of  encoded  variable  names  and  descriptive  phrases  that  define 
each  variable.  The  data  is  stored  in  a  linear  array  for  delivery  to  the 
nuclear  code  itself.  The  test  routines  are  used  to  check  all  input  data. 
Occasionally  a  block  of  data  may  be  standard  information,  such  as  neutron 
delayed -yield  fractions  for  fissionable  species,  in  which  case  the  test 
routine  also  provides  the  standard  data.  Upon  completion  of  the  input/ 
edit  phase,  the  input  data  array  is  written  out  in  a  pseudotape  for  delivery 
to  the  nuclear  code.  Within  the  symbol  table,  the  name  of  the  file  con¬ 
taining  the  nuclear  code  is  stored  for  transfer  of  control. 

Finally,  the  routine  QUIN  is  loaded.  QUIN  scans  the  symbol  cable  and 
prints  comments  and  questions  on  the  console.  Thus,  the  symbol  table 
contains  the  names  of  all  input  data  along  with  additional  details.  The 
test  routine  is  used  to  test  the  validity  of  each  piece  of  input  data.  If  an 
impossible  or  inappropriate  value  is  given,  the  test  routine  prints  a 
comment  to  that  effect.  Such  details  as  decimal-to-binary  conversion 
are  also  contained  in  QUIN. 

Several  nuclear  codes  have  been  prepared  and  are  now  routinely  available. 

A  brief  description  of  each  is  given  below. 

1.  KINET  solves  the  infinite  medium-kinetics  equations  with  a 
variable  number  of  delayed  neutron  groups. 

2.  DIFFUS  is  a  one-dimensional,  multigroup-diffusion  code. 

Its  program  handles  up  to  16  groups,  5  regions  and  5  isotopes  per  region, 
and  also  a  maximum  of  100  space  points. 

3.  CROSEC  is  a  code  that  computes  group-average,  nuclear 
cross -s  ections . 

4.  The  Cross  Section  Library  is  not  a  code  but  rather  a  large 
block  of  nuclear  cross  sections  of  many  isotopes.  The  data  is  encoded  in 
a  form  used  by  CROSEC. 

5.  DNPORT  is  a  double  transport  code.  This  code  analytically 
solves  the  transport  equation  for  one-group  slabs.  Although  it  is  written 
and  debugged,  the  symbol  table  is  not  yet  completed. 
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:^„Stress-AnaIvsis  Program  -  Richard  P.  Parmelee 

During  the  past  year,  I  finished  and  tested  a  two-dimensional  stress - 
analysis  algorithm  for  thin  domes  or  shells,  and  then  extended  the 
algorithm  to  handle  three-dimensional  stress  analysis  of  sculptured  bodies 
such  as  crankshafts,  I  completed  this  three-dimensional  program  in 
March  of  1965  and  since  then  have  solved  test  problems. 

The  object  of  my  research  to  find  a  method  for  stress  analysis  of 
virtually  any  geometric  region,  and  then  implement  this  method  with  a 
computer  program  that  will  permit  a  user  who  knows  nothing  about  com¬ 
puters  to  define  his  problem  and  obtain  accurate  answers.  Although  I 
found  an  algorithm  suitable  for  analyzing  a  large  class  of  problems,  I 
had  problems  with  the  implementation  of  this  algorithm.  The  availlbility 
of  a  time-shared  computer  made  possible  the  coding  and  debugging  of  this 

program  in  a  reasonable  time,  as  well  as  the  guiding  and  error  checking 
of  the  user  input. 

Although  I  did  some  graphical  displa^.work  with  the  ESL  Display  Console, 

the  input,  output,  and  program  control  was,  for  the  most  part,  restricted 
to  a  teletype. 


Stress -Analysis  Coniormal  Mapping  -  Paul  A 


W  ieselmann 


Fundamental  to  the  solution  of  plane-stress -analysis  problems  by  comple: 
variable  methods  is  the  generation  of  a  function  to  conformally  map  the 
unit  circle,  or  upper  infinite  half-plane,  onto  a  region  described  by  the 

physmal  domain  of  the  problem.  In  general,  only  approximations  to 

physical-domain  geometry  can  be  obtained,  such  as  an  approximation  by 
a  polygon.  ^ 

My  work  in  the  past  year  has  attempted  to  implement  analytic  methods  to 
obtain  the  mapping  function,  in  particular  the  Schwarz-Christoffel  integral 
that  n,aps  the  unit  circle  onto  the  interior  of  an  n-sided  polygonal  domain, 
w  1C  IS  restricted  to  be  convex.  Riemann's  mapping  theorem  applied  to 
the  Schwarz-Christoffel  integral  allows  three  of  the  n+z  coefficients  to  be 
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fixed  arbitrarily.  The  program  first  estimates  values  for  the  remaining 
n-1  coefficients,  determines  an  improvement,  and  continues  until  it  attains 
a  predetermined  accuracy.  At  present,  this  interactive  scheme  is  com¬ 
pletely  free  of  manual  intervention,  except  in  a  few  special  cases  which 
require  much  careful  and  tedious  intervention. 

I  shall  free  the  iteration  scheme  of  manual  intervention  and  make  complex- 
plane  numerical-integration  techniques,  especially  in  the  vicinity  of  a 
pole,  more  accurate  in  the  evaluation  of  the  coefficients.  A  result  of 
this  increased  accuracy  will  be  an  improvement  in  the  iteration  scheme, 
since  one  of  its  present  failures  is  due  in  part  to  an  accumulation  of  error 
in  the  determination  of  the  length  of  the  boundary. 

Computer-Aided  Teaching  of  Dynamic  Systems  Behavior  -  Ronald  C. 

Rosenberg 

The  three  main  aspects  of  our  work  in  the  past  year  were:  the  continued 
development  of  the  ENPORT  system  for  use  as  a  simulated  dynamic - 
system  laboratory,  the  extension  of  part  of  its  programming  algorithm  to 
a  hand  procedure  for  student  use  in  analysis,  and  the  introduction  of  a 
simulated  blank  box  as  a  significant  teaching  tool. 

The  ENPORT  system  can  simulate  the  behavior  of  linear,  lumped-parameter 
dynamic  systems  up  to  fifteenth,  order.  Its  input/output  request  forms  have 
been  simplified,  but  have  been  made  more  flexible  by  improved  program¬ 
ming  techniques.  In  addition,  a  local  digital  memory  to  the  teletype- 
station-associated  analog  computer,  which  are  automatically  setup  by  the 
ENPORT  program,  'provide  a  s ignal-versus -time  oscilloscope  display  that 
represents  the  dynamic  system  behavior. 

The  analysis  procedure  starts  with  a  bond  graph  model,  assig'  computing 
causality  and  a  sign  convention,  and  selects  state  variables.  The  system 
equations  may  then  be  written  from  the  graph  in  an  organized  form  and 
reduced  to  a  standard  form  algorithmically.  By  means  of  a  black  box,  a 
user  can  specify  the  input  and  observe  the  output  in  order  to  discover  the 
unknown  system.  -  The  system  has  been  successfully  operated  as  an  on-line 
aid  to  classroom  teaching  in  a  class  given  by  Prof.  H.  M.  Paynter. 
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Automatic  Network  Synthesis  -  Marc  P.  Weinbtreer 

In  the  past  year  I  attempted  to  find  good  search  strategies  in  a  parameter 
space.  A  strategy  achieves  optimum  performance  through  a  choice  of 
several  possible  network  structures  followed  by  a  search  for  a  best 
parameter  combination.  It  is  called  optimum  if  it  requires  the  fewest 
trials  or  the  least  time  on  a  computer.  I  found  a  strategy  which  is 
optimum  under  certain  conditions  and  .-often  substandard  in  other  cases 
where  the  required  conditions  are  not  met.  Instead  of  conducting  a  single 
series  of  trials  in  the  entire  search  space,  this  procedure  conducts  a 
multi-stage  search  in  ever-narrowing  subspaces.  Simple  rules  indicate 
the  allocation  of  the  entire  search  to  successive  stages.  This  procedure 
leads  to  a  large  decrease  in  the  number  of  trials  and  is  more  efficient 
than  g radient -type  methods  or  other  deterministic  strategies  in  several 
nontrivial  cases. 

I  also  studied  the  effect  of  noise  in  the  measurements  or  calculations. 

Since  the  random  search  appears  to  be  less  affected  by  these  errors  than 
gradient-like  methods  in  several  instances,  it  has  a  broad  field  of  applica¬ 
bility.  The  designer  must  be  able  to  override  the  automatic  search  pro¬ 
gram,  whenever  he  wants  to  try  out  his  own  ideas.  Indeed,  he  learns 
more  about  his  networks  during  the  early  stages  of  design.  Since  the 
randomness  of  the  chosen  trial  points  often  yields  good  results  from  the 
very  beginning,  the  designer  may  want  to  cut  off  the  automatic  routine  in 
order  to  investigate  more  closely  the  excellent  early  returns.  Thus,  an 
interaction  between  man  and  machine,  such  as  the  one  afforded  by  time¬ 
sharing,  is  a  very  important  part  of  this  method. 
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A  Computer  Model  of  Kinesthesis  -  George  Piotrowski 


The  purpose  of  this  work  is  to  experiment  with  the  use  of  the  kinesthetic 
in  communicating  between  man  and  computer.  The  model  of  kinesthesis 
will  be  a  radially-extensible  pivoting  arm  which  has  digital  shaft  encoders 
for  position  measurement  and  computer-controlled  torque  motors  to  pro¬ 
vide  a  force  output.  The  arm  will  cover  a  field  of  about  fifteen  inches 
square  and  be  able  to  exert  about  five  pounds  of  force.  The  first  version 
of  the  arm  will  be  as  simple  as  possible  and  have  most  of  the  computa¬ 
tional  functions  carried  out  by  software  instead  of  hardware.  Thus  far, 

I  have  written  a  program  to  simulate  the  device  and  also  a  control  program 
for  the  device;  and  I  am  using  the  two  to  ascertain  the  effect  of  several 
design  parameters,  notably  the  sampling  rate  and  the  accuracy  of  measure¬ 
ment,  on  overall  performance.  With  the  control  program,  a  user  can 
select  various  mechanical  systems  for  simulation.  This  program  has 
been  written  in  a  general  fashion,  so  that  other  algorithms  can  be  added 
very  easily. 

Computer  Solutions  for  Boundary-Value  Problems  -  Coyt  C.  Tillman 

I  am  developing  a  general,  easy-to-use  system  for  the  solution  of 
dimensional  boundary-value  problems.  The  programs  written  so  far  are  a 
command  processor,  together  with  associated  I/O  and  free-storage 
routines,  boundary  shape  description  routines,  and  a  set  of  mathematical 
routines  for  the  interpretation  of  arithmetic  argument  strings.  The 
command  processor,  though  written  with  boundary-value  problems  in 
mind,  will  be  more  generally  usable.  It  admits  to  completely  format- 
free  input  and,  at  the  same  time,  avoids  the  imposition  of  the  stringent 
punctuation  conventions  of  other  fc^rmat-independent  input  interpreters. 

A  macro  facility  provides  for  the  definition  of  new,  composite  commands, 
the  modification  of  the  system's  vocabulary,  and  so  on.  I  am  now  attend¬ 
ing  to  the  requirements  of  the  data  structure  needed  to  represent  the 
diflerence  equations  by  which  boundary  value  problems  are  being  modelled. 
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Algebraic  Expression  Compiler  -  Grover  C.  Gregory 

I  am  developing  a  compiler  for  algebraic  expressions.  It  is  primarily 
intended  for  use  under  the  supervisory  system  that  Coyt  Tillman  has  des¬ 
cribed  in  this  section. 

The  most  recent  work  has  been  on  the  logic  of  the  compiler's  second  pass. 
This  pass  is  responsible  for  producing  the  appropriate  machine  code, 
using  as  its  input  the  tables  generated  during  the  first  pass.  The  logic 
is  similar  to  that  found  in  the  MAD  compiler  which  interprets  lists  to 
decide  how  to  arrange  the  machine  code.  The  aim  has  been  to  retain  the 
flexibility  of  using  such  limits  while  reducing  the  amount  of  core  they 
require.  To  that  end,  a  different  logic  is  used  to  handle  the  triple  inter¬ 
actions.  In  addition,  the  compiler  is  constructed  so  that  it  sees  variables 
and  unsubscripted  functions  as  identical,  a  requirement  of  the  system 
under  which  it  will  operate.  Testing  of  the  whole  compiler,  as  well  as 
modifications  to  allow  handling  of  subscripted  variables  and  Boolean  con¬ 
ditionals,  will  begin  soon. 

Algorithm  to  Aid  Logic  Design  -  Gerald  A.  Mal-ey 

My  main  interest  at  the  beginning  of  this  project  lay  in  the  development  of 
a  computer  algorithm  for  the  design  of  minimum  Nor-Nand  circuits  from 
Boolean  expressions.  Algorithms  do  exist  in  this  field  for  the  limited 
case  where  complements  of  the  input  variables  are  available;  but,  in  the 
design  of  an  actual  computer,  these  complemented  variables  are  seldom 
available.  When  complements  are  available,  only  two  levels  of  logic  are 
i-equireM;  but,  in  the  general  case,  that  of  no  complements,  three  levels 
are  required.  With  the  available  computer  time.  I  was  able  to  develop 
the  required  algorithm  for  this  type  of  implementation. 

The  proposing  of  an  algorithm  for  a  problem  such  as  this  is  relatively 
simple,  but  the  proving  or  disproving  of  such  an  algorithm  can  take  months 
Oil  a  remote  console  I  was  able  to  propose  an  algorithm  and  run  several 
hundred  cases  in  an  hour  or  so.  With  response  times  of  this  order,  i  was 
able  to  find  delinquencies  and  modify  the  algorithm  in  a  matter  of  days, 
rather  than  months.  Iterations  of  this  kind  finally  produced  a  simple 
algorithm  that  I  would  not  have  uncovered  with  normal  machine  methods. 
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Social  Systems  Analysis  -  Ithiel  D.  Pool 

We  are  currently  working  on  a  generalized  editing  system  for  social 
science  data.  Over  past  decades  the  handling  of  accumulated  large  files 
of  survey  records  and  census  records  have  presented  two  major  problems; 
the  adaptation  of  older  files  to  a  computer,  and  the  presence  of  errors  in 
the  data  that  will  mislead  a  computer  program.  We  are  developing  editing 
routines  for  on-line  revision  and  correction  of  these  older  files.  The 
editing  system  needs  to  compare  data  records  with  the  protocol  of 
allowable  forms  and  provide  feedback  with  respect  to  the  editing  decisions. 
The  first  version  ol  the  editing  system  is  now  being  used,  and,  the  ex¬ 
perience  of  its  users  is  very  helpful  in  the  current  design  phase  of  the 
next  version. 

The  social  science  data  surveys  and  census  records  on  which  we  focus  con¬ 
sists  of  both  a  codebook  and  records  that  describe  individual  people  in  ways 
outlined  in  the  codebook.  The  editing  system  reorganizes  the  record,  and 
a  set  of  programs  handles  the  codebook.  The  codebook  requires  that:  1) 
the  data  be  put  into  computer  readable  form;  2)  parts  of  the  data  be  tagged  ' 
by  people  to  indicate  characteristics  of  its  form  and  content;  and  3)  the 
data  be  checked  for  internal  discrepancies  and  the  user  informed  of  these 
on-line  so  that  errors  may  be  corrected. 

In  the  CONCOM  Project  (ARPA  Contract  No.  920F-9717)  we  try  to  simulate 
the  flow  of  messages  in  a  nation.  We  first  represent  the  audience  by  creat¬ 
ing  a  population  of  hypothetical  indi\’iduals  that  would  constitute  a  fair  sample 
of  a  real  population.  To  describe  our  hypothetical  population  of  individuals, 
we  reconstructed  from  aggregate  statistical  data  a  possible  distribution  of 
facts  about  individuals  compatible  with  the  reported  statistics.  In  this 
description  we  used  a  technique,  developed  by  Frederick  Mosteller,  which 
synthesizes  an  n-dimensional  table  from  any  number  of  m-dimensional 
tables  that  constitute  its  marginals.  For  each  of  the  cells  into  which  the 
population  is  divided,  this  technique  determines  its  frequency  and  several 
of  its  media  characteristics.  With  it  we  can  use  not  only  marginals  but 
also  whatever  estimates  of  interactions  are  available. 
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The  next  part  of  the  simulation  produces  estimates  of  probable  media 
exposure  by  the  persons  in  the  hypothetical  population.  An  audience 
estimate  is  usually  gotten  by  some  kind  of  transformation  of  a  circulation 
or  set  estimate:  for  example,  the  number  of  radio  sets,  the  total  daily 
newspaper  circulation,  the  total  magazine  circulation,  total  book  publi¬ 
cation,  and  so  on.  All  circulation  figure.s  then  need  to  be  multiplied  by 
factors  to  provide  audience  estimates.  In  such  total  audience  estimates, 
we  use  the  Mosteller  technique  to  take  account  of  all  marginal  .s  about 
which  we  have  information  and  all  interactions  about  which  we  have 
information. 

To  make  our  simulation  dynamic,  we  must  take  into  account  differences 
in  audience  habits.  For  example,  a  newspaper  read  by  one  million  people 
each  day  might  conceivably  be  read  by  one  million  people  365  times  a 
year  or  it  might  conceivably  be  read  by  365  million  different  people  only 
once.  The  actual  distribution  lies  somewhere  in  between.  We  refer  to 
this  distribution  of  readers  by  frequency  as  the  concept  of  cumulation. 

We  have  developed  programs  that  estimate  the  underlying  distribution  of 
probabilities  for  any  available  cumulation  data  and  generate  cumulation 
curves  from  assumptions  about  the  structure  of  the  medium. 

Our  simulation  also  accounts  for  duplication:  that  is,  the  incidence  of 
shared  audience  among  different  media.  Exposure  probabilitie .s  for  each 
individual  in  the  population  are  thus  assigned  on  the  basis  of  mean  pioba- 
bihties  for  the  cell  of  the  population  to  which  he  belongs,  the  likely  shape 
ot  the  probability  distribution  as  indicated  b\  cumulation  facts,  and  the 
likely  non -random  duplication  among  media.  In  the  final  part  of  the 
simulation,  a  reader  of  a  particular  medium  will  note  a  particular  news 
story  in  it  with  a  probability  p,  which  need  not  be  constant.  In  the  early 
period  of  a  crisis,  when  the  saliency  of  the  events  is  low,  p  may  be  small: 
later,  it  may  gradually  rise,  p  may  also  vary  with  different  groups  in 
the  population.  With  the  simulation,  we  can  thu.s  formulate  a  set  of 
assumptions  about  the  flow  of  news  in  a  total  society  and  reach  a  con¬ 
clusion  ns  to  who  has  been  reached  and  how  often.  At  present,  we  have 
a  working  model  written  in  FORTRAN  which  we  are  running  on  test  data 
and  we  are  translating  the  entire  set  of  programs  into  AED-O  in  order  to 
run  the  simulation  with  on-line  human  decision  making. 
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With  CRISISCOM  Project  support  from  the  CONCOM  Project,  \ve  have 
developed  a  different  siumlation  of  an  aspect  of  crisis  communications, 
namely  the  CRISISCOM  siumlation  (partially  supported  by  the.Naval 
Research  Laboratory).  Over  the  past  year  we  have  tested  CRISISCOM 
by  representing  the  Kaiser  and  the  Tsar  as  they  receive  about  1500 
messages  covering  events  which  occurred  during  the  week  before  World 
War  1.  On-line,  an  experimenter  manipulates  the  initial  perceptions  of 
each  decision  maker,  the  parameters  of  his  selecti\’e  perception  biases, 
his  distortions,  and  his  forgetting  mechanisms.  He  can  e'xamine  the 
cognitive  structure  and  attention  spaces  of  the  decision  makers  and  displa 
the  spaces  in  a  variety  of  ways.  Allen  Kessler  is  the  main  programmer 
of  this  model. 
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Time-sharing  vs.  Ordinary  Computer  Usage  -  Elizabeth  J.  Campbell 

To  compile  actual  statistics  of  CTSS  compared  with  the  usual  method 
of  batch  processing  is  difficult,  because  the  approach  to  the  program 
is  usually  different.  On  the  console,  a  programmer  tends  to  be  more 
adventuresome  and  try  new  approaches  to  a  problem.  The  long  turn¬ 
around  time  of  batch  processing  causes  a  programmer  to  write  a  safe, 
conservative  program.  Comparison  is  also  difficult  because  when  a 
programmer  uses  batch  processing,  several  programs  are  going  at  a 
time;  while  in  CTSS,  it  is  possible  to  concentrate  on  one  program  and 
save  time,  since  a  train  of  thought  is  not  interrupted.  The  following 
examples  illustrate  this  difference. 

In  one  afternoon,  an  inexperienced  programmer  and  I  worked  out  a 
routine  for  Q^,  the  Legendre  polynomial  of  the  2*^^  kind,  to  replace  an 
inaccurate  section  of  a  program  sent  to  us  from  Berkeley.  We  tried 
three  new  methods,  comparing  accuracy  against  tabulated  results. 

This  work  would  have  taken  us  two  to  three  days  with  batch-processing 
procedures. 

In  another  instance,  I  spent  about  three  weeks  trying  to  build  a  Bound- 
State  program  from  pieces  of  an  old  program.  I  decided  not  to  do  it  on 
time-sharing  because  there  would  be  many  changes  due  to  the  differences 
between  FORTRAN  and  MADTRAN.  However,  when  the  delay  proved  too 
great,  I  managed  to  convert  the  whole  program  to  CTSS  in  a  weekend.  I 
decided  at  the  time  that  I  could  have  done  the  whole  job  on  CTSS  in  a 
third  to  a  half  the  batch-processing  time. 

Electron  Scattering  Spectra  -  Elaine  Miller 

I  have  been  working  on  a  set  of  programs  that  involve  problems  associ¬ 
ated  with  the  radiative  degradation  of  electron  scattering  spectra.  One 
of  these  programs  is  used  as  a  source  for  predictions  of  data.  These 
data  are  employed  in  a  program  that  introduces  all  radiative  effects  into 
otherwise  undegraded  spectra.  The  program,  when  used  as  the  source 
of  predictions,  computes  elastic  and  inelastic  yields  throughout  triangles 
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selected  according  to  the  most  reliable  predictions.  Doth  straight  pole 
fit  and  Yang  assumption  are  available  for  the  computation  of  nucleon 
form  factors . 

Since,  when  primary  electron  energy  is  increased  to  20  BeV,  the  yield 
often  exists  only  within  a  very  small  portion  of  the  selected  triangle, 
this  small  portion  is  to  be  expanded  into  a  larger  area.  Using  time¬ 
sharing,  I  tested  various  methods  of  decreasing  mesh  size,  thereby 
increasing  area.  Through  comparison,  elimination,  and  merging  of 
methods,  a  new,  efficient  method  evolved.  For  each  set  of  data,  time¬ 
sharing  will  be  used  to  determine  a  now  triangle,  that  strictly  adheres 
to  the  conditions  imposed  upon  the  original  triangle,  but  uses  a  finer 
mesh  to  obtain  additional  information, 

Bound-State  Wave  Functions  -  Elmer  Bartels 

I  have  revised  and  perfected  a  program,  which  Elizabeth  Campbell 
started  with  the  Bound-State  calculation,  as  the  first  of  a  package  of 
programs  for  the  use  of  physicists.  A  question-and-answer  set  of 
routines  compute  the  normalized  Bound-State  wave  function  and  binding 
energy  for  a  Wood-Saxon  well.  The  user  may  LOADGO  the  proper  set 
of  files  and  then  type  in  well  parameters  and  eigenvalue  quantun,  numbers 
as  they  are  requested.  He  has  the  option  to  guess  the  binding  energy,  or 
the  program  will  supply  the  guess  and  subsequently  compute  the  wave 
function  and  eigenvalue.  Anotlrer  option  in  the  program  stores  from  two 
to  four  Bound-State  wave  functions  for  computing  the  radial  integral 

00 
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Tlic  user  may  either  compute  another  set  of  wave  functions  and  binding 
energies  with  the  same  well  parameters,  or  he  may  return  to  the  begin¬ 
ning  of  the  program  and  insert  all  new  parameters  of  the  Wood-Saxon 
well  without  having  to  LOADGO  the  files  again. 
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Nucleon-Nucleon  Interactions  -  Nancy  Spencer 

This  use  of  CTSS  is  to  investigate  the  applicability  of  the  boundary  condition 
model,  developed  by  Professor  Feshback  and  Dr.  Lomon,  to  the  proton- 
proton  and  neutron-proton  scattering  data  in  the  energy  range  extending 
to  350  MeV.  A  basic  set  of  parameters,  necessary  to  specify  the 
problem,  are  read  into  a  program,  which  calculates  such  data  as  scatter¬ 
ing  length,  partial  cross  section,  total  cross  section,  and  polarization. 

I  then  compare  these  theoretical  results,  predicted  by  the  model,  to 
experimental  data  at  various  energies  and  apply  a  chi-square  criterion 
to  determine  goodness -of-fit.  After  the  program  is  run  off-line,  each 
parameter  is  varied  systematically  by  a  specified  amount  until  an  optimum 
set  of  parameters  appears. 

Since  only  one  parameter  varies  at  a  time,  this  process  does  not  take 
into  account  any  correlation  which  may  exist  between  the  parameters. 

Thus,  with  CTSS,  a  physicist  can  examine  the  results  as  they  are  com¬ 
puted  and  can  interrupt  the  program  and  change  a  parameter  whenever 
he  sees  a  correlation  or  feels  that  one  may  exist.  With  such  an  inter¬ 
action  between  a  physicist  and  a  computer,  the  development  model, 
which  gives  a  precision  fit  to  the  data,  can  proceed  rapidly. 
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Magnet  Design  Parameters  -  James  E.  Spencer 

The  main  purpose  of  this  problem  is  to  produce  a  computer  program 
capable  of  calculating  all  first  and  second-order  focusing  coefficients 
for  a  serial  array  of  magnets.  The  design  parameters  for  these 
magnets  are  arbitrary,  so  that  any  geometrical  form  of  boundary  can 
be  studied.  Functional  forms  of  greater  than  second  order  (for  the 
geometrical  form  of  boundaries)  will  not  affect  the  first  and  second- 
order  focusing  coefficients,  but  they  are  included  so  the  program  can 
be  used  as  a  general  ray-tracing  subroutine.  A  realistic  field  descrip¬ 
tion  is  being  used,  including  fringing  fields  at  entrance  and  exit  bound¬ 
aries  which  take  into  consideration  the  effects  of  finite  u'idths  and 
curvatures  on  the  iso -induction  lines. 

Analytic  expressions  for  the  fringing  fields  in  the  median  plane  are 
obtained  from  experimental  data  by  curve  fitting.  A  knowledge  of  the 
field  off  the  ntedian  plane  is  then  obtained  using  Taylor  e.xpansions 
through  the  third  order.  The  resulting  expressions  for  the  field  are 
then  used  with  equations  of  motion  to  trace  rays  of  any  initial  conditions 
through  the  system. 

Once  the  performance  of  this  program  has  been  determined  by  making 
specific  calculations  on  known  cases,  other  cases  of  interest  will  be 
investigated  with  a  view  towards  learning  the  utility  of  these  coefficients 
in  the  actual  design  of  magneta. 
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Mathematical  Research  and  CTSS  -  Louis  N.  Hovvard 

f 

This  project  has  been  largely  directed  toward  exploring  some  of  the 
potentialities  of  CTSS  for  assisting  a  mathematician  in  "everyday" 
aspects  of  his  research.  In  general,  most  pure  mathematicians,  and 
many  applied  mathematicians,  almost  never  use  a  computer.  This  is 
probably  because  computers  have  seemed  most  appropriate  for  large 
computations  on  completely-formulated  and  numerically-analyzed  prob¬ 
lems;  and  mathematicians  are  seldom  interested  in  such  problems. 
However,  it  seems  that  there  are  various  special  aspects  of  mathemati¬ 
cal  research  in  which  the  ability  to  rapidly  make  relatively-simple, 
tedious  computations,  not  necessarily  of  a  strictly  numerical  nature, 
could  be  of  considerable  use;  for  example  in  seeking  counter-examples, 
.testing  conjectures,  or  trying  to  obtain  an  overall  numerical  or  geomet¬ 
rical  understanding  of  some  formula  or  relationship.  CTSS,  used  with 
certain  basic  programs  designed  to  be  fairly  self-explanatory,  seems 
to  offer  the  possibility  of  providing  such  service  to  mathematicians 
without  requiring  much  knowledge  of  programming. 

We  have  experimented  with  several  programs  oriented  toward  this  kind 
of  work.  One  example  is  a  very  flexible  and  almost  self-explanatory 
program  for  constructing  rough  typewriter  graphs  of  contour  lines  for 
arbitrarily-given  real  functions  of  two-variables.  This  is  quite  easy 
for  the  uninitiated  to  use,  and  has  been  found  valuable  in  obtaining  a 
qualitative  "feeling"  for  the  behaviour  of  such  functions.  Another 
example  is  a  program  for  making  algebraic  computations  with  matrices 
whose  entries  are  elements  of  an  arbitrary  finite  field.  The  field  is 
specified  by  giving  its  characteristic  p,  number  of  elements,  and  an 
irreducible  polynomial  over  the  integers  modulo  p.  A  number  of  matrices 
can  be  handled  simultaneously  and  operations  with  them  are  carried  out 
by  giving  various  simple  commands.  A  third  example  is  a  program  for 
constructing  the  Lagrange  interpolation  polynomial  through  an  arbitrary 
set  of  points  and  carrying  out  various  operations  with  this  polynomial  on 
the  basis  of  simple  commands.  Other  similar  programs  have  been 
explored  for  computations  in  group  theory,  differential  equations  (specially 
adapted  to  time -sharing),  and  number  theory.  Finally,  a  program  still 
under  development  is  intended  to  provide  a  simple  and  flexible  means  for 
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carrying  out  operations  and  sequences  of  operations  that  occur  in 
elementary  numerical  analysis.  If  a  suitable  program  can  be  developed, 
we  hope  to  have  students  use  it  for  a  numerical  analysis  class. 

While  the  greatest  use  of  CTSS  by  mathematicians  is  likely  to  be  with 
the  type  of  small-scale  everyday  problems  mentioned  previously,  some 
mathematicians  will  occasionally  be  interested  in  relatively-large-scale 
computational  problems.  Two  such  problems  in  which  time-sharing 
appears  to  be  particularly  useful  have  been  explored,  partially  at 
Project  MAC  and  partially  at  the  Computation  Center.  One  is  an 
elaborate  program  for  studying  the  structure  of  finite  groups,  developed 
by  Dr.  C.  C.  Sims;  the  other,  a  set  of  programs  for  the  study  of  prob¬ 
lems  in  hydrodynamic  stability,  developed  by  L.  N.  Howard  and  Prof. 
M.  Landahl  of  the  Aeronautics  Department.  In  both  of  these  cases,  the 
possibility  of  guiding  computation  on  the  basis  of  previous  results  has 
proved  to  be  very  useful. 

In  general,  these  experiments  suggest  that  a  time-sharing  system  can 
be  effectively  and  usefully  employed  in  some  parts  of  mathematical 
research;  though  for  general  applicability,  further  simplification  of 
programming,  and  simplifying  at  least  a  part  of  the  command  system  , 
of  CTSS,  should  be  brought  about.  The  principal  limitation  of  the 
specialized  programs  mentioned  previously  is  that,  in  the  case  of 
complicated,  flexible  programs  which  actually  do  relatively  little  com¬ 
putation,  like  the  contour-line  program  and  the  numerical-analysis 
program,  the  amount  of  time  lost  in  swapping  seems  excessive. 
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On-Line  Simulation  System 
The  OPS-3  System 

Computer  Techniques  in  Business  Problems 
Console-Operated  Statistical  Routines 
An  Automated  Slock  Exchange 
Computer-Planning  System 
Industrial  Dynamics 
On-Line  Array  Processor 
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On-Line  Simulation  System  -  Martin  Greenberger  and  Malcolm  M.  Jones 

The  availability  of  time -sharing  systems  makes  development  of  an  on-line 
simulator  both  feasible  and  desirable.  From  our  experience  with  GPSS, 
SIMSCRIPT,  and  a  preliminary;  on-line  simulator  called  OPSIM,  we  have 
drawn  up  specifications  for  a  flexible,  interactive  simulation  system.  Our 
on-line  simulator  uses  OPS,  a‘ multi-purpose  system  with  compound 
operators  (KOP).  One  KOP  may  call  another  and  transmit  parameters 
by  way  of  a  calling  sequence;  one  may  also  activate  or  schedule  the  occur¬ 
rence  of  another  by  means  of  special  schedule  operators.  A  special  com¬ 
pound  operator,  called  AGENDA,  contains  a  list  of  future  activities,  which 
are  automatically  removed  from  AGENDA  after  it  has  been  activated.  A 
user  may  at  any  time  inspect  and  modify  AGENDA. 

The  execution  of  a  compound  operator  may  result  in  the  alteration  of  state 
\ariables  stored  in  the  common  system,  or  it  may  result  in  changes  to 
AGENDA,  Among  the  operators  that  change  AGENDA  are  WAIT  and  DELAY 
WAIT  is  conditional  upon  the  value  of  a  state  variable,  and  DELAY  is  un¬ 
conditional  and  specifies  that  a  certain  period  of  time  elapse  before  execu¬ 
tion  of  the  KOP  is  continued.  Thus,  activities  listed  in  the  AGENDA  have 
either  a  time  of  execution  or  a  Boolean  condition  associated  with  them. 

Two  other  operators  allow  activities  scheduled  in  the  AGENDA  to  be  can¬ 
celled  or  rescheduled.  In  addition,  the  LOCAL  operator  lists  variables 
to  be  sa\-ed  when  an  activity  is  interrupted  by  DELAY  or  WAIT. 

Among  the  simulations  in  progress  or  completed  are: 

1.  A  simple  multi-server  queuing  model 

2.  An  inventory  distribution  system  I 

3.  A  process -control  computer  system 

4.  A  three-level  priority  interrupt  scheduling  system 
(similar  to  the  CTSS  scheduling  algorithm) 

5.  A  hospital  recepfion  desk 

6.  An  automated  stock  exchange. 
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The  OPS -3  System  -  Martin  Greenberger,  James  H.  Morris,  Jr.  ,  and 
David  N.  Ness 

OPS-3  is  a  multi-purpose  on-line  system  within  which  a  user  can  write 
and  test  his  own  programs  regardless  of  the  nature  of  his  particular  pro¬ 
ject.  Since  the  system  is  comnletely  open,  a  user  can  add  his  own  sub¬ 
routines  and  use  them  in  conjunction  with  system  subroutines,  as  long  as 
they  obey  customary  conventions  regarding  common  storage.  He  is  free 
to  run  his  programs  within  the  system  during  the  formative  and  testing 
stages,  and  then  run  them  independently  during  the  production  stage. 

The  system  contains  operators  for  symbolic  vector  and  matrix  calcula¬ 
tions,  statistical  analysis,  live  data-base  reorganization,  dynamic  loading 
of  routines,  FORTRAN-like  program  writing  with  instant  execution,  and 
on-line  simulations.  While  sitting  at  a  console,  the  user  may  access 
information  about  the  system.  The  information  is  stored  in  regular  CTSS 
disc  files  and  may  be  retrieved  from  within  the  OPS-3  system  by  a  single 
key-word  mechanism.  OP.S-3  is  now  being  used  in  a  number  of  projects. 
These  include: 

1.  an  on-line  simulation  system  within  which  a  user  may  change 
the  course  of  the  simulation  or  examine  partial  results, 

2.  a  real-time  simulation  of  the  stock  exchange  that  involves  both 
simulated  and  console -initiated  orders, 

3.  a  polynomial  manipulation  system, 

4.  a  system  for  statistical  analysis  in  the  social  sciences, 

5.  a  general  debugging  system  for  MAD  programs  in  which 
programs  are  first  run  interpret! vely  and  then  compiled  into  regular  BSS 
subroutines . 

OPS-3  is  currently  available  as  the  CTEST6  command  at  Project  MAC. 
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Computer  Techniques  in  Business  Problems  -  Donald  C.  Carroll, 

Paul  Clermont,  Thomas  J.  Johnson,  Lalit  Kanodia,  and  Francis  J.  Russo 

This  project  investigates  the  use  of  heuristic  programming  and  interactive 
problem-solving  techniques  in  complex  business  problems,  such  as 
scheduling  and  sequencing  production  in  a  job  shop.  We  developed  a 
simulation  model  of  a  hypothetical  shop,  and,  by  testing  some  heuristics 
procedures  for  the  sequencing  and  routing  of  jobs  under  varied  and  ex¬ 
tensive  conditions,  determined  the  general  steady-state  behavior  of  this 
model.  The  first  part  of  our  research  studied  fully-automatic  decision¬ 
making,  in  which  the  heuristics  for  sequencing  and  routing  require  no 
human  guidance;  the  second  part  is  to  study  the  rule  of  a  human  pattern- 
recognizer  and  decision-maker  who  is  closely  coupled  with  the  automatic 
system.  In  the  latter  research,  we  are  investigating  the  establishment 
of  manning  levels  for  the  shop,  overtime  decisions,  and  load-leveling. 

These  decisions  determine  the  structure  under  which  an  automatic  decision 
system  must  operate;  consequently,  we  call  them  structuring  decisions, 

A  complete  set  of  programs  that  study  non-interactive  decision  making  is 
now  available.  It  includes  a  highly-parametric  input  generator  (MADTRAN), 
a  variety  of  simulation  models  that  encompass  such  features  as  multiple- 
channel  work  stations,  multiple-component  orders,  alternate-routing 
possibilities,  and  machine  substitution  (FAP),  and  finally  an  output  analysis 
program  (MADTRAN).  Three  pieces  of  research  are  complete.  One,  a 
study  of  heuristic  sequencing  rules  for  single  and  multiple  components 
orders,  demonstrated  the  general  superiority  of  the  so-called  COVERT 
rule  for  sequencing.  Another,  which  dealt  with  heuristics  for  selecting 
alternate  routes  for  jobs,  produced  a  superior  heuristic  derived  from 
COVERT.  The  last  study,  on  the  effect  of  feed-back  phenomena  in  heavily- 
loaded  shops,  led  to  several  generalizations  about  the  control  of  inputs 
and  processing  rates  to  obtain  stable  operation.  A  by-product  of  these 
studies  of  the  minutiae  of  the  decision  system  is  a  comprehensive  under¬ 
standing  of  the  behavior  of  these  models  with  respect  to  their  aggregate 
behavior. 

An  extension  of  these  studies  to  heuristics  for  machine  substitution  is 
almost  complete;  however,  the  major  effort  is  in  interactive 
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decision-making.  Toward  this  end,  we  are  elaborating  some  of  the  basic 
programs  to  include  such  structural  parameters  as  manning  levels  by 
shift,  and  we  are  developing  generalized  information  retrieval  and  display 
programs  for  the  ESL  Display  Console.  A  rudimentary  version  of  our  inter 
active  system  should  be  available  in  the  fall. 

Console-Operated  Statistical  Routines  -  James  R.  Miller  and 
Stephen  Whitelaw 

A  battery  of  statistical  routines  have  been  written,  tested,  and  debugged 
for  use  with  a  CTSS  console,  primarily  for  the  testing  of  hunches  and 
limited  hypotheses.  Data  are  typed  directly  into  the  console,  analyses 
are  performed,  and  results  are  printed  out  immediately.  Many  of  our 
routines  are  non-parametric  in  nature,  and  are  oriented  toward  the 
analysis  of  typical  problems  in  the  social  sciences.  They  include: 

1.  Standard  two-dimensional  contingency  analyses. 

2.  Special  three-dimensional  contingency  analyses. 

3.  Rank-order  intercorrelation  and  partial  correlation  analyses. 

4.  Binomial  tests. 

5.  Single  and  multiple  percentage  difference  tests. 

6.  Homogeneity  test  on  samples  classified  strictly  according  to 
nomial  categories. 

7.  The  Mann-Whitney  U  test  of  the  difference  of  the  sample  medians 
We  have  also  written  several  parametric  routines: 

r 

1.  Linear  intercor relation  and  partial  correlation  analyses. 

2.  Simple  and  multiple  regression  analyses. 

3.  A  T-test  of  the  difference  between  two  sample  means. 

4.  A  one-way  analysis  of  variance. 

5.  A  test  for  normality  of  sample  data. 

We  shall  expand  the  range  and  capability  of  our  current  library  by  improv¬ 
ing  existing  routines  and  by  adding  new  ones.  When  S.  Whitelaw  has 
adapted  the  routines  to  tlie  OPS  system,  a  user  will  be  able  to  call  routines 
by  name  from  the  console,  with  the  name  of  the  data  matri.x  as  argument, 
to  compound  routines,  and  to  save  data  and  results  for  future  use. 
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An  Automated  Stock  Exchange  -  Stephen  Whitelaw 

With  the  OPS  system,  methods  of  automating  large  securities  markets, 
such  as  the  New  York  Stock  Exchange  have  been  explored.  This  research 
has  included  the  programming  of  subroutines  or  operators  to  automatically 
handle  all  functions  of  a  large  stock  exchange;  entry  of  orders,  consumma¬ 
tion  , of  orders,  modification  of  completed  transactions,  and  automatic 
bookkeeping.  Additional  operators  embody  the  decision  rules  of  specialists 
who  mainta.in  a  market  by  committing  capital. 

Several  sets  of  decision  rules  have  been  tested  on  different  simulated 
markets.  The  performance  of  these  rules  is  measured  by  price  and  net 
return  criteria.  The  market  model  accepts  simulated  orders  from  an 
order  generator.  By  treating  a  common  file  as  an  intermediate  storage 
area  for  incoming  orders,  reports  of  completed  transactions,  requests 
for  quotations,  and  output  of  information,  it  can  also  accept  orders 
entered  from  other  consoles.  The  system  gives  faster  response  times 
than  those  found  in  non-automated  markets,  and  a  programmed  specialist 
is  able  to  control  the  market  under  very  adverse  conditions. 

Currently  the  system  is  being  reprogrammed  to  function  in  two  modes: 
mode  1  provides  for  continuous  simulation  and  output  of  consummated 
transactions,  and  mode  2  for  the  entry  of  console  orders  into  the  simulation. 
Since  the  console  remains  undisturbed  by  this  entry,  it  can  function  both 
as  a  ticker  tap that  reports  the  results  of  the  simulated  market,  and  also 
as  a  broker's  device  that  enters  orders  into  the  same  market.  This 
functioning  should  make  an  effective  management  game  and  extend  the 
use  of  the  automated  exchange  as  a  research  tool.  The  simulated  market 
will  be  synchronized  with  the  computer  clock  to  prevent  the  simulation 
from  getting  ahead  of  real  time. 
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Computer-Planning  System  -  James  C.  Emery 

This  project  studies  a  management  computer  planning  system.  Because 
planning  requires  a  system  that  can  be  augmented  and  modified  extensively, 
careful  attention  has  been  paid  to  general  techniques  that  provide  for  such 
evolution.  With  such  a  model,  a  planner  will  be  able  to  make  sequential, 
hierarchical  search  for  improved  plans.  First,  from  an  aggregate  model 
he  will  construct  a  large-scale  plan,  which  is  generated  through  a  sequential 
search  of  alt  rnatives.  In  this  search,  the  human  planner  proposes  alter¬ 
natives  and  the  computer  transforms  a  proposed  plan  into  its  predicted 
consequences.  Each  probe  for  an  improved  plan  provides  the  planner  with 
additional  information  for  the  next  probe.  If  the  process  converges,  the 
aggregate  plan  emerging  from  it  is  then  amplified  by  more  detailed  models 
within  the  constraints  provided  by  large-scale  plans.  This  hierarchical 
process  continues  until  plans  are  sufficiently  detailed  for  execution. 

The  system,  specifically  designed  for  budget  planning,  accepts  a  fairly 
wide  class  of  functional  relationship  among  variables.  It  will  also  deal 
with  hierarchically-defined  variables,  such  as  manufacturing  cost  defined 
as  the  sum  of  direct  labor,  direct  material,  and  overhead  cost.  Since 
functional  relationships  and  the  hierarchical  definition  of  variables  are 
represented  primarily  in  the  form  of  list  structures,  the  system  has  the 
advantages  of  both  flexibility  and  economy  of  storage.  During  the  past 
year,  I  have  developed  the  general  structure  of  the  system  and  have  also 
begun  to  specify  storage  organization  and  the  algorithms  for  the  retrieval 
and  manipulation  of  stored  data  and  to  define  the  set  of  subroutines  needed 
for  the  generation  of  displays  for  budget  planning. 
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Industrial  Dynamics  -  Jay  W.  Forrester 

The  Industrial  Dynamics  Group  uses  time-sharing  and  the  DYNAMO  Com¬ 
piler  for  analysis  and  design  of  feedback  processes  in  industrial  organiza¬ 
tions.  These  systems  contain  multiple-loop  feedback  structures  with 
highly-nonlinear  relationships.  Sixtieth-order  models  of  systems  have 
been  simulated  which  contain  several  hundred  variables.  Such  mathema¬ 
tically  intractable  systems  contain  important  modes  of  behavior  which 
are  totally  absent  in  the  simpler  linear  systems.  The  DYNAMO  Compiler 
is  very  efficient!  for,  in  only  a  few  seconds  of  computer  time,  it  can 
simulate  a  system  of  a  hundred  or  more  equations  through  two  hundred 
time  steps.  The  Compiler  contains  a  graphical  plotting  routine  which 
makes  time  plots  of  systems  responses  on  remote  printing  stations. 

Orie  use  of  the  system  has  explored  the  dynamics  of  corporate  growth. 

Many  of  the  important  growth  processes  depend  on  the  interactions  between 
positive  feedback  loops,  negative  feedback  loops,  and  the  coupling  between 
these  created  by  nonlinear  elements.  Another  use  has  developed  illustra¬ 
tive  material  for  a  student  seminar  on  the  principles  underlying  the  behavior 
of  feedback  systems.  These  principles  are  generalizations  about  system 
behavior  which  have  not  previously  been  organized  and  illustrated.  Some 
of  them  are  implicit  in  the  mathematics  of  linear  feedback  systems,  but 
have  not  been  clearly  identified  and  illustrated:  others  exist  only  inmultiple- 
loop  nonlinear  systems  and  could  not  have  been  previously  identified  with¬ 
out  the  extensive  empirical  results  obtained  by  computer  simulation. 

On-line  Array  Processor  -  Mayer  E.  Wantman 

Operators  have  been  developed  on  the  OPS-3  on-line  system  that  provide 
a  simple,  flexible  way  of  manipulating  and  analyzing  numerical  data  within 
the  computer.  They  can  be  used  for  general  processing  arrays,  fitting 
of  regression  models,  transformation  of  variables,  rank-ordering  of  data, 
calculation  of  roots  of  polynomials  of  arbitrary  degree,  error  analysis, 
generalized  desk  calculation,  and  MAD-type  programming.  They  are; 

1.  Assignment  Operator.  The  operator  SET  specifies  and  executes 
MAD-type  assignment  statements.  Operands  may  be  vectors  and  matrices, 
as  well  as  elements.  In  addition,  SET  possesses  embedded  functions. 
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calls,  matrix  multiplication  and  transposition,  a  power  generator,  and  a 

variety  of  logical  operations.  Lines  are  of  the  form 

SET  ALPHA=(BETA  +  GAMMA)=;=DELTA.  P.  (A  +  B) 

2.  Regression  Operator.  The  operator  FIT  creates  a  multiple 
linear  or  simple  polynomial  model  derived  from  given  arguments.  Sampl 
points  may  be  weighted  and  residuals  saved.  As  an  example,  multiple 
linear  regression  of  a  dependent  variable  (stored  in  the  vector  named 
ALPHA)  on  independent  variables  (stored  in  the  vectors  named  BETA, 
GAMMA,  and  DELTA)  would  be  performed  by  the  line 

FIT  ALPHA  TO  BETA  GAMMA  DELTA. 

A  cubic  fit  of  ALPHA  to  BETA  would  be  specified  by 
FIT  ALPHA  TO  BETA  DEGREE  3. 

Rp°t-Solving  Operator.  The  operator  ROOT  uses  the  Newton- 
Raphson  method  of  successive  approximations  for  the  determination  of  the 
root  of  a  polynomial  equation  in  one  variable.  Since  the  user  specifies 

the  starting  value,  he  can  obtain  all  roots  by  setting  different  initial  values 
Exponents  need  not  be  integral. 


—Bering  Operator.  The  operator  RANK  takes  a  vector  of 
numbers  and  creates  a  second  vector  which  is  an  ordered  version  of  the 
first.  A  vector  which  describes  the  mapping  performed  may  also  be 


calculated. 


5.  ^rting  Operator.  The  operator  SORT  sorts  by  intervals  on  a 
specified  vector.  The  intervals  may  be  specified  either  e.xplicitly  or  by 
range  and  spacing.  If  any  samples  fall  outside  the  sorting  range,  an 
appropriate  message  is  printed. 


APPENDICES 


Appendix  A 

Project  MAC  Memoranda 

Appendix  B 
M.I.  T.  Thesis 

Appendix  C 
External  Publication 

Appendix  D 

Project  MAC  Technical  Reports 


APPENDIX  A 


PROJECT  MAC  MEMORANDA 


MEMORANDUM 
MAC- M/ No. 

SUBJECT 

AUTHOR 

DATE 

172 

On-Line  Documentation  of  the 
Compatible  Time-Sharing  System 

J,  Winett 

7/21/64 

173 

SQZ/PADBIN  (Ref:  MAC-M-157) 

R,  Hussey 

7/20/64 

174 

Abstracts  of  CTSS  Console  Com¬ 
mands  (CC239)  (174-2.  1/4/65; 
174-3,  3/7/65 

E.  Kliman 

8/  3/64 

175 

Access  to  the  MAC  System  via 
Alternate- Use- TWX  ("TWXyPrime 

R.  G.  Mills 

") 

7/24/64 

176 

String  Manipulation  in  the  New 
Language  (AIP  memo  71) 

D.  Bobrova 

7/28/64 

177 

CTSS  Library  Abstracts  (CC240) 
(177-2,  11/23/64;  177-3,  3/22/65 

E.  Kliman 

8/6/64 

178 

Crunch  Command  in  CTSS 

R.  Bayles 

8/7/64 

179 

BEFAP  Command  within  CTSS 

R,  Bayles 

8/  10/ 64 

180 

The  ARCHIV  Command 

A.  Scherr 

10/  19/64 

181 

The  Structure  of  On-Line  Infor¬ 
mation  Processing  Systems 

E.  Glaser 

J.  Dennis 

10/1/64 

182 

Preliminary  Notes  on  the  Hard 
Matter  for  the  MAC  635 

E.  Glaser 

8/  28/ 64 

183 

Proposed  Instructions  on  the 

GE  635  for  List  Processing  and 
Push  Down  Stacks  (AIP  memo  72) 

M.  Levin 

9/18/64 

184 

On  Memory  Addressing  in  the 

GE  635 

M.  Minsky 
M.  Levin 

9/22/64 

185 

MAC  and  Its  Users 

U.  Neisser 

9/29/64 

186 

A  Note  on  Asynchronous  Parallel 
Processing  (MSG  memo  5) 

H.  Witsenhausen 

7/64 

187 

Nesting  and  Recursion  of  Proce¬ 
dure  in  Segmented  Memory 
(MSG  memo  6) 

J.  Dennis  10/28/64 

E.  C.  Van  Horn 

188 

Automatic  Scheduling  of  Priority 

J.  Dennis 

10/  64 

Processes 


179 


memorandum 

MAC- M- No.  SUBJECT 


189 


n  Example  of  Intersphere  Com¬ 
munication  and  Asynchronous 
Parallel  Processing:-  Typewritei 
Console  Message  Handling  by 
Protected  Service  Routines 
(MSG  memo  8) 


author 

J.  Dennis 


DATE 

3l2.blb4: 


190 

191 

192 

193-2 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 


General  Comments  about  PDP-6 
Coding  and  Hardware  (PDP-6 
memo  1) 

D.  Edwards  9/25/64 

lECO  6  (PDP-6  memo  2) 

D.  Edwards 

Unrecognizable  Sets  of  Numbers 
(AIP  memo  73) 

TYPSET  and  RUNOFF,  Memor¬ 
andum  Editor  and  Type- out 
Commands  (CC  244-2) 

M.  Minsky  11/2/64 
S.  Papert 

J.  Saltzer  11/6/64 

On  the  Control  Specifications  of 
Computation 

D.  Kuck 

11/64 

p,  A  Context  Editor  for  Card 
3/15/65) 

(CC245;  CC234-1) 

R.  Dailey 

11/20/64 

,  Some  Aspects  of  the  Structure  of 
Computation 

D.  Kuck 

R.  Kain 

11/13/64 

Some  Time-Sharing  Ideas 

Relocation  in 

the  AED-O  System  (9442-M-llO) 

13.  Kuck 

T.  Fox 

11/13/64 

10/26/64 

Facilities  of 

the  AED-O  Language  (9442-M-112) 

R.  Coe 

11/18/64 

Processes,  Spheres  of  Protection 
(MSg"9)"'’""''""‘ 

E,  Van  Horr 

(no  date) 

New  Oprating  System  for  the 
isSE  Display  Console  (9443-M-118 
(Cont'd  12/10/64) 

R.  Bayles 

11/23/64 

Proposal  to  Improve  the  Rotation 

(^4lrM/u4V'®'^ 

R.  Stotz 

11/27/64 

A  Compacting  Free  Storage  Pro¬ 
gram  (9442-M-113 

B.  Wolman 

3/5/65 

AED-O  Compiler  for  Batch  P-o- 
cessing  (9442-M-115 

H.  Spencer 

12/2/64 

ISO 


memorandum 


MAC-M-No. 

SUBJECT 

AUTHOR  DATE 

205 

MADBUG;  a  mad  Debugging 
System  (CC247) 

R.  Fabry  11/25/64 

206 

CTSS  LISP  Notice:  Supplement  to  T.  Hart  12/7/64 

A.  I.  Memo  No.  67  (AIP  memo  74)  /  U  b4 

207 

Internal  Memos  for  AED  Users 
(9442-M-116) 

C.  Feldmann  12/8/64 

208 

AED  Flash  10,  New  CTEST2 
Command  {9442-M-I19) 

C.  Feldmann  12/ 10/ 64 

209 

The  Naval  Mine  Warfare  Game 

I.  Sutherland  12/18/64 

210 

Graphical  Output  Device 

A.  Rutchka  12/21/64 

211 

New  TWX-Prime  Lines 

R.G.  Mills  12/23/64 

2 1  Z* 

Translation  of  Communication 
with  Automata  bv  Carl  Ariam 

O.  Selfridge  6/10/65 

213 

AED  FLASH  11;  AED  BUG  Usage 
(9442-M-12n  ^ 

T.  Fox  12/28/t)4 

215 

Television  Camera-to-Computer 
Adapter:  PDP- 6  Device  770 
(AIP  memo  75) 

M.  Minsky  1/9/65 

216 

New  B-Core  System  for  Program- 
tning  the  ESL  Display  Console 
(9442-M-  122) 

C.  Lang  4/30/65 

217 

Operating  Manual  for  the  ESL 
Display  Console  (9442-M- 129) 

R.  Stotz  3/9/65 

J.  Ward 

218 

Social  Science  Data  Bank 
Development 

I.  D.  Poole  1/21/65 

219 

The  COMIT  Feature  in  LISP  II 
(AIP  memo  76)  (LISP  II  Proiect 
Memo  2) 

D.  Bobrow  2/18/65 

220 

The  COGO  Language 

D,  Rods  2/25/65 

221 

The  INFO  Command 

J.  Wi'nett  2/25/65 

222 

AED  f  LASH  12:  Multi- Entry 
Procedures  (9442-M- 124) 

D.  Ross  2/ 25 '•65 

224 

Plotting  Subroutine, 

X  Y  jp  JLXJ  r 

T.  Stockham.  3/3/65 

Jr. 

In  older  listing's  of  MAC- M- memos 


this  memo 


was  numbered  as  MAC-M-214. 


MEMORANDUM 


MAC-M-No. 

SUBJECT 

AUTHOR 

DATE 

?.25 

AED  FLASH  13:  Argument  Check¬ 
ing  Procedures  for  AED  (9442-M- 
125)  (225-2,  4/  15/65) 

J.  Walsh 

3/5/65 

226 

AED  FLASH  14:  Availability  of 
AED  JR  Systems  (9442-M-126) 

D.  Ross 

3/9/65 

111 

The  Editing  Routine 

D,  Griffel 

3/8/65 

228 

A  General  Cross- Tabulation 

System  (COM/COM  Sim.  memo  23) 
(See  Additional  Information  on 
Cross- Tabs:  New  Output) 

W,  Selles 

4/24/65 

229 

Proposed  Datanet- 30  Program 

S,  Dunten 

3/18/65 

230 

Mind,  Matter  and  Models 
(AIP  memo  77) 

M,  Minsky 

3/26/65 

231 

Description  of  SNOBOL  Commands 
(CC235-2) 

L,  Pouzin 

4/9/65 

232 

Preliminary  Description  of 
DAEMON,  the  Disk  Dump  and  Re¬ 
load  Program  Package,  and  its 
Operation  (CC  252) 

M,  Bailey 

4/26/65 

233 

Trial  Launching  of  the  OPS- 3 

System 

M,  Jones 

D,  Ness 

4/26/65 

234 

Users'  Manual  on  the  Use  of 
Structural  Language 

R.,  Logcher 

9/ 15/65 

235 

BLODI  Command 

O,  Wright, 
Jr, 

5/4/65 

236 

AED  FLASH  15:  Execution  Time 
Procedure  Calls  (9442-M-131) 

J.  Walsh 

4/22/65 

237 

A  Generalized  File  Structure  and 
Input/ Output  System  (CC  241-1) 

R.  C.  Daley  5/7/65 

R.  J.  Creasy 

R.  M.  Graham 

238 

Some  Facilities  of  the  MAC 
Document  Room 

Document 

Room 

(A.  Bowen) 

5/10/65 

239 

KREADQ,  SETRDQ,  MODRDQ, 

A  Set  of  Format-Free  Input 
Subroutines 

D.  Ness 

6/ 1  6/ 65 

240 

Topics  in  Model  Theory 

M.  Levin 

5/25/65 

241 

PDP-6  LISP  Input- Output  for  the 
Dataphone  (AIP  memo  79) 

W.  Martin 

6/1/65 
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MEMORANDUM 


MAC-M-No. 

SU3  JECT 

AUTHOR  DATE 

242 

PDF- 6  LISP  Input- Output  for  the 
Display  (AlP  memo  80) 

W.  Martin  6/28/65 

243 

Specifications  for  a  Dataphone- 
Driven  Remote  Display  Console 
for  Project  MAC 

R.  H.  Stotz  6/24/65 

J.  E.  Ward 

U.  F.  Gronemann 

244 

Proposal  tor  Commands  for 
Dynamic  Resource  Allocation 
(CC  254) 

G.  Schroeder  4/17/65 

245 

DIRECT;  A  Command  Subroutine 
for  Controlling  Programs 

T.  Crystal  6/23/65 

E.  River 

246 

Computer  Experiments  in  Finite 
Algebra 

W,  D.  Maurer 

6/14/65 

247 

AED  FLASH  16:  Use  of  AED  JR  to 
Parse  FORTRAN  11  {9442-M-134) 

H,  Spencer  (//25/65 

248 

Use  of  MAC DMP  (AIP  memo  83) 

P.  Samson  7/ 1/65 

249 

MAC  PDP-6  DECtape  FHe 

Structure  (AIP  memo  82) 

P.  Samson  7/  1  '  65 

■  250 

PDP-6  TECO,  July  1965  Version 

P.  Samson  7/23/65 

251 

AED  FLASH  17:  Features  of  the 
Revised  AED-JR  {9442-M-136) 

J.  F.  Walsh  7/1/65 

R,  Feldmann 

252 

AED  FLASH  18:  A  Subroutine  for 
Generating  Coons'  Surfaces 
(9442-M-137) 

C.  Lang  6/ 30/ 65 

253 

FMS  Tape  Distribution 
{9442-M-127-3) 

R.  Feldmann  7/ 1  /  65 

D.  E,  Walker 

254 

Derivation  of  a  Mean  Pseudo- 
Associative  Search  Time 
(MSG  memo  1  3) 

E.  C.  Van  Horn 

8/16/65 

255 

Surfaces  for  Computer-Aided 

Design  of  Space  Figures 
(9442-M-139) 

S.  A.  Coons  7/21/65 

256 

Three- Dimension  Pseudo  Pen 
Sub-routine  for  use  with  the  ESL 
Display  Console  (9442-M-138) 

R.  Polansky 

7/13/65 

257 

Syntax  and  Display  of' Mathemati¬ 
cal  Expressions  (XiP  memo  85) 

W.  Martin  7/29/65 

1 

1 


MEMORANDUM 

MAC-M-No. 

SUBJECT 

AUTHOR 

DATE 

.  .  .  REVISIONS  .  .  . 

168-2 

A  Subroutine  Trace  Program 
fo'r'tITSS 

B.  Wolman 

3/5/65 

174-3 

CTSS  Console  Commands 

E.  Kliman 

3/7/65 

177-3 

CTSS  Library  Subroutines 

E.  Kliman 

3/22/65 

193-2 

TYPSET  and  RUNOFF,  Memo¬ 
randum  Editor  and  Type-out 
Commands  (CC-244-2) 

J.  Saltzer 

1/11/65 

195-1 

ED,  a  Context  Editor  for  FAP 
and  MADTRN  Program 

R,  Dailey 

3/15/65 

I 


I 


4 
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APPENDIX  B 

M.  I.  T.  THESIS 

Amstutz,  A,  E.  ,  A  Management  Oriented  Behavox  al  Theory  of  Inter- 

actions  within  Consumer  Product  Markets,  M.  I.  T.  Sloan  School  of 
Management,  Ph.  D.  Thesis,  June  1965 

Armstrong,  A.  E,  ,  A  Braille  Telecommunication  Terminal,  M.  I.  T. 

Department  of  Mechanical  Engineering,  S,  M.  Thesis,  June  1965 

Beckreck,  L.  N.  ,  SEPOL:  A  Soil  Engineering  Problem-Oriented 

Language,  Department  of  Civil  Engineering,  M.  S.  Thesis,  June  1965 

Beusch,  J.  U.  ,  Dynamic  Behavior  and  Control  of  Communication  Networks, 
M.  I,  T.  ,  Department  of  Electrical  Engineering,  Ph.  D,  Thesis, 

May  1965 

Brach,  J,  R,  ,  A  Problem-Oriented  Language  for  Project  Management, 

M.  I.  T.  ,  Department  of  Civil  Engineering,  S.  M.  Thesis,  January  1965 

Carroll,  D.  C,  ,  Heuristic  Sequencing  of  Single  and  Multiple  Compound 

Jobs,  M.I.  T.  ,  Sloan  School  of  Management,  Ph.  D.  Thesis,  June  1965 

Chang,  W.  Y.  ,  Instability  Analysis  of  Transverse  Waves  in  Solid-State 

Plasmas ,  M.  I.  T.  ,  Department  of  Electrical  Engineering,  B.  S.  Thesis, 
June  1965 

Chiodi,  W.  ,  A  Mechanism  for  the  Transformation  of  S,  M.  I.  T.  ,  Depart¬ 
ment  of  Electrical  Engineering,  S.,M.  Thesis,  June  1964 

Diamond,  D.  S.  ,  A  Step  toward  Coniputer  Aided  Preparation  of  Magazine 

Advertising  Copy,  M.  I.  T.  Sloan  School  of  Management,  B.  S.  Thesis, 
June  1965 

Dvorak,  A.  Jr.,  An  Input-Output  Program  for  Electronic  Circuits  Using 
a  CRT,  M.  I.  T.  ,  Department  of  Electrical  Engineering,  Thesis, 

June  1965 
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Efimba,  R.  E,  ,  Analysis  of  a  Bridge  Deck  as  a  Plane  Grid.  M.  I.  T.  , 
Department  of  Civil  Engineering,  M.  S,  Thesis,  June  1965 

Ell'S,  E.  ,  A  MAD-based  Programming  Eanguage.  M.  I.  T,  ,  Depart¬ 
ment  of  Electripal  Engineering,  B.S.  Thesis,  June  1965 

•  -An  Information  Retrieval  System  for  Reference  Filing.  M.  I.  T. 
Department  of  Electrical  Engineering,  B.S.  Thesis,  June  1965 

Fmkelstem,  S.  C.  ,  An  Analysis  of  a  Space  Beam.  M.  I.  T.  ,  Department 
of  Civil  Engineering,  M.  S.  Thesis,  June  1964 

Fiske.  R.  S.  ,  A  Real-Time  System  for  Computer  Monitoring  of  the  Stock 
Exchange,  M.  I,  T.  ,  Sloan  School  of  Management,  B.  S.  Thesis, 

June  1965 

Foster,  A,  T.  ,  Special  Input-Output  Facilities  for  the  IBM  1620, 

M.  I.  r.  ,  Department  of  Civil  Engineering,  M.  S.  Thesis,  May  1964 

Goeke,  R,  F,  ,  jntimduction  to  Machine  Language  Programming  on  a  Time 
Sha_ij3_d_IBM  7094,  M,  I.  T.  ,  Department  of  Electrical  Engineering, 
B.S.  Thesis,  June  1965 


Gonzalez-Zubieta,  R.  ,  On  Some  Aspects  of  Integer  Einear  Programming, 
Pii.  D.  Thesis,  Sloan  School  of  Management,  June  1965  (Also 
released  by  Operations  Research  Center  as  Technical  Report 
No,  16,  June  1965) 


Goodman,  R.  V.,  A  Preliminary  Method  of 
Structural  Design  Language.  M.  I.  T.  , 
S.M.  Thesis,  June  1965 


Analysis  for  a  Computer-Aided 
Department  of  Civil  Engineering, 


Groneman,  U.  F.  ,  Coding  Color  Pictures.  M.  I.  T.  ,  Department  of  Elec¬ 
trical  Engineering,  Ph.  D.  Thesis,  February  1964 

Guttmann.  E,  G.  ,  hvestigation  of  Threshold  Logic  in  High-Speed  Display 

Systems,  M.  I.  T.  ,  Department  of  Electrical  Engineering.  S.M.  Thesis, 
June  1965 
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\ 
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Guzman  Barron-TorreS ,  G.  ,  Optimum  Design  of  Prestressed  Concrete  Highv^^ 
Bridges  using  Digital  Computers,  M.  1.  T.  ,  Department  of  Civil 
Engineering,  S.  M.  Thesis,  May  1965 

Hadley,  R.  A.  ,  Investigation  of  Automated  Direct  Design  of  Rigid  Steel 
Planar  Frames,  M.  I.  T.  ,  Department  of  Civil  Engineering,  M.  S. 

Thesis,  February  1965 

Haines,  E.  C.  Jr.,  The  TREET  List- Processing  Language,  M.  I.  T.  , 

Department  of  Electrical  Engineering,  M,  S.  Thesis,  September  1964 
(Also  produced  by  Mitre  Corporation  as  Report  SR- 133,  April  1965) 

Hamann,  W.  C.  ,  Computer  Aided  Design  of  Slender  Structural  Members, 

M.  I.  T.  ,  Department  of  Mechanical  Engineering,  M,  S,  Thesis, 

May  1964 

Hamilton,  M.  L.  ,  and  A.  D.  Weiss,  An  Approach  to  Computer- Aid^ 

Preliminary  Design,  M.  I.  T.  ,  Department  of  Mechanical  Engineering 
and  Naval  Architecture,  S,  M.  Thesis,  August  1964  (Also  as 
ESL-TM-228,  January  1965) 

Hershdorfer,  A.  M.  ,  Optimal  Routing  of  Urban  Traffic,  M.  I.  1 .  ,  Depart¬ 
ment  of  Civil  Engineering,  Ph.  D.  Thesis,  June  1965 

Hoag,  J.  D.  ,  Instabilities  in  Transverse  Waves  Along  P  for  a  Beam  Type 
Distribution,  M.  I.  T.  ,  Department  of  Electrical  Engineering, 

B.S, 'Thesis,  October  1964 

Ishii,  H,  ,  Instabilities  in  InSb,  M.  I.  T.  ,  Department  of  Electrical  Engineer¬ 
ing,  S,  M.  Thesis,  May  1965 

Johnson,  W.  F.  ,  Use  of  Traffic  Volume  Data  in  Evaluation  of  Highway 

User  Costs,  M.  I.  T.  ,  Department  of  Civil  Engineering,  S.  M.  Thesis, 
June  1964 

Kanodia,  L.  S.  ,  Queue  Balancing  -  Another  Approach  to  the  Theory  of 
Schedulmg,  M.  I,  T,  Sloan  School  of  Manag 'ment,  M.  S,  Thesis, 


June  1965 


187 


Kanstroom,  E.  ,  An  Information  Retrieval  Program  Compatible  with  the 
M.  I.  T.  Time-Shared  Computer  System.  M.  I.  T,  Sloan  School  of 
Management,  B.S.  Thesis,  June  1965 
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